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Abstract 

The objective of Phase II was to design, build and demonstrate a three dimensional laser 
fluorescence anemometer for use in the Langley 16- by 24- Inch Water Tunnel. Innovative optical 
design flexibility combined with compact and portable data acquisition and control systems have 
been incorporated into the instrument. This.will allow its use by NASA in other test facilities. The 
final instrument and support systems differ in several significant aspects from the design envisaged 
during our Phase II proposal preparations. Our original mirror traverse alignment concept has 
been replaced by a more versatile fiber optic system. This facilitates normal and off-axis beam 
alignment, removes mirror losses and improves laser safety. This added optical flexibility will also 
enable simple adaptation for use in the adjacent jet facility. New proprietary concepts in 
transmitting color separation, light collection and novel prism separation of the scattered light have 
also been designed and built into the system. Off-axis beam traverse and alignment proved much 
more complex than initially conceived. This led to the requirement for a specialized, programmable 
traverse controller and the inclusion of an additional traverse for the off-axis arm. To meet this 
challenge, an “in-house” prototype unit was designed and built and traverse control software 
developed specifically for the water tunnel traverse applications. A specialized data acquisition 
interface was also required. This was designed and built for the Laser Fluorescence Anemometer 
system. 

Introduction 

At present, significant efforts are being made to effect design changes which will improve 
aircraft agility, maneuverability and performance. But, although significant progress has been 
made in computational aerodynamics, reliable design changes still cannot be made without recourse 
to experiment. Attempts to extend tactical flight envelopes still require extensive preflight ground 
based model testing. Unfortunately, conventional wind tunnel testing is expensive and time 
consuming and most facilities were built before present-day optical methods for quantitative flow 
field measurements were envisaged. Consequently, the non-intrusive detailed documentations of 
lee-side vortex flow- fields which are often required to support design evaluation and optimization 
are few. 

However, in the past, qualitative water tunnel simulations have guided many practical 
designs and, since most of these facilities have been built with excellent optical access, they are 
ideally suited for use in advanced flow field diagnostics. Since the performance of most lifting and 
maneuvering bodies is governed by extensive transitional and turbulent viscous wakes and vortical 
lee-side flows, non-intrusive optical measurement techniques are required. Consequently, water 
tunnels offer the opportunity to obtain inexpensive, detailed flow field measurements to support 
"cut and try" designs and basic research. 

To realize this capability, a two dimensional laser fluorescence anemometer was built and 
tested in the Ames-Dryden Water Tunnel during Phase I. The instrument was used in an 
experimental study of vortex flow fields designed to determine the mechanisms and feasibility of 
controlling vortex breakdown by introducing relatively low rates of jet blowing along the vortex 
core. The objective of Phase II was to build a three dimensional instrument for studies in the 
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Langley 16- by 24- Inch Water Tunnel. 


Background 

When a slender delta wing is at an angle of attack to an oncoming stream, the upper and 
lower surface boundary layers flow outward and separate from the leading edges to form two free 
shear layers that roll up into a pair of vortices above the wing. Increasing angle of attack 
strengthens the vortices until the induced wing pressure field and associated adverse streamwise 
pressure gradients cause vortex breakdown. The flow is further complicated as the leading edge 
vortices mix with the wake from the trailing edge downstream of the wing. The phenomenon of 
vortex breakdown (or vortex bursting) can have a significant influence on control surface 
performance and unsteady loading. The inherent unsteadiness of the breakdown process 
compounds the problem as it continually moves the breakdown region back and forth along the 
vortex axis. This creates serious time dependent flow problems and asymmetrically disposed 
breakdown positions above the wing that are aggravated with side-slip. 

Wide variations of breakdown patterns have been observed, and with increasing swirl the 
patterns change from spiral to near axisymmetric (Ref. 1). Spiral breakdown most commonly 
occurs over delta wings. In this breakdown process, the filament of fluid along the axis does not 
spread out symmetrically from a fixed stagnation point but, instead, takes on a spiral form around 
an unsteady "stagnation point" which varies in both space and time. Axisymmetric breakdown 
over delta wings, although rare, can also occur (Ref. 2). In this case, the vortex has a roughly 
axisymmetric breakdown pattern with a characteristic bubble which can have single or multiple tails 
(Ref. 3). 

Unfortunately, the parameters and conditions that result in vortex breakdown are poorly 

understood, because reliable quantitative experimental data are difficult to obtain. With limited 
experimental information to guide flow field modeling, numerical studies of vortex breakdown and 
control have met with only limited success (Ref. 4). There have been two principal reasons for 
this. In the first place, flow field unsteadiness associated with breakdown produces directional 
intermittency. This leads to large uncertainties in mean and unsteady flow measurements obtained 
with conventional pitot and hot wire probes. Secondly, and perhaps more important, is the fact 
that vortex breakdown is known to be extremely sensitive to any form of introduced disturbance. 
Probes, due to their blockage, may drastically alter the breakdown position. For these reasons, 
almost complete reliance has been placed on flow visualization techniques to determine flow field 
characteristics. In the past, air or hydrogen bubbles have been used as tracers to visualize flow 
patterns in water tunnels. For steady flows, streak lines can be identified with streamline patterns. 
However, in more complex flows of practical interest, the use of bubbles for flow visualization has 
distinct drawbacks. First of all, their introduction acts as a fluid lubricant which alters the apparent 
fluid viscosity and so its turbulent structure. Secondly, light refraction at the gas/water interfaces 
will destroy laser beam coherence and make it impossible to obtain laser velocimeter measurements 
in the regions where the tracer is present. But, with the advent of the laser fluorescence 
anemometer, there are now opportunities to determine accurate quantitative flow field velocity 
measurements of the vortex bursting process. 
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The Laser Fluorescence Anemometer 

The principle of operation of the laser fluorescence anemometer is shown schematically in 

Fig. 1. The mean velocity and turbulence measurements are made with a dual-beam velocimeter 
utilizing a Bragg cell that enables moving interference fringes to be generated in the focal volume 
so that instantaneous velocity magnitude and direction measurements can be achieved from the 
frequency shift (f^) around the incident and modulated laser beam interference frequency (f 0 ). 

i.e. U = X(f D - f 0 ) / 2Sin (0/2) where X is the wavelength of the incident laser light. Mean and 
fluctuating concentration measurements are achieved by observing the intensity of fluorescent light 
emitted from the focal volume at a different wavelength (X c ). At correct levels, tracer fluorescence 
is linearly proportional to the trace material concentration and, therefore, fluorescent intensity is 
directly proportional to the concentration of fluid from the seeded flow in the focal volume. The 
cathode current from a second photomultiplier tube is coupled to a high-gain current to voltage 
converter to produce a continuous voltage proportional to the instantaneous concentration. 

Since fluorescence is such a complex phenomenon dependent upon many parameters, only 
those of particular importance to the present application have been considered. In this context, the 
principal requirements were linearity combined with adequate sensitivity (i.e., signal/noise ratio), 
frequency response, and spatial resolution. Since the fluorescent intensity of a particular organic 
dye can be a strong function of the solvent, which in this case was room temperature water, other 
dye-solvent-temperature combinations would produce different ( and possibly increased) 
fluorescence. However, the fluorescent output was more than adequate for the present tests. In 
addition, the relationship between fluorescent intensity and concentration is of course exponential 
but, at the extremely low dye concentrations used in these experiments, a linear approximation 
could be made without introducing significant errors. 

The data in Fig. 2 show this linearity of the present technique, which employed dysodium 
fluorescein dye (a sodium salt of fluorescein which has been used for flow visualization for many 
years) as the trace material in water. Since, for a given dye concentration, the measured 
fluorescence is a function of laser beam intensity and collection optics, the ordinate of Fig. 2 has 
been plotted in arbitrary (voltage) units. The high signal-to-noise ratio that can be obtained for dye 
concentrations as low as 0.04 ppm (by weight) is illustrated in Fig. 3. The rise time (~ 50 |i.sec), 
which corresponds to the time taken to chop the laser beam, shows the adequate frequency 
response of the system. 

Since there is usually an overlap of the absorption and emission curves for most dye- 
solvent combinations it is possible that fluorescent photons emitted from the probe volume could 
be reabsorbed by the dye molecules that are between the probe volume and the collection optics. 
This effect could produce a range dependent signal. However, fluorescent intensity measurements 
previously obtained across an entire test section when filled with stagnant water at maximum dye 
concentration showed that these effects and those of possible beam attenuation were negligible. 
Thus, unlike absorption techniques which measure integrated ("line of sight") properties, the 
receiving optical arrangement primarily governs the spatial resolution of the fluorescence 
technique. In the present experiments, off-axis light collection and multimode fiber aperture size 
resulted in a maximum focal volume length dimension of approximately 0.5 mm although smaller 
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spatial resolution can be achieved by appropriate choice of collection optics without affecting the 
LDV, since the velocimeter interference fringes are moving. 

The instrument delivered to NASA LaRC comprises three primary elements namely: the 
optical, traverse control, and data acquisition systems. The Laser Optical system (Fig. 4) uses a 6 
watt Argon-ion laser. The transmitting optical arrangement (Fig. 5) is straightforward with a few 
unique features addressing the common problem of beam distortion or thermal blooming at higher 
laser powers. Frequency shifting is done before the color separation prisms using a single 
Acousto-Optic modulator made of a selected flint glass which can handle the full laser power with 
minimal distortion. (For significantly greater laser powers we also found a water cooled Bragg cell 
which was more than required for this application.) This is followed by the color separation 
prisms, the first of which are fused silica for power handling capacity. A final prism of dense flint 
provides maximum dispersion once the laser beam has been split into at least eight beams. Final 
color selection is made using right angle prisms. The lines used for this experiment were 
514.5nm, 488nm, 476.5nm. Although the emission spectrum of the fluorescence is centered 
about 5 15nm, there was sufficient higher wavelength emission for the edge filter to select 
wavelengths above 525nm. Other laser lines could have been selected if needed for fluorescence 
excitation or for separation from emission lines. 

Pure fused- silica core single-mode polarization-preserving fibers are used for light 
transmission; two fibers per color. The use of optical fibers not only avoids the tedium of mirror- 
traverse alignment, but also greatly simplifies the transmission of light to the third axis. The pure 
fused silica core fibers seem immune from the progressive transmission losses which are found in 
other fibers. Polarization preserving fibers provide greater modal stability when the fibers are 
flexed or manipulated. For mechanical protection the fibers are armored and contained within a 
conduit. Upon exiting the fibers the beams are collimated at 4.4 mm dia. with a separation of 
60mm. 

Forward-scattered light is collected with a single 80mm diameter lens and focused into a 
200 |im multi-mode optical fiber which conducts it to the color separation and signal detection box 
(Fig. 6). The fluorescence signal is split off with an edge filter. For maximum efficiency, a prism 
separation scheme is used rather than di-chroic filter and interference filters for the LDV signals. 

Experience has shown that accurate positioning is vital to a successful test program. 
Accurate positioning is complicated in an air- window- water environment due to the difference in 
refractive indices on either side of the water tunnel window. Refractive index problems are 
particularly acute for the third component beams which are transmitted at 45 degrees to the normal 
incidence of the four beam axis. Fig. 7 illustrates the problem, which may not be intuitively 
obvious. In order to traverse the measurement volume horizontally some distance in water, the 
two orthogonal optical components, axial(x) and model vertical(z), must be traversed some lesser 
amount in air. In order for the third or off-axis pair of beams to intersect the measurement volume 
at the required angle of thirty degrees, the beams must strike the window at an angle of 
approximately forty-five degrees. In addition, linear horizontal movement of the focal point of the 
third pair of beams requires lens movement on a sloping line. As the lens gets closer to the 
window, it will have to rise. The length of movement on this sloping path will be only about half 
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of the movement of the focal point in water. Thus, two lens systems must move on different paths 
at different speeds in order to maintain a coincident focal point. To sort this all out, two three-axis 
traverse tables were installed for computer controlled, algorithm driven positioning of the LDV 
probe volume and forward scatter collecting optics. A fourth traversing axis supports the off-axis 
transmitted beams. Position is maintained by a custom designed eight axis traverse controller with 
micro-stepping drives, optical encoder feedback, and limit switch safety stops. Details of the 
Traverse Control System are given in Appendix A. 

A Laser Velocimeter Data Acquisition System (LVDAS) has been designed. This 
instrument processes one to three channels of LDV data and digitizes up to four channels of analog 
data, one of which represents the concentration of dye. The instrument ensures coincidence and 
multiplexes the data to the computer. Velocities and analog channel values are displayed as well as 
data rates. Details are given in Appendix B. 

Fig. 8 shows the modified data handling system for the 3D laser fluorescence anemometer. 
The continuous though not necessarily non-zero output from the high-gain current to voltage 
converter is fed directly into an analog to digital converter to provide 12 bits of digital information 
at 50 kHz. In water flows, this was more than sufficient to provide essentially real time point 
concentration data in digital form. But data from the three component LV system were not 
continuous wave since particle arrival times in the focal volume were random. However, 
whenever valid and essentially coincident data were received on all LV channels, a necessary 
requirement for shear stress measurement, these velocities, along with the instantaneous 
concentration voltage, was recorded. From a series of such readings, mean and turbulent velocity 
and concentration profiles were determined along with the turbulent shear stress and 
velocity/concentration cross correlations. These latter cross products provide new information on 
turbulent mixing rates in complex flows. In addition, we are able to determine details of the 
concentration/turbulent intermittencies from ensemble averages generated for selected instantaneous 
concentration levels. This will shed quantitative light on the turbulent structure and entrainment of 
fluid originating from different points in the flow field. 

Experimental Details 

Test Facility 

The NASA Langley 16- by 24- Inch Water Tunnel is shown in Fig. 9. The tunnel has a 
vertical test section with an effective working length of about 4.5 ft. The velocity in the test section 
can be varied from 0 to 0.75 ft/sec., resulting in unit Reynolds numbers from 0 to 7.73 x 10 4 ft ~ l 
based on a water temperature of 78°F. The normal test velocity yielding smooth flow is 0.25 
ft/sec, resulting in unit Reynolds numbers of 2.58 x 10 4 ft _1 at 68°F. The model support system 
has deflection ranges of ±33° and ±15° in two planes of rotation. Rotation is accomplished Via 
electronic remote control, and visual indicators allow the user to set angles within about ±0.25°. 

The fluorescence seeding method for this investigation used fluorescein dye injected into 
the jet flow field from inside the model. Naturally-occuring particles in the water were used for 
seeding for the LDV part of the instrument. A representative size distribution provided by NASA 
is shown in Fig. 10. 
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Model 


The test model selected by NASA was a non-axisymmetric afterbody propulsion model and 
is shown in Fig. 11. This model is a scaled-down version of a model to be tested in the NASA 
Langley 16-foot Transonic Wind Tunnel as part of a computational fluid dynamics code validation 
study. It can be used to simulate nozzle exit velocity ratios typical of those in the wind tunnel 
study. It consists of a generic forebody with a non-axisymmetric boattail and nozzle. Water is 
injected into the interior of the model and exhausted through flow-conditioning foam ahead of the 
throat and exits through the nozzle. Fluorescein dye is introduced upstream of the model in the 
water supply tube for the jet. 

Sample Test Results 

Unfortunately, optical beam refraction problems caused by complex test section wall 
deformations («0. 10") under hydrostatic loading impeded laser beam alignment although limited 
data were obtained at zero angle of attack and a nominal jet exit to free stream velocity ratio of 1 .5. 
For these test conditions, the cross flow velocity components (v,w) were negligible as expected. 
However, the model did provide a flowfield in which the basic instrument concepts could be 
verified. Fig. 12 shows the measured axial velocity distributions. It shows the extent of the jet 
and the location of the velocity defects in the model wall wakes. Clearly, the jet is highly skewed, 
but this has since been corrected by subsequent modification of the internal model flow treatment 
devices. The mean concentration profile (F 

g. 13), which is much more symmetric, defines the extent of the jet at this axial station. Some 
insight into the unsteady features of the flow can be obtained from the fluctuation measurements. 
Figs. 14 and 15 show there is significant mixing in the outer shear layers between the jet and 
freestream flows. Peak fluctuation levels in the regions of maximum mean gradients indicate that 
small scale mixing is dominant. A measure of the level of streamwise mixing in the jet can be 
determined from Fig. 16. As expected the u’c’ cross correlation function is positive as faster 
moving jet fluid is associated with higher concentration whereas fluid originating from the slower 
moving freestream has lower or zero concentration. 

Concluding Remarks 

A Laser Fluorescence Anemometer which comprises a three component laser doppler 
velocimeter system with a fourth channel to measure fluorescent dye concentration has been 
installed in the NASA Langley Water Tunnel. The system includes custom designed optics, data 
acquisition, and traverse control instruments and a custom software package. 

Feasibility studies clearly demonstrate how water tunnels can be used in conjunction with 
advanced optical techniques to provide non-intrusive detailed flow field measurements of complex 
fluid flows with a minimum of expense. The measurements show that the laser fluorescence 
anemometer will provide new insight into the structure, entrainment and mixing of vortical and 
shear layer flows. 
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Fig. 4a Schematic of Plan View LFA System 
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Fig. 6 Receiving Optics 
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Fig. 8 Data Acquisition and Traverse Control Systems 















Fig. 9 Langley 16- by 24- Inch Water Tunnel 
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Appendix A. The Traverse Control System 


The traverse control system is made up of four subsystems, see Fig. Al. The first 
subsystem is the main data taking computer, the HP 9000-330. The second subsystem, the TCS8 
(Traverse Control System 8 Axis), receives high level traverse commands from the HP 9000. The 
full duplex serial communications that links these two subsystems allows the HP 9000 to monitor 
the position and status of each axis in the system, see Appendix A.3 TCS8’s Serial Interface 
Command Descriptions. The TCS8 can also function as a stand alone traverse controller. 

Through the use of the TCS8’s front panel, an operator can execute all of the commands that the 
HP 9000 can, plus the operator can control all axes in jog mode, see Appendix A.1 TCS8’s Front 
Panel Descriptions and Appendix A.2 TCS8’s Local Command Descriptions. The third 
subsystem, the MDS (Motor Drive System), is controlled solely by the TCS8. The TCS8 
translates the high level commands from the HP 9000 and its front panel into low level indexer 
commands, see The Compumotor AX Drive User Manual previously delivered. The TCS8 also 
receives encoder pulses from the traverses via the MDS. This allows the TCS8 to display realtime 
position information on its front panel. The fourth and final subsystem of the traverse control 
system is the slide, motor, encoder, and limit switches that make up each axis. A drawing of each 
cable that is used to connect the traverse control system is included in Appendix A.4 Traverse 
Control System Cables. 

The TCS8 

The TCS8 is a microprocessor controlled system designed to interface an operator with a 
traverse system. The operator can utilize the TCS8 through the front panel, see Appendix A. 1 
TCS8’s Front Panel Descriptions and Appendix A.2 TCS8’s Local Command Descriptions, and/or 
with one or two host computers over serial interfaces, see Appendix A.3 TCS8’s Serial Interface 
Command Descriptions. The TCS8 stores all the critical parameters of motion, for each of the 
eight axes that it controls, in non-volatile memory. The critical parameters of motion being: 
position, encoder counts per unit travel, encoder counts per motor revolution, velocity, and 
acceleration. All of these parameters may be viewed, set, and saved. The TCS8 has three modes 
of motion. They are absolute, relative, and jogged. With absolute movements, the operator 
specifies the final location; with relative movements, a distance is specified; and with jogged 
movements the operator presses a jog key on the front panel of the TCS8 until the desired location 
is obtained. 

The M[otor Drive System 

Each of the two MDS ’ s have 4 indexer/drivers contained within them. The TCS8 
communicates with the indexers in the MDS’s over a closed loop serial daisy chain. When two 
MDS’s are used, as in this setup, the first MDS in the chain must be set to 8 and the second set to 
4. By setting the first MDS to 8, the operator is opening the closed loop serial daisy chain 
allowing the second MDS to be included in the chain. The 4/8 switch is located on the back panel 
of each MDS, see Fig. A2. This figure also shows the location of all the motor, limit, and encoder 
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connections. Channels XI, X2, Yl, and Y2 of the TCS8 control axis 1 through 4 on the first 
MDS and channels Zl, Z2, Al, and A2 control axis 1 through 4 on the second MDS. The TCS8 
Encoders connector on the back of each MDS has a corresponding connector of the back of the 
TCS8, see Fig. A3 Schematic of TCS8 Back Panel. The interconnecting cable is detailed in 
Appendix A. 4 Traverse Control System Cables. 

Positioning Resolution 

The indexer/drivers that are used in the MDS can drive the motors at 12,800 
steps/revolution. The encoder used on each axis are 1000 pulses/revolution with quadrature 
encoding. Quadrature encoding adds a factor of 4 to the number of pulses/revolution to make this 
number 4000 pulses/revolution. This number, 4000 pulses/revolution, is well within the limit of 
12,800 steps/revolution set by the indexer. The final factor in the product of the resolution of an 
axis is the number of threads/inch of the lead screw. All of the axes of the traverse system, except 
the auxiliary axis, have lead screws of 5 threads/inch, the auxiliary axis has a lead screw of 10 
threads/inch. So the positioning resolution of the axes with a 5 threads/inch lead screw is 0.00005 
inches and the auxiliary axis has a resolution of 0.000025 inches. 
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Fig. A1 Langley Traverse Control System 





MOTOR 1 


LIMIT 1 


ENCODER 1 


TCS8 

ENCODER 



Fig. A2 Schematic of Motor Drive System Back Panel 
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1. POSITION DISPLAY WINDOWS 

There are eight windows corresponding to the eight axes that the TCS8 is capable of controlling. 
The position of each axis is continuously updated, by monitoring its encoder, and displayed in a 
fixed format of a sign, two digits, a decimal point, and four digits. 

2. POWER KEY 

The power key is used to store the current configuration to nonvolatile memory before turning off 
power to the TCS8. Pressing the power key turns the displays off and saves the current 
configuration. Pressing it again turns the displays back on. This key can be used to implement a 
screen saver function. 

3. JOG CONTROL KEYS 

These keys are used to control up to eight axes in a jog mode. The mode (slaved, one’s only, or 
two’s only) can be set through the jog menu. When the operator presses a jog key, the respective 
axis will begin to move. The direction that the axis moves is determined by the operator pressing 
either a plus or minus jog key. A plus jog key will turn the lead screw in a clockwise direction 
(away from the motor), a minus jog key will turn it in the counter-clockwise direction (towards the 
motor). By releasing the jog key the operator stops motion on that axis. Motion will also stop, if 
the axis reaches the limit for the direction it is moving, or if the indexer determines that the axis has 
stalled. 
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4. SCROLL KEYS 

These keys are used to scroll items through the MENU, COMMAND, and CHANNEL windows. 
All of the menus, their commands, and channel variations will be detailed in another appendix. 

5. COMMAND WINDOWS 

These three windows (MENU, COMMAND, and CHANNEL) are used, in tandem with their 
respective scroll keys, to formulate a command to be executed by the TCS8. 

6. EXECUTE KEY 

This key is used to execute the command currently formulated in the MENU, COMMAND, and 
CHANNEL windows. 

7. DATA WINDOW 

Many of the TCS8’s commands require some added data, e.g. the distance to move or an axis’ 
encoder counts per unit. Data for these commands is entered from the numeric key pad on the 
lower right of the TCS8 into the DATA window. Only a valid real number can be entered into the 
DATA window. If the operator enters an invalid real number the character that is invalid will flash 
until the operator presses backspace or a valid character. 

8. STOP KEY 

The stop key, when pressed, will stop motion on all axes. The TCS8 will not loose track of the 
position of any axis. A move command started by the host computer and stopped by the stop key 
will finish normally with the position being reported. The position reported is the instantaneous 
position when the stop key was pressed. The final position of the axis being moved could be 
different than what was reported thus the host computer should read the position again after a panic 
stop. 

9. STATUS WINDOW 

The STATUS window reflects the result of all commands. For commands that are not 
instantaneous, this window displays a busy status and then when the command completes it 
displays a ready status. The results of all view commands are displayed in the STATUS window. 
The STATUS window also displays the activity over the COM interfaces. For example, when the 
command for viewing position is sent over the COM1 interface, the STATUS window will display 
“COM1 VP” and when the command is completed the window will display “COM1 vp”. 

10. NUMERIC KEY PAD 

The numeric key pad is used to enter a number into the data window. The user may backspace in 
the window or clear (shift- backspace) the window. Only a valid real number can be entered into 
the data window. If the operator enters an invalid real number the character that is invalid will 
flash until the operator presses backspace or a valid character. 
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Appendix A.2 TCS8’s Local Command Descriptions 

This appendix describes the command set that can be executed from the front panel of the 
TCS8. Using the up and down keys under the MENU, COMMAND, and CHANNEL windows, 
the operator can formulate a command and then execute it by pressing the EXECUTE key. Some 
commands require extra information to be entered into the DATA window through the use of the 
numeric key pad. Each description includes a list of related commands that should be refer to to 
enhance the operator’s understanding of the command. Also where applicable, the default setting 
is given. 

MENU: MOVE 
COMMAND: TO ZERO 

CHANNELS: ALL, ONE’S, TWO’S, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, 
Zl, Z2, Al, A2 

DESCRIPTION : The MOVE TO ZERO command is an easy way to move some or all of the 
axes to the zero position. This command can also be accomplished with the MOVE ABSOLUTE 
command and a zero in the DATA window. Before using this command, selected axes must be 
initialized with the INIT DRIVE ON command. This command can be canceled by pressing the 
STOP key. When the STOP key is pressed, all axes will stop immediately. If an axis encounters a 
limit before reaching zero, the rest of its movement is aborted. 

RELATED COMMANDS: MOVE ABSOLUTE, MOVE RELATIVE, INIT Drive ON 


MENU: MOVE 


COMMAND: ABSOLUTE 


CHANNELS: X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 

DESCRIPTION: The MOVE ABSOLUTE command requires a position to be entered in the 
DATA window. This position and the current position of the axis is used to calculate the relative 
distance the axis must move. Before using this command, selected axes must be initialized with 
the INIT DRIVE ON command. This command can be canceled by pressing the STOP key. 

When the STOP key is pressed, all axes will stop immediately. If an axis encounters a limit before 
reaching the position entered in the DATA window, the rest of its movement is aborted. 

RELATED COMMANDS : MOVE TO ZERO, MOVE RELATIVE, INIT Drive ON 
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MENU: MOVE 


COMMAND: RELATIVE 


CHANNELS: X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 


DESCRIPTION : The MOVE RELATIVE command requires a distance to be entered in the 
DATA window. This position is used to calculate the relative distance the axis must move. Before 
using this command, selected axes must be initialized with the INIT DRIVE ON command. This 
command can be canceled by pressing the STOP key. When the STOP key is pressed, all axes will 
stop immediately. If an axis encounters a limit before moving the distance entered in the DATA 
window, the rest of its movement is aborted. 

RELATED COMMANDS : MOVE TO ZERO, MOVE ABSOLUTE, INIT Drive ON 


MENU: JOG 


COMMAND: MODE 


CHANNELS: SLAVED, ONE’S, TWO’S 

DESCRIPTION: The JOG MODE command sets the way the JOG keys operate. When 
SLAVED is the setting, both the one and two axis of the X, Y, Z, or A coordinate will move the 
same amount. When ONE’S is the setting, only the one axes of the X, Y, Z, or A coordinate will 
move. And finally, when TWO’S is the setting, only the two axes of the X, Y, Z, or A coordinate 
will move. The current mode is marked with an asterisk. 

RELATED COMMANDS: none 


DEFAULT: SLAVED 
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MENU: SET 


COMMAND: CPU 

CHANNELS: ALL, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 


DESCRIPTION : The SET CPU command allows the user to change the counts per unit travel. 
The CPU for an axis is determined by multiplying the encoder resolution (counts/revolution) by the 
lead screws resolution (revolutions/unit of travel). A units conversion can be added here to change 
for example from inches to centimeters. When the CPU for an axis is changed, the position is 
automatically converted. This command requires a value to be entered in the DATA window. 

RELATED COMMANDS: SET CPR, SET POSITION 


DEFAULT: XI 

20000 

X2 

20000 

Yl 

20000 

Y2 

20000 

Zl 

20000 

Z2 

20000 

Al 

40000 

A2 

40000 


MENU: SET 
COMMAND: CPR 

CHANNELS: ALL, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 


DESCRIPTION : The SET CPR command allows the user to change the encoder counts per 
motor revolution. The CPR for an axis is determined by dividing the encoder resolution 
(counts/revolution) by the lead screws resolution (revolutions/unit of travel). The encoder counts 
per motor revolution, that is entered in the DATA window, must be a positive integer. 

RELATED COMMANDS: SET CPU 


DEFAULT: XI 

4000 

X2 

4000 

Yl 

4000 

Y2 

4000 

Zl 

4000 

Z2 

4000 

Al 

4000 

A2 

4000 
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MENU: SET 


COMMAND: POSITION 

CHANNELS: ALL, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 


DESCRIPTION : The SET POSITION command allows the user to change the current position 
of an axis. The new position must be entered in the DATA window be for executing the command. 

RELATED COMMANDS: SET CPU 


MENU: SET 
COMMAND: VELOCITY 


CHANNELS: ALL, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 

DESCRIPTION : The SET VELOCITY command allows the user to change the maximum speed 
at which an axis will travel. The range of valid velocities is 0.002 to 50.000 revolutions per 
second. The default is 5 revs/sec. An axis may stall at velocities higher than the default. The new 
velocity must be entered in the DATA window be for executing the command. 

RELATED COMMANDS: SET ACCEL. 


DEFAULT: XI 

5.000 

X2 

5.000 

Yl 

5.000 

Y2 

5.000 

Zl 

5.000 

Z2 

5.000 

Al 

5.000 

A2 

5.000 
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MENU: SET 


COMMAND: ACCEL. 

CHANNELS: ALL, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, Zl, Z2, Al, A2 

\ 

DESCRIPTION : The SET ACCEL, command allows the user to change the maximum 
acceleration for an axis. The range of valid accelerations is 0.01 to 999.99 revolutions per second 
per second. The default is 5 revs/sec/sec. The new acceleration must be entered in the DATA 
window be for executing the command. 

RELATED COMMANDS: SET VELOCITY 


DEFAULT: XI 

5.00 

X2 

5.00 

Yl 

5.00 

Y2 

5.00 

Zl 

5.00 

Z2 

5.00 

Al 

5.00 

A2 

5.00 


MENU: SET 
COMMAND: CmtsOn 


CHANNELS: ALL, ONE'S, TWO'S, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, 
Z1,Z2, Al, A2 

DESCRIPTION : The SET CmtsOn command allows the user to turn the motor currents on. 
The motor current must be on for an axis to be moved. The information in the DATA window is 
ignored. 

RELATED COMMANDS: SETCmtsOff 
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MENU: SET 


COMMAND: CmtsOff 

CHANNELS: ALL, ONE’S, TWO’S, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, 
Zl, Z2, Al, A2 


DESCRIPTION : The SET CmtsOff command allows the user to power down motors when 
they will not be used for long periods of time. The information in the DATA window is ignored. 

RELATED COMMANDS : SET CmtsOn 


MENU: SET 
COMMAND: INITS ON 

CHANNELS: ALL, ONE’S, TWO’S, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, 
Zl, Z2, Al, A2 

DESCRIPTION : The SET INITS ON command allows the user to initialize the indexers 
without turning on the power to the motors. The information in the DATA window is ignored. 

RELATED COMMANDS: INIT Drive ON 
MENU: VIEW 
COMMAND: Cnt/Unit 

.CHANNELS: XI, X2, Yl, Y2, Zl, Z2, Al, A2 

DESCRIPTION : The VIEW Cnt/Unit command displays the current setting of the encoder 
counts per unit travel parameter for the selected axis in the STATUS window. The information in 
the DATA window is ignored. 

RELATED COMMANDS: SET CPU 
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MENU: VIEW 


COMMAND: Cnt/MRev 

CHANNELS: XI, X2, Yl, Y2, Zl, Z2, Al, A2 

DESCRIPTION : The VIEW Cnt/MRev command displays the current setting of the encoder 
counts per motor revolution parameter for the selected axis in the STATUS window. The 
information in the DATA window is ignored. 

RELATED COMMANDS: SETCPR 


MENU: VIEW 

COMMAND: VELOCITY 

CHANNELS: XI, X2, Yl, Y2, Zl, Z2, Al, A2 

DESCRIPTION : The VIEW VELOCITY command displays the current setting of the velocity 
parameter for the selected axis in the STATUS window. The information in the DATA window is 
ignored. 

RELATED COMMANDS: SET VELOCITY 


MENU: VIEW 
COMMAND: ACCEL. 

CHANNELS: XI, X2, Yl, Y2, Zl, Z2, Al, A2 

DESCRIPTION : The VIEW ACCEL, command displays the current setting of the acceleration 
parameter for the selected axis in the STATUS window. The information in the DATA window is 
ignored. 

RELATED COMMANDS: SET ACCEL. 
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MENU: VIEW 


COMMAND: INIT 


CHANNELS: none 


DESCRIPTION : The VIEW INIT command uses the STATUS window to display a 
one(initialized) or a zero(uninitialized) for each axis. The STATUS window has eight characters; 
left to right respectively reflecting the status of: XI, X2 Al, A2. The information in the DATA 
window is ignored. 

RELATED COMMANDS: SET INITS, INIT Drive ON 


MENU: VIEW 


COMMAND: CURRENTS 


CHANNELS: none 


DESCRIPTION : The VIEW CURRENTS command uses the STATUS window to display a 
one(current on) or a zero(current off) for each axis. The STATUS window has eight characters; 
left to right respectively reflecting the status of: XI, X2 ..., Al, A2. The information in the DATA 
window is ignored. 

RELATED COMMANDS: SET CmtsOn, SET CmtsOff, INIT Drive ON, INIT Drive OFF 


MENU: VIEW 


COMMAND: PlusLMT 


CHANNELS: none 


DESCRIPTION : The VIEW Plus LMT command uses the STATUS window to display a 
one(on limit) or a zero(not on limit) for each axis. The STATUS window has eight characters; left 
to right respectively reflecting the status of: XI, X2 ..., Al, A2. The information in the DATA 
window is ignored. 

RELATED COMMANDS: none 
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MENU: VIEW 


COMMAND: Minus LMT 


CHANNELS: none 


DESCRIPTION : The VIEW Minus LMT command uses the STATUS window to display a 
one(on limit) or a zero(not on limit) for each axis. The STATUS window has eight characters; left 
to right respectively reflecting the status of: XI, X2 Al, A2. The information in the DATA 
window is ignored. 

RELATED COMMANDS: none 


MENU: VIEW 
COMMAND: HOME 
CHANNELS: none 


DESCRIPTION : The VIEW HOME command uses the STATUS window to display a one(on 
limit) or a zero(not on limit) for each axis. The STATUS window has eight characters; left to right 
respectively reflecting the status of: XI, X2 ..., Al, A2. The information in the DATA window is 
ignored. 

RELATED COMMANDS: none 


MENU: VIEW 
COMMAND: STALL 
CHANNELS: none 


DESCRIPTION : The VIEW STALL command uses the STATUS window to display a 
one(stalled) or a zero(not stalled) for each axis. The STATUS window has eight characters; left to 
right respectively reflecting the status of: XI, X2 ..., Al, A2. The information in the DATA 
window is ignored. 

RELATED COMMANDS: none 
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MENU: INIT 


COMMAND: DEFAULT 
CHANNELS: none 

DESCRIPTION : The INIT DEFAULT command restores the initial factory defaults (CPU, 
CPR, VELOCITY, ACCELERATION, BAUD RATE, BITS/CHAR, PARITY,STOP 
BITS .HANDSHAKE) of the TCS8. After executing this command, execute the command INIT 
Drive ON to initialize the indexers. The information in the DATA window is ignored. 

RELATED COMMANDS: SET CPU, SET CPR, SET VELOCITY, SET ACCEL. 


MENU: INIT 
COMMAND: Drive ON 

CHANNELS: ALL, ONE’S, TWO'S, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, 
Zl, Z2, Al, A2 

DESCRIPTION : The INIT Drive ON command initializes the selected axes for movement. 
After executing this command the currents are on to the motors. The information in the DATA 
window is ignored. 

RELATED COMMANDS: SET CPU, SET CPR, SET VELOCITY, SET ACCEL., SET 
CnrtsOn, SET CnrtsOff, INIT DEFAULT 


MENU: INIT 
COMMAND: Drive OFF 

CHANNELS: ALL, ONE'S, TWO’S, X1&X2, Y1&Y2, Z1&Z2, A1&A2, XI, X2, Yl, Y2, 
Zl, Z2, Al, A2 

DESCRIPTION : The INIT Drive OFF command is an alias for SET CmtsOff. 

RELATED COMMANDS: SET CmtsOff 
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MENU: C0M1/C0M2 


COMMAND: BaudRate 

CHANNELS: 19.2K, 9600, 4800, 2400, 1200, 300, 110 

DESCRIPTION : The COM1/COM2 BaudRate command set the baud rate for the selected 
communication channel. The information in the DATA window is ignored. The current baud rate 
is marked with an asterisk. 

RELATED COMMANDS: none 
DEFAULT: 9600 


MENU: COM1/COM2 
COMMAND: Bit/Char 
CHANNELS: SEVEN, EIGHT 

DESCRIPTION : The COM1/COM2 Bit/Char command set the bits per character for the selected 
communication channel. The information in the DATA window is ignored. The current number of 
bits per character is marked with an asterisk. 

RELATED COMMANDS: none 
DEFAULT: EIGHT 


MENU: COM1/COM2 
COMMAND: Parity 
CHANNELS: NONE, EVEN, ODD 


DESCRIPTION : The COM1/COM2 Parity command set the parity for the selected 
communication channel. The information in the DATA window is ignored. The current parity is 
marked with an asterisk. 

RELATED COMMANDS: none 


DEFAULT: EVEN 
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MENU: C0M1/C0M2 


COMMAND: StopBits 
CHANNELS: 1,15,2 

DESCRIPTION : The COM 1/COM2 StopBits command set the stop bits for the selected 
communication channel. The information in the DATA window is ignored. The current number of 
stop bits is marked with an asterisk. 

RELATED COMMANDS: none 
DEFAULT: 1 


MENU: COM1/COM2 
COMMAND: HandShak 
CHANNELS: NO, YES 

DESCRIPTION : The COM1/COM2 HandShak command set the handshake for the selected 
communication channel. The information in the DATA window is ignored. An asterisk marks 
whether there is handshaking or not. 

RELATED COMMANDS: none 


DEFAULT: YES 



Appendix A.3 TCS8’s Serial Interface Command Descriptions 

This appendix describes the command set that can be executed through the serial interfaces 
of the TCS8. Each description includes a code section that outlines the characters that must be sent 
to execute the command. The vertical bar in this section is used as a separator and is not sent as 
part of the command code. The symbol “CRLF” stands for the two characters carriage return and 
line feed. Also where applicable, the default setting is given. 


COMMAND: CHANGE SERIAL CONFIGURATION 


CODE: CS COM;CATEGOR Y ; ATTRIBUTE; 


PARAMETERS: 


COM: 

1/COM 1 


2/COM2 

CATEGORY: 

0/BAUDRATE 

ATTRIBUTE: 

0/19.2K 


1/9600 


2/4800 


3/2400 


4/1200 


5/300 


6/110 

CATEGORY: 

1(BITS PER CHARACTER) 

ATTRIBUTE: 

0/SEVEN 


1/EIGHT 

CATEGORY: 

2(PARITY) 

ATTRIBUTE: 

0/NONE 


1/EVEN 


2/ODD 

CATEGORY: 

3(STOP BITS) 

ATTRIBUTE: 

0/ONE 


1/ONE AND A HALF 


2/TWO 

CATEGORY: 

4(HANDSHAKE) 

ATTRIBUTE: 

0/NO 


1/YES 


DESCRIPTION : This command must be executed with extreme caution and thought. If the 
user changes an attribute of the same COM port that he is sending the command, he must change to 
that attribute on the host computer before sending the next command. The best way to change the 
serial configuration of a COM port is to utilize the front panel commands. 

DEFAULT: 9600 baud, EIGHT bits/char, EVEN parity, ONE stop bit, handshaking YES 

EXAMPLE: To change the baudrate of COM1 to 2400 the user must send CS1;0;3; 
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COMMAND: MOVE TO ABSOLUTE POSITION AND REPORT FINAL POSITION 

CODE: MA CHANNEL:POSITION,CHANNEL:POSITION,...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

POSITION: Real number free format 

DESCRIPTION : This command moves selected channels to absolute positions. 

EXAMPLES: 

To move all channels to zero the user may send MA0:0,CRLF or MA12345678:0,CRLF 
To move channel X 1 to zero the user must send MA 1 :0,CRLF 

To move channels XI and X2 to zero the user may send MA12:0,CRLF or MA1:0,2:0,CRLF or 
MAI :0,CRLF and MA2:0,CRLF 


COMMAND: MOVE TO RELATIVE DISTANCE AND REPORT FINAL POSITION 

CODE: MR CHANNEL:DISTANCE,CHANNEL:DISTANCE,...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4m 

5/Z1 

6/Z2 

7/A1 

8/A2 

POSITION: Real number free format 


DESCRIPTION : This command moves selected channels relative distances. 

EXAMPLES: 

To move all channels one unit the user may send MR0:1,CRLF or MR12345678:1,CRLF 
To move channel X 1 one unit the user must send MR 1 : 1 ,CRLF 

To move channels XI and X2 one unit the user may send MR 12: 1,CRLF or MR1:1,2:1,CRLF or 
MR 1 : 1 ,CRLF and MR2: 1 ,CRLF 
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COMMAND: SET ACCELERATION 

CODE: SA CHANNEL:ACCELERATION, CHANNEL-ACCELERATION, ...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

ACCELERATION: Real number free format between 
0.01 and 99.99 inclusive. 

DESCRIPTION: This command sets the acceleration for selected channels. 

DEF AULT : All channels 5.00 revolutions/second/second 

EXAMPLES: 

To set the acceleration for all channels to 4.00 revolutions/second/second the user may send 

SA0:4.00,CRLF or SA12345678:4.00,CRLF 

To set the acceleration for channel XI to 4.00 revolutions/second/second the user must send 

SA1:4.00,CRLF 

To set the acceleration for channels XI and X2 to 4.00 revolutions/second/second the user may 

send SA12:4.00,CRLF or SA1:4.00 ,2:4.00, CRLF or SA1:4.00,CRLF and SA2:4.00,CRLF 


COMMAND: VIEW ACCELERATION 

CODE: VA CHANNELICHANNEL...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/yi 

4/Y2 

5/Z1 

6172 

7/A1 

8/A2 

DESCRIPTION : This command views the acceleration for selected channels. The TCS8 
transmits each of the accelerations requested back to the host computer separated by carriage return 
line feeds. 

EXAMPLES: 

To view the acceleration for all channels the user may send VA0CRLF or VA12345678CRLF 
To view the acceleration for channel XI the user must send VA ICRLF 

To view the acceleration for channels XI and X2 the user may send VA12CRLF or VA ICRLF and 
VA2CRLF 
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COMMAND: SET VELOCITY 

CODE: SV CHANNEL:VELOCITY,CHANNEL:VELOCITY,...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

VELOCITY : Real number free format between 

0.001 and 50.000 inclusive, 

DESCRIPTION : This command sets the velocity for selected channels. 

DEF AULT : All channels 5.000 revolutions/second 

EXAMPLES: 

To set the velocity for all channels to 4.00 revolutions/second the user may send SV0:4.00,CRLF 

or S V12345678:4.00,CRLF 

To set the velocity for channel XI to 4.00 revolutions/second the user must send SV1:4.00,CRLF 

To set the velocity for channels XI and X2 to 4.00 revolutions/second the user may send 

SV12:4.00,CRLF or SV1:4.00 ,2:4.00,CRLF or SV1:4.00,CRLF and SV2:4.00,CRLF 


COMMAND: VIEW VELOCITY 

CODE: VV CHANNELICHANNEL...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

DESCRIPTION : This command views the velocity for selected channels. The TCS8 transmits 
each of the velocities requested back to the host computer separated by carriage return line feeds. 

EXAMPLES: 

To view the velocity for all channels the user may send VV0CRLF or VV 12345678CRLF 
To view the velocity for channel XI the user must send VV1CRLF 

To view the velocity for channels XI and X2 the user may send VV12CRLF or VV1CRLF and 
VV2CRLF 
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COMMAND: SET ENCODER COUNTS PER UNIT TRAVEL 

CODE: SU CHANNEL: CPU,CHANNEL:CPU,.. . ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

CPU: Non-zero real number free format. 

DESCRIPTION : This command sets the encoder counts per unit travel for selected channels. 
DEFAULT: X1,2,Y1,Y2,Z1,Z2 20000 counts/inch and A1,A2 40000 counts/inch 
EXAMPLES: 

To set the encoder counts per unit travel for all channels to 5000 the user may send 
SU0:5000,CRLF or SU12345678:5000,CRLF 

To set the encoder counts per unit travel for channel XI to 5000 the user must send 
SU 1:5000, CRLF 

To set the encoder counts per unit travel for channels XI and X2 to 5000 the user may send 
SU12:5000,CRLF or SU1:5000 ,2:5000, CRLF or SU1:5000,CRLF and SU2:5000,CRLF 


COMMAND: VIEW ENCODER COUNTS PER UNIT TRAVEL 

CODE: VU CHANNELICHANNEL...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

DESCRIPTION : This command views the encoder counts per unit travel for selected channels. 
The TCS8 transmits each of the encoder counts per unit travel requested back to the host computer 
separated by carriage return line feeds. 

EXAMPLES: 

To view the encoder counts per unit travel for all channels the user may send VUOCRLF or 
VU12345678CRLF 

To view the encoder counts per unit travel for channel XI the user must send VU1CRLF 
To view the encoder counts per unit travel for channels XI and X2 the user may send VU 12CRLF 
or VU1CRLF and VU2CRLF 
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COMMAND: SET ENCODER COUNTS PER MOTOR REVOLUTION 

CODE: SR CHANNEL: CPR, CHANNEL: CPR,. . . ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

Am 

5/Zl 

6/Z2 

7/A1 

8/A2 

CPU : Non-zero integer free format. 

DESCRIPTION : This command sets the encoder counts per motor revolution for selected 
channels. 

DEFAULT: X1,2,Y1,Y2,Z1,Z2 and A1,A2 4000 counts/inch 
EXAMPLES: 

To set the encoder counts per motor revolution for all channels to 500 the user may send 
SR0:500,CRLF or SR12345678:500,CRLF 

To set the encoder counts per motor revolution for channel X 1 to 500 the user must send 
SR1:500,CRLF 

To set the encoder counts per motor revolution for channels XI and X2 to 500 the user may send 
SR12:500,CRLF or SRI :500 ,2:500, CRLF or SRI :500,CRLF and SR2:500,CRLF 


COMMAND: VIEW ENCODER COUNTS PER MOTOR REVOLUTION 

CODE: VR CHANNELICHANNEL...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Zl 

6/72 

7/A1 

8/A2 

DESCRIPTION : This command views the encoder counts per motor revolution for selected 
channels. The TCS8 transmits each of the encoder counts per motor revolution requested back to 
the host computer separated by carriage return line feeds. 

EXAMPLES: 

To view the encoder counts per motor revolution for all channels the user may send VR0CRLF or 
VR 1 234567 8CRLF 

To view the encoder counts per motor revolution for channel XI the user must send VR1CRLF 
To view the encoder counts per motor revolution for channels XI and X2 the user may send 
VR12CRLF or VR ICRLF and VR2CRLF 
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COMMAND: SET POSITION 

CODE: SP CHANNEL:POSITION,CHANNEL:POSITION,...ICRLF 


PARAMETERS: CHANNEL: 


POSITION: 


0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

real number. 


DESCRIPTION : This command sets the position for selected channels. 


EXAMPLES: 

To set the position for all channels to 1 .5 the user may send SPO: 1 .5,CRLF or 
SP12345678: 1.5, CRLF 

To set the position for channel XI to 1.5 the user must send SP1:1.5,CRLF 

To set the position for channels XI and X2 to 1.5 the user may send SP12:1.5,CRLF or SP1:1.5 

,2:1.5, CRLF or SP1:1.5,CRLF and SP2:1.5,CRLF 


COMMAND: VIEW POSITION 

CODE: VP CHANNELICHANNEL...ICRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

DESCRIPTION : This command views the position for selected channels. The TCS8 transmits 
each of the positions requested back to the host computer separated by carriage return line feeds. 

EXAMPLES: 

To view the position for all channels the user may send VPOCRLF or VP12345678CRLF 
To view the position for channel XI the user must send VP 1 CRLF 

To view the position for channels XI and X2 the user may send VP12CRLF or VP1CRLF and 
VP2CRLF 
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COMMAND: SET CURRENT TO MOTOR WINDINGS 

CODE: SC CHANNEL:ON/OFF,CHANNEL:ON/OFE,.JCRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

ON/OFF: 1/ON 

O/OFF 

DESCRIPTION : This command sets the current to the motor windings for selected channels on 
or off. 

EXAMPLES: 

To set the current to the motor windings for all channels on the user may send SC0:1,CRLF or 
SC12345678:1,CRLF to set them off the user may send SC0:0,CRLF or SC12345678:0,CRLF 


COMMAND: VIEW CURRENT TO MOTOR WINDINGS 

CODE: VC CHANNELICHANNEL. JCRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/72 

7/A1 

8/A2 

DESCRIPTION : This command views the current to the motor windings for selected channels. 
The TCS8 transmits each response of on/off (1/0) back to the host computer separated by carriage 
return line feeds. 

EXAMPLES: 

To view the current to the motor windings for all channels the user may send VCOCRLF or 
VC12345678CRLF 

To view the current to the motor windings for channel XI the user must send VC1CRLF 
To view the current to the motor windings for channels XI and X2 the user may send VC12CRLF 
or VC1CRLF and VC2CRLF 
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COMMAND: SET INITIALIZATION OF INDEXER/DRIVERS 

CODE: SI CHANNELICHANNEL.JCRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

DESCRIPTION : This command sends the current value of the acceleration, velocity, and the 
encoder counts per motor revolution to the indexer/driver for the selected channels. This command 
must be sent before any move commands may be sent. 

EXAMPLES: 

To initialize all channels the user may send SIOCRLF or SI 1 234567 8CRLF 
To initialize channel XI the user must send SI1CRLF 

To initialize channels XI and X2 the user may send SI12CRLF or SI1CRLF and SI2CRLF 


COMMAND: VIEW INITIALIZATION OF INDEXER/DRIVERS 

CODE: VI CHANNELICHANNEL.JCRLF 

PARAMETERS: CHANNEL: 0/ALL CHANNELS 

1/X1 

2/X2 

3/Y1 

4/Y2 

5/Z1 

6/Z2 

7/A1 

8/A2 

DESCRIPTION : This command returns “1” if the indexer/driver has been initialized since the 
TCS8 was turned on and “0” if it has not. The TCS8 transmits each of the responses back to the 
host computer separated by carriage return line feeds. 

EXAMPLES: 

To check the initialization of all channels the user may send VIOCRLF or VI12345678CRLF 
To check the initialization of channel XI the user must send VI1CRLF 

To check the initialization of channels XI and X2 the user may send VI12CRLF or VI1CRLF and 
VI2CRLF 
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Appendix A.4 

Traverse Control System Cables 
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MOTOR 

DRIVE 

SYSTEM 


& 



DWG# LFA-TAA-002 - 
DWG# LFA-TAA-003 - 
^ — DWG# LFA-TAA-004 - 
DWG# LFA-TAA-002 - 
S-DWG# LFA-TAA-003 - 
DWG# LFA-TAA-004 - 
DWG# LFA-TAA-002 
r — DWG# LFA-TAA-003 
N— DWG# LFA-TAA-004 


\ 


SERIAL 



~ INTERFACES 


S-DWG# LFA-TAA-002 
S-DWG# LFA-TAA-003 
DWG# LFA-TAA-004 
S- DWG# LFA-TAA-002 
S-DWG# LFA-TAA-003 
DWG# LFA-TAA-004 
S-DWG# LFA-TAA-002 
^ — DWG# LFA-TAA-003 
’—DWG# LFA-TAA-004 
S-DWG# LFA-TAA-002 
I s — DWG# LFA-TAA-003 
V '— DWG# LFA-TAA-004 



COMPLERE INC. 


NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 
TRAVERSE SYSTEM CABLE ROUTING DIAGRAM 


DRAWING DATE 

JULY 8, 1991 


FILENAME 

LANGLEY TRAVERSE DIAGRAM 


DESIGN ENGINEER 


DRAWING NUMBER 


TODD A. AMBUR 


LFA-TAA-001 




COMPLERE INC. 
MOTOR DRIVE SYSTEM 


COMPUMOTOR 
STEPPER MOTOR 


SOCKET SIGNAL PIN SIGNAL PIN SIGNAL SOCKET SIGNAL 



Lh 

O 


AMPHENOL CONNECTOR AMPHENOL CONNECTOR 
2060434 2060444 

AMPHENOL SOCKETS AMPHENOL CABLE CLAMP 
66360-2 206070-1 

AMPHENOL PINS 
66361-2 


BELDEN CABLE 
9418 


AMPHENOL CONNECTOR 
206044-1 

AMPHENOL CABLE CLAMP 
206070-1 

AMPHENOL PINS 
66361-2 


AMPHENOL CONNECTOR 
206043-3 

AMPHENOL CABLE CLAMP 
206070-1 

AMPHENOL SOCKETS 
66360-2 


SIGNAL 

DESCRIPTION 

A+ 

Motor Winding 

A- 

Motor Winding 

B+ 

Motor Winding 

B- 

Motor Winding 

SHIELD 

Motor Case Ground 

m 

No Connection 


COMPLERE INC. 

NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 

MOTOR DRIVE SYSTEM TO COMPUMOTOR STEPPER MOTOR 

DRAWING DATE 

JULY 3, 1991 

FILENAME 

LANGLEY MOTOR 

DESIGN ENGINEER 

TODD A. AMBUR 

DRAWING NUMBER 

LFA-TAA-002 





COMPLERE INC 
MOTOR DRIVE SYSTEM 



AMPHENOL CONNECTOR AMPHENOL CONNECTOR 
206705-1 206708-1 

AMPHENOL PINS AMPHENOL CABLE CLAMP 

66103-2 206966-1 

AMPHENOL SOCKETS 
66105-2 


SIGNAL 

DESCRIPTION 

A+ 

Quadrature Encoder Signal 

A- 

Logical Complement of A+ 

B+ 

Ouadrature Encoder Signal 

B- 

Logical Complement of B+ 

Z+ 

Once Per Revolution 

~T- ~ 1 

Logical Complement of Z+ 

SHIELD ™ 

Case Ground 

GROUND 1 

Logical Ground 

NC 

No Connection 


DYNAMICS RESEARCH 
ENCODER 


SOCKET SIGNAL 


PIN SIGNAL 




BELDEN CABLE 
9504 


AMPHENOL CONNECTOR 
206708-1 

AMPHENOL CABLE CLAMP 
206966-1 

AMPHENOL SOCKETS 
66105-2 




AMPHENOL CONNECTOR 
206705-2 

AMPHENOL CABLE CLAMP 
206966-1 

AMPHENOL PINS 
66103-2 


COMPLERE INC. 


NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 
MOTOR DRIVE SYSTEM TOTCS8 ENCODER SIGNALS 


DRAWING DATE 


FILENAME 


JULY 3, 1991 


LANGLEY ENCODER 


DESIGN ENGINEER 


DRAWING NUMBER 


TODD A. AMBUR 


LFA-TAA-003 

















COMPLEREINC 
MOTOR DRIVE SYSTEM 


PIN SIGNAL SOCKET SIGNAL 



AMPHENOL CONNECTOR AMPHENOL CONNECTOR 

206061-1 206060-1 

AMPHENOL PINS AMPHENOL CABLE CLAMP 
66103-2 206062-1 

AMPHENOL SOCKETS 
66105-2 


"sTGNaC 

DESCRIPTION 

Home 

Home Switch Signal 

cw 

End of Travel Limit Signal 
Clockwise 


TCW 

End of Travel Limit Signal 
Counter Clockwise 


comM6N 

Logical Ground 


LINEAR INDUSTRIES 
LIMIT SWITCHES 


SOCKET SIGNAL PIN SIGNAL 



BELDEN CABLE 
9418 


AMPHENOL CONNECTOR 
206060-1 

AMPHENOL CABLE CLAMP 
206062-1 

AMPHENOL SOCKETS 
66105-2 


AMPHENOL CONNECTOR 
206153-2 

AMPHENOL CABLE CLAMP 
206062-1 

AMPHENOL PINS 
66103-2 


COMPLERE INC . 

NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 

MOTOR DRIVE SYSTEM TO LINEAR INDUSTRIES LIMIT SWITCHES 

DRAWING DATE 

JULY 3, 1991 

FILENAME 

LANGLEY LIMIT SWITCH 

DESIGN ENGINEER 

TODD A. AMBUR 

DRAWING NUMBER 

LFA-TAA-004 


- 




COMPLEREINC. 
MOTOR DRIVE SYSTEM 

SOCKET SIGNAL 


SIGNAL 


SOCKET SIGNAL 



10 I C4GND 


11 


12 I NC 


13 1 NC 


14 C1A+ 


15 Cl GND 


16 C2 +5V 


17 I C2B+ 


18 I C2 SHIELD 


20 I C3 GND 


21 


22 C4B+ 


23 I C4 SHIELD 


24 NC 


DB-25 FEMALE 




C1B+ 


BEimBBli 


DB-25 MALE 


6 FT. DB-25 STRAIGHT THROUGH 
EXTENSION CABLE 


C2A+ 


C2GND 


C3+5V 


C3 B+ 


8 C3 SHIELD 


C4A+ 


10 C4 GND 


11 


12 NC 


13 NC 


14 Cl A+ 


15 Cl GND 


16 C2+5V 


17 C2B+ 


18 C2 SHIELD 


19 I C3A+ 


20 C3 GND 


C4+5V 


22 C4B+ 


23 I C4 SHIELD 


24 NC 


DB-25 FEMALE 


SIGNAL 


NC 


Cl B+ 



10 I C4GND 


11 


12 NC 


13 NC 


14 C1A+ 


15 Cl G 


16 NC 


17 


18 C2 SHIELD 


19 1 C3A+ 


20 C3GND 


21 INC 


22 C4B+ 


23 C4 SHIELD 


24 I NC 


DB-25 MALE 


Cl- CHANNEL 1 
C2- CHANNEL 2 
C3- CHANNEL 3 
C4- CHANNEL 4 


COMPLERE INC. 

NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 
MDS TO TCS8 ENCODER SIGNALS 


DRAWING DATE 


FILENAME 


JULY 3, 1991 


ENCODER SIGNALS 


DESIGN ENGINEER 


DRAWING NUMBER 


TODD A. AMBUR 


LFA-TAA-005 





























































































































































HP 9000-330 


SOCKET SIGNAL PIN SIGNAL 



DB-9 FEMALE DB-9 MALE 


Ln 

4 * 


SIGNAL 

DESCRIPTION 

TxD 1 

Transmit Data 

RxD 

Receive Data 

“RTS 

Ready to Send 

TT3 “1 

Clear to Send 

DCD 

Data Carrier Detect 

“DTK 1 

Data Terminal Ready 

USE 

Data Set Ready 

" 

Ring Indicator 

GND 

Logical Ground 

' 'NC 

No Connection 

































































TCS8 

MDS1/MDS2 


SOCKET SIGNAL 


1 

TxD 

2 

RTS 

3 

GND 

4 

NC 

5 

NC 

6 

RxD 

7 

CTS 

8 

NC 

9 

NC 


-<f-. 

-«r- 

-«r~ 


-<f-: 


PIN 

SIGNAL 

1 

TxD 

2 

NC 

3 

GND 

4 

NC 

5 

NC 

6 

RxD 

7 

NC 

8 

NC 

9 

NC 




DB-9 FEMALE 


DB-9 MALE 


Ln 


SIGNAL 

DRSPR TPTTOM 

TxD 

Transmit Data 

RxD 

Receive Data . _ 

RTS 

Ready to Send 

CTS 

Clear to Send 

GND 

T ogical Ground 

_NC 

Nn Connection 


MDS 

HOST 






BELDEN CABLE 9925 


PIN 


SIGNAL 


SOCKET SIGNAL 


1 

TxD 

2 

NC 

3 

GND 

HEXB 

NC 


NC 

6 

RxD 

7 

NC 

8 

NC 

9 

NC 







1 

RxD 

2 

NC 

3 

GND 

4 

NC 

5 

NC 

6 

TxD 

7 

NC 

8 

NC 

9 

NC 


DB-9 MALE 


DB-9 FEMALE 


COMPLERE INC. 

NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 
TCS8 TO MDS SERIAL 

DRAWING DATE 

JULY 3, 1991 

FILENAME 

TCS8 SERIAL 

DESIGN ENGINEER 

TODD A. AMBUR 

DRAWING NUMBER 

LFA-TAA-007 






































































MDS 

CHAIN 

SOCKET SIGNAL 


1 

RxD 

2 

NC 

3 

GND 

4 

NC 
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NC 

6 

TxD 
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NC 

8 

NC 

9 

NC 


HH 


PIN 

SIGNAL 

1 

RxD 

2 

NC 

3 

GND 

4 

NC 

5 

NC 

6 

TxD 

7 

NC 

8 

NC 

9 

NC 


PIN 


SIGNAL 


MDS 

HOST 

SOCKET SIGNAL 


z 







1 

TxD 

2 

NC 

3 

GND 

4 

NC 

5 

NC 

6 

RxD 

7 

NC 

8 

NC 

9 

NC 



1 

RxD 

2 

NC 

3 

GND 

4 ...... 

NC 

5 

NC 

6 

TxD 

7 

NC 

8 

NC 

9 

NC 


DB-9 FEMALE 


DB-9 MALE 


BELDEN CABLE 9925 


DB-9 MALE 


DB-9 FEMALE 


Ln 

o\ 


SIGNAL 

DESCRIPTION 

TxD 

Transmit Data 

RxD 

Receive Data 

GND 

I logical Ground 

_fcj£ 

NnConnectinn 


COMPLERE INC . 

NASA LaRC LASER FLUORESCENCE ANEMOMETER 
LFA TRAVERSE SYSTEM 
MDS TO MDS SERIAL 

DRAWING DATE 

JULY 3, 1991 

FILENAME 

CHAIN SERIAL 

DESIGN ENGINEER ' 

TODD A. AMBUR 

DRAWING NUMBER 

LFA-TAA-008 





























Appendix B Laser Velocimeter Data Acquisition System 

The LVDAS acquires simultaneous digital data, analog data, and time information data. 

The data are sampled, multiplexed, buffered, and then transferred to the facility’s host computer 
for further data reduction, analysis, and presentation. 

Four 16 bit word parallel input ports are provided to accept the digital output of LV counter 
processors and/or other instrumentation. 

New applications in laser velocimetry have brought about the need for a more advanced 
laser velocimeter data acquisition system. These new applications require high data rates that are 
not hindered by on-line time dependent data sorting and real time graphic data presentation. The 
new Laser Velocimeter Data Acquisition System (LVDAS) was designed specifically to meet these 
advanced requirements. 

High data acquisition rates are achieved by providing a separate latched input for each laser 
velocimeter digital input and a separate converter for each laser fluorescence analog input. The 
system will allow for a data acquisition rates of approximately 100,000 samples per second 
simultaneously on each of the laser velocimeter and laser fluorescence inputs. 

A 32 bit time of day (TOD) 10MHz counter is used to tag arrival times to acquired digital 
LDV data as they become available on each of four digital inputs. When a data valid sync pulse is 
sensed for a particular channel, the LVDAS latches the current TOD into a 32 bit time of arrival 
register (TOA). A separate TOA register is available for each digital input so that particle arrival 
times of measured velocity information U,V, W can be monitored for coincidence. The latched 
times of arrivals have a resolution of 100 ns and maximum time of over 7 minutes. 

All of the acquired digital velocity data with corresponding time of arrival data can be 
processed and stored. However, if coincident data is required, then the arrival time of the various 
channels can be conditionally accepted if they all occur within a finite window of time. These 
coincident events can then be assigned interarrival times which represent elapsed time since the 
previous event. 

The coincidence control logic allows for 3 channel coincidence. The coincidence time is 
adjustable from 0. 1 |is to 1 s. In addition to the laser velocimeter inputs, three additional data 
words are generated internally. They are the inter arrival time, the coincidence time, and status 
words. The inter arrival and coincidence time is provided by a clock whose resolution and 
maximum time is 100 ns and 500 seconds respectively. The status word contains information 
about coincidence and the order in which the laser velocimeter data arrived. 

During data acquisition, it is important that the user obtain some visual feedback about the 
data being acquired. This is necessary so that the user can make informed decisions about both the 
quality and quantity of data received. The user is either reassured about the quality of the data or 
can make alterations and improvements in technique while on line. To help achieve this, the 
instantaneous velocities are used to generate real time histograms from which probability density 
distributions are determined for all velocity components. 

Additionally, the laser velocimeter data acquisition system has the capability of reducing the 
raw laser velocimeter data. Each laser velocimeter output contains the information required to 
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calculate the instantaneous velocities. From the instantaneous velocity determinations, the average 
velocities, turbulence levels, and the turbulence cross correlations are all be calculated. 

The coincidence control logic will allow for up to 4 channel coincidence. The coincidence 
time can be adjustable to any resolution or duration within the capability of the time of arrival 
registers. When coincident criteria are met, the analog inputs can be sampled and converted to 
provide concurrent data with the digital data. A single time of arrival is latched for all four of the 
analog to digital inputs since they are all sampled and converted simultaneously. A final time of 
arrival is latched for external events. These might be derived from such sources as oscillating 
models or model surfaces, rotating helicopter blades, rotating engine fans, or flow sensors. 

All digital Macrodyne data, optional digital data, analog to digital data, and time of arrival 
data can be sent by the LVDAS to other computers via two serial and two parallel input/output 
ports. One parallel port will be used for the HP series 9000 model 330 computer while the other 
can be used by the facility host computer. The serial ports can be used by PC type computers such 
as IBMs or MACs. Software has been developed for on-line data acquisition and display. A 
program listing is enclosed. 
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LDVWT 


100 Main 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550. 

560 
570 
580 
590 
600 
610 
620 . 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 


NASA LANGLEY RESEARCH CENTER Property of COMPLERE INC. 

16 BY 24 INCH WATER TUNNEL Proprietary software 

Copyright February 25, 1992 

LASER FLUORESCENCE ANEMOMETER Developed by: T. Kevin McDevitt 


! PROGRAM DESCRIPTION: 

! This program provides the capability to acquire simultaneous Laser Doppler Velocimeter (LDV) , Laser 

! Fluorescence Anemometer (LFA) , and Analog Voltage Data at user selectable traverse controlled probe volume 

! positions within the water ttnnel flow, 

! The LVDAS (Laser Velocimeter Data Acquisition System) is used to sample the LDV, LFA, and Analog Voltage 

! data simultaneously with a coincidence criterion being applied to LDV incoming data. The LVDAS also generates 

! interarrival times and coincidence time. 

! The measured LDV data provides the necessary frequency information from which three components of flow 

! velocities can be determined. These velocities are measured directly in "LASER” coordinates. Coordinate 

I system transformations are applied to these measured velocities to obtain velocities in "TUNNEL" and "MODEL" 

! coordinates. 

! The TCS8 (Traverse Control System) is used to precisely move the LDV probe volume within the tunnel and 

! about the model. The TCS8 provides three axes plus one auxiliary axis of traverse capability for both the 

! transmitting and receiving side optical packages. The Tx and Rx side traverses can be moved independently to 

! achieve laser alignment or they can be moved together to maintain laser alignment. 

! The TCS8 will give the traverse positions in TCS8 coordinates where one inch of commanded movement will 

! yield one inch of movement on the traverse slides. However, this will not yield one inch of movement of the 

! probe volume crossover point within the water filled tunnel because of the differences of refraction in air, 

\ glass, and water. Therefore, coordinate system transformations are applied to TCS8 positions to obtain 

I positions in "TUNNEL" and "MODEL" coordinates. 

! During data acquisition, real time histograms will be displayed of the LDV and analog data. After the data 

! has been acquired, the averages, standard deviations, and shear stresses will be calculated and displayed in 

! profile plots where the data is plotted versus traverse position. The reduced data is also sent to the 

! printer in tabular form. The reduced data as well as the raw data are stored along with the tunnel conditions 

! on the hard disc for archival purposes and also to allow for further data reduction, data plotting, or data 

! transfer to other computers. 


! PROGRAM OPERATION: 

! The following power up sequences should be completed before this program is run: 

! 1. Turn on the "Motor Drive System" boxes. 

! 2. Turn on the "TCS8" traverse control system, 

! 3. Turn on the "LVDAS" Laser Velocimeter Data Acquisition System. 

! 4. Turn on the HP series 9000 model 330 computer. 

! This program will automatically be loaded and run when the computer is turned on. If it is not loaded then 

! you can type in the following commands to load and then run it. 

! LOAD "LDVWT:, 1400, 0,0" 

! RUN 

! When the program is ready for user operation, it will display three things on the CRT. These are the main 

! menu, TCS8 traverse positions, and new sets of histogram & profile graphs. If they do not appear on the CRT 

! then you should perform the following actions to reinitialize the systems. 

! 1. Press shift reset on the HP series 9000 model 330 computers keyboard. 

! 2. Press reset on the back of the TCS8 . 

! 3. Press reset on the front (or back) of the LVDAS. 

! 4. LOAD "LDVWT:, 1400,0,0" 

! 5. RUN 


! PROGRAM VARIABLES: 

! 

I Mass Storage Variables: 

! 

! System$ Tells the program where to read/store system data related files. 

! Data$ Tells the program where to read/store raw & reduced data related files. 

! File$ File name for tunnel conditions data or raw & reduced data. 


! Menu Variables: 
I 


! 

j 

! 

t 

! 


MenuS (*) 
Menu 

Key 

Busy 

Ready 


String array where each element describes its corresponding menu subroutine^ function. 
Used as an index to the string array MenuS (*) . Indicates which of the menus has been 
selected as the current menu. 

Used as an index to the string array MenuS (*) . Indicates which one of eight menu 
subroutines in the menu is to be executed. 

Tells the Menu Status subprogram to display the current menu selection in inverse video. 
Tells the Menu Status subprogram to display the current menu selection in normal text. 


Traverse Position Variables: 


Tcsl(*) TCS8 transmitting side traverse positions (XI, Yl, Zl, Al) in TCS coordinates. 
Tcs2(*) TCS8 receiving side traverse positions (X2, Y2, Z2,A2) in TCS coordinates. 
Tun ( *) Traverse positions (X,Y,Z) in TUNNEL coordinates. 
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Mod ( * ) 


Traverse positions (X,Y,Z) in MODEL coordinates. 


900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
12 90 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
14 50 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 


Auto Move Traverse Position Variables: 

Pos(*) Array of preprogrammed auto move positions. 

Pname${*) Names for the variables in Pos(*). 

Pimage$(*) Image formats for the variables in Pos(*). 

Punits$(*) Units for the variables in Pos(*). 

Npos Number of preprogrammed auto move positions in Pos(*). 

Paxis Specifies which axis is to be traversed for the profile. Also defines axis for plots. 

Traverse Coordinate System Transformation Variables: 

Index(*) Array of indexes of refraction for air, glass, and water. 

Index (1) : Index of refraction for Air. 

Index (2) : Index of refraction for Glass. 

Index (3) : Index of refraction for Water. 

Theta Tx Side Off Axis Angle. 

Fs Focal length for sending side onaxis and offaxis lenses. 

Fr Focal length for receiving side offaxis lens. 

Bs Beam spacings for .sending side onaxis and offaxis beam pairs. 

Br Beam spacing for receiving side offaxis. 

Ts Angle of offaxis sending side beam pair. 

Tr Angle of offaxis receiving side beam pair. 

Ta Sending side offaxis auxiliary traverse angle. 

Tcs2tunl(*) Sending side coordinate system transformation matrix for converting Tcsl(*) to Tun{*) . 

Tun2tcsl(*) Sending side coordinate system transformation matrix for converting Tun(*) to Tcsl(*). 

Tcs2tun2(*) Receiving side coordinate system transformation matrix for converting Tcs2{*) to Tun(*). 

Tun2tcs2(*) Receiving side coordinate system transformation matrix for converting Tun(*) to Tcs2(*). 

Tun2mod(*) Coordinate system transformation matrix for converting Tun(*) to Mod(*) . 

Mod2tun(*) Coordinate system transformation matrix for converting Mod(*) to Tun(*) . 

Velocity Coordinate System Transformation Variables: 

Index(*) Array of indexes of refraction for air, glass, and water. 

Index(l) : Index of refraction for Air. 

Index (2) : Index of refraction for Glass. 

Index(3) : Index of refraction for Water. 

Thetal (*) Angles between LASER & TUNNEL UVW laser beams in Air (N=Indexl) . 

Theta3 (*) Angles between LASER & TUNNEL UVW laser beams in Water (N=Index3) . 

Tun21dv(*) Coordinate system transformation matrix for converting from TUNNEL to LASER. 

Ldv2tun(*) Coordinate system transformation matrix for converting from LASER to TUNNEL. 

Traverse 4 Velocity Coordinate System Transformation Variables: 

,Alpha(*) Angles of attack, yaw, and roll. 

Alpha (1) : Angle of Attack. 

Alpha (2) : Angle of Yaw. 

Alpha (3) : Angle of Roll. 

Mod2tun(*) Coordinate system transformation matrix for converting positions & velocities from MODEL to TUNNEL, 

Tun2mod(*) Coordinate system transformation matrix for converting positions 4 velocities from TUNNEL to MODEL. 

Tunnel Condition Variables: 

Array(*) Array of tunnel conditions, laser parameters, graph scales, etc. 

Name$(*) Names for the variables in Array(*). 

Image$ (*) Image formats for the variables in Array (*>. 

Units$ ( *) Units for the variables in Array (*). 

Misc. Tunnel Condition Variables: 

Date Date. 

• Time Time. 

Run Run Number. 

File File Number. 

Mach Mach Number. 

Temp Room Temperature (deg. F) . 

Uedge Freestream Velocity. 

Ujet_ue Jet exit velocity normalized by Uedge. 

LVDAS Variables: 

Table(*) Lookup table of frequencies. 

Atime The maximum desired acquisition time (seconds). 

Ctime The maximum desired coincidence time (seconds) . 

At_exp Exponent for interarrival times. 

Ct_exp Exponent for coincidence times. 

Nreads Number of desired samples. 

Nsam Number of acquired samples. 

Coin(*) Coincidence criteria. 
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Cmask 
Raw (*) 


Coincidence mask for U,V, W selection. 

Array of raw data acquired from the LVDAS. 


1700 
1710 
1720 
1730 
1740 
1750 
* 1760 

1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
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2110 
2120 
2130 
2140 
2150 
21 60' 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
22 90 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
r 2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 


Instantaneous Velocity and Voltage Variables: 

Ui(*) Read from LVDAS as the instantaneous U frequency data# then converted into U velocities. 

Vi<*) Read from LVDAS as the instantaneous V frequency data# then converted into V velocities. 

Wi(M Read from LVDAS as the instantaneous W frequency data# then converted into W velocities. 

Ai(*) Read from LVDAS as the instantaneous A voltage data. 

Bi{*> Read from LVDAS as the instantaneous B voltage data. 

Ii<*) Read from LVDAS as the raw interarrival time data# then converted into interarrival times 

Ci ( *) Read from LVDAS as the raw coincidence time data, then converted into coincidence times 

Valid(*) Validation words. Initially all ones, then some set to zero during histogram clipping. 

Histogram Clipping Variables: 

Umin The minimum acceptable U frequency (MHz) . 

Umax The maximum acceptable U frequency (MHz) . 

Vmin The minimum acceptable V frequency (MHz) . 

Vmax The maximum acceptable V frequency (MHz). 

Wmin The minimum acceptable W frequency (MHz) . 

Wmax The maximum acceptable W frequency (MHz) . 

Clip Clip: 1 turn histogram clipping on; 0 turns it off. 

Frequency to Velocity Conversion Variables: 

Beam_spc(*) Beam spacing at lens. 

Focl_len ( * ) Focal length. 

Bearn_sep(*) Beam separation angle in degrees (full angle). 

Wave_len(*) Wave length. 

Frng_spc(*) Fringe Spacings. 

Brg frq(*) Bragg Frequencies. 

Mix_frq(*) Mixing Frequencies. 

Mea;_sgn(*) Measured Frequencies' Signs. 

Brg_sgn(M Bragg Frequencies' Signs. 

Mix_sgn(*) Mixing Frequencies' Signs. 

Summation Variables: 

Sum(l#l) Summation of all of the valid Ui. 

Sum(2#l) Summation of all of the valid Vi. 

Sum(3#l) Summation of all of the valid Wi. 

Sum(4#l) Summation of all of the valid Ai . 

Sum(5,l) Summation of all of the valid Bi. 

Sum(6#l) Summation of all of the valid Ii. 

Sum (7# 1) Summation of all of the valid Ci. 

Sum(l#2) Summation of all of the valid Ui*Ui. 

Sum(2#2) Summation of all of the valid Vi*Vv. 

Sum(3#2) Summation of all of the valid Wi*Ww. 

Sum(4#2) Summation of all of the valid Ai*Ai. 

Sum ( 5# 2) Summation of all of the valid Bi*Bi. 

Sum(6,2) Summation of all of the valid Ii*Ii. 

Sum(7,2) Summation of all of the valid Ci*Ci. 

Sum(l# 3) Summation of all of the valid Ui*Vi. 

Sum (2# 3) Summation of all of the valid Vi*Wi. 

Sum (3,3) Summation of all of the valid Wi*Ui. 

Sum (4,3) Summation of all of the valid Ai*Bi. 

Sum ( 5# 3) Summation of all of the valid Ui*Ai. 

Sum (6, 3) Summation of all of the valid Vi*Ai. 

Sum(7,3) Summation of all of the valid Wi*Ai. 

N ( * ) Number of valid samples for the above summations. 

Reduced Data Variables: 

U Average U frequency or velocity. 

V Average V frequency or velocity. 

W Average W frequency or velocity. 

A Average A voltage. 

B Average B voltage. 

I Average interarrival time. 

C Average coincidence time. 

Ul Standard deviation for U frequency or velocity. 

VI Standard deviation for V frequency or velocity. 

Wl Standard deviation for W frequency or velocity. 

Al Standard deviation for A voltage. 

Bl Standard deviation for B voltage. 

II Standard deviation for interarrival time. 

Cl Standard deviation for coincidence time. 

Ulvl Velocity: Velocity Shear Stress. 

Vlwl Velocity: Velocity Shear Stress. 
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Wlul 

Velocity: Velocity 

Shear 

Stress. 

Albl 

Voltage :Voltage 

Shear 

Stress. 

Ulal 

Velocity: Volt age 

Shear 

Stress. 

Vial 

Velocity: Voltage 

Shear 

Stress . 

Wlal 

Velocity : Voltage 

Shear 

Stress. 


! Data Plotting Symbol Variables: 

! 

] Symbols (*) Array of Symbol arrays. Each symbol arrays contains a distinct geometric symbol. 

! Symbol (*) Array of goordinates which when connected produce a distinct geometric symbol. 

! Dot<*) Array of coordinates which produce a dot. The dot symbol is added to all symbols. 

! Noc The number of coordinates in a symbol. 

! S Used to index the Symbols array. 

! Histogram and Profile Graph Variables: 

Array containing the plot’s scales. 

Array containing the plot's CRT position. 

Array containing the number of X divisions for the plot’s X axis. 

Array containing the number of Y divisions for the plot's Y -axis. 

String array containing labels for the X axis. 

String array containing labels for the Y axis. 

String array containing labels for the Plots. 

String array containing image formats for the X axis labeling. 

String array containing image formats for the Y axis labeling. 

String array containing labels for each symbol in a profile plot. 

Used as an index to the above arrays. Specifies one of nine plots. 

Used to save the entire graphics contents of the CRT. 

OPTION BASE 1 

COM /Data/ INTEGER Raw <1000, 10) , Valid (1000) , REAL Table (0 : 32766) , Ui (1000) , Vi (1000) , Wi < 1000 ) ,Ai (1000) , 
Bi (1000) , Ii (1000) , Ci (1000) 

COM /Array/ Name$ (100,4) [10] , Image$ (100,4) [10], Unit s$ (100, 4) [10], REAL Array<100,4) 

COM /Pos/ P name$ (25,1) [10 ] , P image $ (25,1) [ 10 ] , PunitsS (25,1) (10 ] , REAL Pos (25, 1) , Npos 
COM /Graph/ Wndw { 9, 4) , Vwprt (9, 4) , Xdiv (9) , Ydiv (9) # XlabelS ( 9) [80) , YlabelS (9) [80 ] , TitleS ( 9) [80], 

XimageS (9) [80 ] , YimageS ( 9) [ 80 ] , Legend$ < 9, 5) [80] 

COM Run, File, Paxis 

DIM MenuS (5,8) [80 ], SystemS [20 ], DataS [ 20 ], File$ [ 50 ], L$ [160 ] 

INTEGER Gsave ( 1280, 1024) , At_exp, Ct_exp, Cmask, Nsam, N (10, 3) 

REAL Atime, Ct ime, Sum (10, 3) , Symbols (5,0:20,3), Apos, Bpos 

DIM Tcs2tunl (4, 4) ,Tun2tcsl (4, 4) , Tun2mod (3 , 3) , Tun21dv(3,3) ,Tun<4) , Tcsl (4) 

DIM Tcs2t un2 (4,4) ,Tun2tcs2(4, 4) ,Mod2tun (3, 3) , Ldv2tun (3, 3) ,Mod<4) , Tcs2 (4) 

DIM Beam_spc ( 3) , Focl__len ( 3) , Mea_sgn (3) , Mix_frq (3) , Mix__sgn (3) , Frng_spc (3) , Thetal (3,3) 

DIM Beam_sep ( 3) , Wave_len ( 3) , Brg__f rq (3) , Brg_sgn (3) , Index (3) , Coin (3 ) , Theta 3 (3,3) , Alpha (3) 
l Perform trigonometric operations in degrees. 

DEG 


! Wndw (* ) 

! Vwprt (*) 

! Xdiv {*) 

! Ydiv (*) 

! XlabelS ( *) 

! YlabelS ( *) 

! TitleS < * ) 

! XimageS (M 

! YimageS (*) 

! LegendS (*) 

! G 

! Gsave(*) 


Here: 


On_key : 


Keys: 


! Clear the CRT and direct printed output to it. 


CLEAR SCREEN 
GCLEAR 

PRINTER IS CRT 

! Perform any necessary setup and initialization routines. 


GOSUB Lvds_set_up 
GOSUB File_set_up 
GOSUB Tcs8_set_up 
GOSUB Menu_set_up 
GOSUB Grph_set_up 


! Initialize the HP to LVDAS interface. 

! Select mass storage device for system & data files. 

! Initialize the HP to TCS8 interface. 

! Initialize the user driven menus and display the main menu. 

I Initialize the CRT and plot the nine empty plots for profiles and histograms. 

! The main program, while continually displaying the time of day, will wait hear for menu key selection. 

Date=TIMEDATE 

Time=Date 


DISP TIMES (Time) , DATES (Date) 
GOTO Here 


ON 

KEY 

1 

GOSUB 

Keyl 

! If 

the 

user 

f unct ion 

key 

#1 

is 

ever 

pressed 

then 

execute 

the 

’’Keyl" 

subroutine 

ON 

KEY 

2 

GOSUB 

Key2 

i If 

the 

user 

function 

key 

#2 

is 

ever 

pressed 

then 

execute 

the 

"Key 2" 

subroutine 

ON 

KEY 

3 

GOSUB 

Key3 

! If 

the 

user 

function 

key 

#3 

is 

ever 

pressed 

then 

execute 

the 

"Key 3" 

subroutine 

ON 

KEY 

4 

GOSUB 

Key4 

! If 

the 

user 

funct ion 

key 

#4 

is 

ever 

pressed 

then 

execute 

the 

"Key 4" 

subroutine 

ON 

KEY 

5 

GOSUB 

Key5 

! If 

the 

user 

function 

key 

#5 

is 

ever 

pressed 

then 

execute 

the 

"Key 5" 

subroutine 

ON 

KEY 

6 

GOSUB 

Key6 

! If 

the 

user 

function 

key 

#6 

is 

ever 

pressed 

then 

execute 

the 

"Key 6" 

subroutine 

ON 

KEY 

7 

GOSUB Key7 

! If 

the 

user 

function 

key 

#7 

is 

ever 

pressed 

then 

execute 

the 

"Key 7" 

subroutine 

ON 

KEY 

8 

GOSUB 

Key8 

! If 

the 

user 

function 

key 

#8 

is 

ever 

pressed 

then 

execute 

the 

"Key 8" 

subroutine 


RETURN 

! Subroutine Keyl, Key2, Key3, Key4 , Key 5, Key6, Key7, Key8 descriptions: 

! When one of the special user function keys is pressed, the main program will execute one the 

! following eight subroutines. Each of these subroutines performs essentially the same basic 

! function in that it subsequently executes one of the menu subroutines. The particular menu 

! subroutine to be executed will depend on the current menu selected and the current key pressed. 

! Before the selected menu subroutine is executed, the corresponding menu entry at the top of 

! the CRT is redisplayed in inverse video. This indicates that the menu selection has been 

! acknowledged and that any resultant actions are still in progress. When the highlighted menu 

! subroutine has completed the current TCS8 traverse positions will be read and updated on the CRT 

! display. The corresponding menu entry displayed at the top of the CRT is redisplayed in normal 
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3590 
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3650 
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3670 

3680 

3690 

3700 

3710 

3720 Key7: 

3730 

3740 

3750 

3760 

3770 

3780 Key8 : 
3790 
3800 . 

3810 

3820 

3830 

3840 Menuls 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 Mlkl : 
4070 


! text to indicate the completion of the menu subroutine. The user can then select another special 

! function key. 

! Variables: 


! 

; 

! 

! 


Menu 

Key 

Menu$ ( *) 

Busy 

Ready 


Indicates which of the menus has been selected as the current menu. 

Indicates which one of eight menu subroutines in the menu is to be executed. 

String array where each element describes its corresponding menu subroutine’s function. 
Tells the Menu Status subroutine to display the current menu selection in inverse video. 
Tells the Menu Status subroutine to display the current menu selection in normal text. 


Key=l 

CALL Menu_status (Menu, Key, Busy, Menu$ (*) ) 

ON Menu GOSUB Mlkl,M2kl,M3kl, M4kl,M5kl,M6kl,M7kl 
CALL Menu_status (Menu, Key, Ready, Menu$ (*) ) 

CALL Tcs8read (@Tcs8, Mod (*) , Tun (*) ,Tcsl(*) ,Tcs2<*) ,Tcs2tunl<*> ,Tcs2tun2(*) ,Tun2mod<*) ) 


RETURN 


Key=2 

CALL Menu_status (Menu, Key, Busy, Menu$ (*) ) 

ON Menu GOSUB Mlk2,M2k2,M3k2,M4k2,M5k2,M6k2,M7k2 
CALL Menu_status (Menu, Key, Ready, Menu$ (*) ) 

CALL Tcs8read (@Tcs8,Mod(*) , Tun(*) ,Tcsl (*) ,Tcs2 (*) ,Tcs2tunl (*) ,Tcs2tun2 (*) ,Tun2mod (*) ) 

RETURN 

Key-3 

CALL Menu_status (Menu, Key, Busy, Menu$ ( *) ) 

ON Menu GOSUB Mlk3,M2k3,M3k3, M4k3,M5k3,M6k3,M7k3 
CALL Menu_status (Menu, Key, Ready, Menu$ (*) ) 

CALL Tcs8read ( @Tcs8, Mod ( * ) , Tun (*) , Tcsl (*) ,Tcs2 (*) ,Tcs2tunl (*) ,Tcs2tun2 (*) ,Tun2mod(*) ) 

RETURN 

Key=4 

CALL Menu_stat us (Menu, Key, Busy, MenuS ( *) ) 

ON Menu GOSUB Mlk 4 , M2k4 , M3k4 , M4k4 , M5k4 , M6k4 , M7k4 
CALL Menu_status (Menu, Key, Ready, Menu$ ( *) ) 

CALL Tcs8read (@Tcs8, Mod (*) , Tun(*) ,Tcsl (*) , Tcs2 ( *> , Tcs2tunl (*) , Tcs2tun2 (*) ,Tun2mod (*) ) 
RETURN 


Key=5 

CALL Menu_status (Menu, Key, Busy,Menu$ ( *) ) 

ON Menu GOSUB Mlk5,M2k5,M3k5,M4k5,M5k5,M6k5,M7k5 
CALL Menu^status (Menu, Key, Ready, Menu$ (*) ) 

CALL Tcs8read (@Tcs8, Mod (*) ,Tun(*) ,Tcsl (*) ,Tcs2(*) ,Tcs2tunl(*) ,Tcs2tun2(*) ,Tun2mod(*) ) 

RETURN 

Key=6 

CALL Menu_status (Menu, Key, Busy,Menu$ (*) } 

ON Menu GOSUB Mlk 6, M2k6, M3k6, M4k6, M5k 6, M6k6, M7k6 
CALL Menu_status (Menu, Key, Ready, Menu$ (*) } 

CALL Tcs8read (@Tcs8, Mod ( * ) ,Tun(*> ,Tcsl (*) ,Tcs2<*) ,Tcs2tunl (*) ,Tcs2tun2 (*) ,Tun2mod (*) ) 

RETURN 

Key=7 

CALL Menu_status (Menu, Key , Busy, MenuS (*) ) 

ON Menu GOSUB Mlk7, M2k7, M3k7, M4k7 , M5k7, M6k7, M7k7 
CALL Menu_status (Menu, Key, Ready, MenuS (*) ) 

CALL Tcs8read (0Tcs8, Mod (* ) , Tun ( *) , Tcsl ( *) , Tcs2 (*) , Tcs2tunl ( *) , Tcs2tun2 (*) , Tun2mod (*) ) 
RETURN 


Key=8 

CALL Menu_status (Menu, Key, Busy, MenuS ( *) ) 

ON Menu GOSUB Mlk8,M2k8,M3k8,M4k8,M5k8,M6k8,M7k8 
CALL Menu_status (Menu, Key, Ready, MenuS (*) ) 

CALL Tcs8read (@Tcs8, Mod (* ) , Tun (*) , Tcsl (*) ,Tcs2(*) ,Tcs2tunl(*) , Tcs2tun2 ( * ) , Tun2mod ( *) ) 

RETURN 

! Descriptions of the "Main Menu" subroutines Ml K1 , . . . , M1K8 : 

! The eight subroutines M1K1 , . . . , M1K8 together implement the "Main Menu". The following will be 

! displayed at the top left of the CRT display when the "Main Menu" is selected: 

i 

1 M1K1: Laser Alignment 

! M1K2 : Pre Run 

! M1K3: Post Run (Dump Graphics) 

! M1K4: Set Auto Move Positions 

! M1K5: Move traverse 

! M1K6: Take data 

! M1K7: Auto move and take 

! M1K8; Display Histograms 

! 

! Ml K1 will change the current active menu from the "Main Menu" to the "Laser Alignment Menu". 

! M1K2 will change the current active menu from the "Main Menu" to the "Pre Run Menu". M1K3 will 

I transfer the graphics contents of the CRT to the printer. This provides a hard copy of the profile 

! plots. M1K4 has the user enter predefined traverse positions for a profile plot. M1K5 moves the 

! traverse to a user selectable position. M1K6 acquires LVDAS data at the current TCS8 traverse 

! position. M1K7 acquires LVDAS data at each of the pre programed TCS8 traverse positions set up by 

! Ml K4 . M1K8 repeatedly displays five channels of real time histograms until the user presses any 

! key on the keyboard. 

i 

! Change the current active menu from the "Main Menu" to the "Laser Alignment Menu". 

Menu-2 
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4090 
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4850 
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CALL Menu__disp (Menu,Menu$ (*) > 

RETURN 

! Change the current active menu from the "Main Menu** to the "Pre Run Menu". 


Menu=3 


CALL Menu_disp (Menu, Menu$ (*) ) 

RETURN 

! Transfer the graphics contents of 

KEY LABELS OFF 

PRINTER IS CRT/WIDTH 132 

DISP "" 

FOR L=1 TO 9 

P RINT TABXY ( 1 , L) ,• RPT$ (" " , 120 ) 
NEXT L 

PRINTER IS PRT 
PRINT USING 
DUMP GRAPHICS 
PRINT USING "#,Q" 

PRINTER IS CRT 

CALL Menu_disp (Menu, Menu$ ( * ) ) 


the CRT to the printer. This provides a hard copy of the plots. 

! Turn off the key labels so that they won’t be printed. 

! Clear the CRT's display line so that they won't be printed. 
! Clear the CRT's menu lines so that it won't be printed. 


Move to the top of the next page on the printer. 
Dump the entire CRT's contents to the printer. 
Move to the top of the next page on the printer. 

Redisplay the menus. 


RETURN 

! Have the user enter predefined traverse positions for a profile plot. 
CALL Enter_value ("number of traverse positions", Npos, "K") 

REDIM Pos (Npos, 1 ) , PnameS (Npos, 1) , Pimage$ (Npos, 1) , Punits$ (Npos, 1) 

MAT Pimage$= ("M4D.4D") 

MAT Punit s$= ("in") 


FOR K=1 TO Npos 

Pname$ (K, 1) ="PosU"&VAL$ (K) 
NEXT K 


GSTORE Gsave ( * ) 

CALL Change ("VALUES", Pos (*) , Pname$ (*) , Pimage$ (*) , Punits$ (*) ) 
GLOAD Gsave (*) 

CALL Menu_disp (Menu, Menu$ ( * ) ) 

RETURN 


! Moves the traverse to a user selectable position. 

GOSUB Read_calc_fill 

CALL Tcs8read (0Tcs8, Mod ( * ) , Tun ( *) ,Tcsl (*) ,Tcs2(*) ,Tcs2tunl (*) ,Tcs2tun2(*) ,Tun2mod(*) ) 

CALL Ente revalue (CHR$ (NUM ("X") +Paxis-1) ,Mod(Paxis) , "K") 

ON KBD CALL Do_nothing 
DISP "Moving" 

Movement =Mod ( Paxi s) 

CALL Tcs8move (@Tcs8, Mod ( * ) , Tun (*) ,Tcsl (*) ,Tcs2(*) ,Mod2tun(*) , Tun2tcsl <*) , Tun2tcs2 (*) , "Tx & Rx", "MODEL", 
"ABSOLUTE", Pax is, Movement ) 

CALL Tcs8read (@Tcs8,Mod(M ,Tun(*) ,Tcsl (*) ,Tcs2 (*) ,Tcs2tunl (*) ,Tcs2tun2(*) ,Tun2mod(*> ) 

GOSUB Calc 
GOSUB Fill 
DISP "" 

OFF KBD 
RETURN 

! Acquire LVDAS data at the current TCS8 traverse position. 

! DISP "Press any key to TAKE DATA" 

! CALL Rt_histo (@Lvdas, Symbols (*), 1 ) 

Cmask=Coin ( 1) *l+Coin(2) *2+Coin(3) *4 
Nsam=MIN (Nreads, 1000) 

Date=TIMEDATE 

Time=Date 

CALL Lvdas_take (QLvdas, At ime, Ctime, At_exp, Ct_exp, Cmask, Nsam) 

IF Nsam>l THEN 

SELECT Paxis ! Select the non-scanning axes for printing their position values at each point 

CASE 3 

D$="X" 

F$="Y" 

Apos=Mod (1) 

Bpos=Mod (2) 

CASE 2 

D$="X" 

F$="Z" 

Apos=Mod(l) 

Bpos=Mod (3) 

CASE 1 

D$="Y" 

F$="2" 

Apos=Mod (2) 

Bpos=Mod (3) 

END SELECT 

OUTPUT PRT USING "K, K";CHR$ (27) &"&k2S"&CHR$ (27) &"&19D», RPT$ ("=", 140) 

PRINTER IS PRT 
Run$=VAL$ (Run) 

File$=VAL$ (File) 

PRINT " RUN " &Run$&" FILE "&File$ 

A$=DATE$ (TIMEDATE) ! Acquire the date and time for printing at each point. 
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B$-TIME$ (TIMEDATE) 

PRINT USING 4890;A$,B$,D$,Apos,F$, Bpos 

IMAGE 6X, 11A, 3X, 8A, 3X, A, 3D.4D, 3X, A, "=", 3D.4D 

PRINTER IS CRT 

CALL Data_reduce (At_exp, Ct_exp, Nsam) 

CALL Pt_histo (Symbols (*) , Run, File, Mod (Paxis) , Nsam) 

CALL Data_clip (Nsam, Umin, Umax, Vmin, Vmax, Wmin, Wmax) 

CALL Data_sum (Sum (*), N (*), Nsam) 

CALL Data_calc(N(*) , Sum(* ) , U, V, W, A, B, 10, CO, U1 , VI, Wl, Al, B1 , II, Cl, Ulvl, Vlwl, Wlul, Albl, Ulal, Vial, Wlal) 
B=U/U ! Replacement of B, Bl, and Albl is being made by velocity ratios since 

B1=V/U ! the second analog channel B is not being used 

Albl=W/U 

CALL Data_print (Pax is, Mod (Pax is) , Nsam, "MHz" , U, V, W, A, B, 10, C0,U1, V1,W1, Al, Bl, II, Cl, Ulvl, Vlwl, Wlul, Albl, 
Ulal, Vial, Wlal, Uedge) 

CALL Data_fconvert (Array (*) ) 

CALL Data_sum (Sum (*), N (*), Nsam) 

CALL Data_calc <N(*),Sum(*) ,U,V,W,A, B, 10, CO, Ul, VI, Wl, Al, Bl, II, Cl , Ulvl , Vlwl , Wlul, Albl, Ulal, Vial , Wlal) 

B=U/Uedge 

Bl=V/Uedge 

Albl=W/Uedge 

CALL Data_print ( Paxis, Mod (Paxis) , Nsam, "LDV", U, V, W, A, B, 10, CO, Ul, VI , Wl, Al, Bl, II, Cl, Ulvl, Vlwl, Wlul, Albl, 
Ulal, Vial, Wlal, Uedge) 

CALL Data_trnsfrm(Ldv2tun(*) ,U,V,W, Ul, VI, Wl, Ulvl, Vlwl, Wlul, Ulal, Vial, Wlal) 

B=U/Uedge 

Bl=V/Uedge 

Albl=W/Uedge 

CALL Datajorint (Paxis, Mod (Paxis) ,Nsam, "TUN" , U, V, W, A, B, 10, CO , Ul , VI , Wl, Al, Bl, II , Cl, Ulvl , Vlwl, Wlul, Albl, 
Ulal, Vial, Wlal, Uedge) 

CALL Data_trnsfrm (Tun2mod (*) , U, V, W, Ul , VI , Wl , Ulvl , Vlwl , Wlul, Ulal , Vial , Wlal ) 

B=U/Uedge 

Bl=V/Uedge 

Albl=W/Uedge 

CALL Data_print (Paxis, Mod (Paxis) , Nsam, "MOD", U, V, W,A, B, 10, CO, Ul, VI, Wl , Al , Bl , II , Cl , Ulvl , Vlwl, Wlul, Albl, 
Ulal, Vial, Wlal, Uedge) 

CALL Data_plot ( Array (*) , Symbols (*) , 6, Mod ( Paxis) , U, V, W, 1/Uedge, N (1, 1) , N(2, 1) ,N(3,1) ) 

CALL Datajplot (Array (*), Symbols (*), 7, Mod (Paxis) ,U1,V1,W1, l/Uedge,N(l, 2) , N (2, 2) , N ( 3, 2) ) 

CALL Datajplot ( Array (*) , Symbols (* ) ,8, Mod ( Paxis) , Ulvl, Vlwl , Wlul, 1/ Uedge* 2, N (1, 3) , N (2, 3 ) , N (3, 3) ) 

CALL Data_plot (Array (M , Symbols (*) , 9, Mod (Paxis) , A, A, Al, 1, N (4, 1) , N (4, 1) , N (4, 2) ) 

OUTPUT PRT USING "K, K";CHR$ (27) &"&k2S"&CHR$ (27) &"&19D", RPT$ (" =", 140) 

GOSUB Store_file 
File=File+l 
END IF 
RETURN 

! Acquire LVDAS data at each of the pre programed TCS8 traverse positions set up by M1K4. 

Quit=0 

ON KBD GOSUB Quit 
FOR J=1 TO Npos 

CALL Tcs8read (@Tcs8, Mod ( * ) ,Tun(*) ,Tcsl (*) ,Tcs2 (*) ,Tcs2tunl (*) ,Tcs2tun2(*) ,Tun2mod<*) ) 

Mod (Paxis) =Pos (J, 1) 

GOSUB Mlk5a 
GOSUB Mlk6 
IF Quit THEN 5360 
NEXT J 
OFF KBD 
GOSUB On_key 

CALL Menu__disp (Menu,Menu$ (*) ) 

RETURN 

! Repeatedly displays five channels of real time histograms until the user presses any key on the keyboard. 
DISP "Press any key to return to main menu" 

CALL Rt_histo (QLvdas, Symbols ( *) , 1 ) 

RETURN 

! Descriptions of the "Laser Alignment Menu" subroutines M2K1, . . . , M2K8 : 

! The eight subroutines M2K1, . . . , M2K8 together implement the "Laser Alignment Menu". The 

! following will be displayed at the top left of the CRT display when the "Laser Alignment Menu" is 

! selected: 


! 

! 


M2K1 will change the current active menu from the "Laser Alignment Menu" to the "Main Menu". 
M2K2 selects whether the transmitting, receiving, or both sides of the traverse are to be moved. 
M2K3 selects the TCS, TUNNEL, or MODEL coordinate systems for traverse movements. M2K4 
specifies movements to be relative to the currents position or to absolute positions. M2K5 has the 
user enter a movement for the X axis and then the movement is performed. M2K6 has the user enter 


M2K1 : 

Return to main menu 

M2K2 : 

Sides : 

Tx & Rx 

M2K3: 

Coordinates: 

MODEL 

M2K4: 

Mode : 

ABSOLUTE 

M2K5 : 

Move X 


M2K6: 

Move Y 


M2K7 : 

Move Z 


M2K8: 

Move A 
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5960 

5970 
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6000 

6010 

6020 M2k5: 
6030 M2k6: 
6040 M2k7 : 
6050 M2k8 : 
6060 
6070 
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6090 
6100 
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6130 
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6180 Menu3: 
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6240 
6250 
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6300 
6310 
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6340 
6350 
6360 
6370 
6380 
63 90 

6400 M3kl: 
6410 


! a movement for the Y axis and then the movement is performed. M2K7 has the user enter a movement 

! for the Z axis and then the movement is performed. M2K8 has the user enter a movement for the A 

! axis and then the movement is performed. 

! 

! Change the current active menu from the "Laser Alignment Menu" to the "Main Menu". 

Menu=l 

CALL Menu_disp (Menu,Menu$ (*) ) 

RETURN 

! Select whether the transmitting, receiving, or both sides of the traverse are to be moved. 

SELECT TRIMS (Menu$ (Menu, Key) [20] > 

CASE "Tx & Rx" 

MenuS (Menu, Key) [20]="Tx" 

CASE "Tx" 

Menu$ (Menu, Key) [20]="Rx" 

CASE "Rx" 

Menu$ (Menu, Key) [201 ="Tx & Rx" 

END SELECT 

CALL Menu_disp (Menu, Menu$ (*) ) 

RETURN 

! Selects the TCS, TUNNEL, or MODEL coordinate systems for traverse movements. 

SELECT TRIMS (MenuS (Menu, Key) [20]) 

CASE "MODEL" 

MenuS (Menu, Key) [20] ="TUNNEL" 

CASE "TUNNEL" 

MenuS (Menu, Key) [20]="TCS" 

CASE "TCS" 

MenuS (Menu, Key) [20]="MODEL" 

END SELECT 

CALL Menu^disp (Menu, MenuS (*) ) 

RETURN 

! Specifies movements to be relative to the currents position or to absolute positions. 

SELECT TRIMS (MenuS (Menu, Key) [20]) 

CASE "ABSOLUTE" 

MenuS (Menu, Key) [20] ="RELATIVE" 

CASE "RELATIVE" 

MenuS (Menu, Key) [20 ] ^"ABSOLUTE" 

END SELECT 


CALL Menu_disp (Menu, MenuS (*) ) 

RETURN 

! The subroutines M2K5 thru M2K8 all execute the same code. The code will have the user enter a 

! movement for the X,Y,Z, or A depending on what the value of "Key" is. The user specified movement 

! for the selected axis will then be performed. 

t 

SideS =TRIM$ (MenuS (Menu, 2) [20]) 

Coor$=TRIM$ (MenuS (Menu, 3) [20] ) 

Mode$=TRIM$ (MenuS (Menu, 4) [20] ) 

CALL Enter_value (ModeSfi" Movement ", Movement, "4D.5D") 

ON KBD CALL Do_nothing 
DISP "Moving" 

CALL Tcs8read (@Tcs8,Mod (*) , Tun (*) ,Tcsl(*) ,Tcs2(*) , Tcs2t uni (* ) , Tcs2tun2 ( * ) ,Tun2mod(*)) 

CALL Tcs8move ( @Tcs8, Mod ( * ) , Tun (*) , Tcsl (*) , Tcs2 (*) , Mod2tun (*) , Tun2tcsl ( *) , Tun2tcs2 (*) , SideS, CoorS, ModeS, 
Key-4 , Movement ) 

CALL Tcs8 read (@Tcs8, Mod (* ) , Tun (* ) , Tcsl (*) , Tcs2 ( *) , Tcs2tunl (*) , Tcs2tun2 (*) , Tun2mod (*) ) 


DISP "» 
OFF KBD 
RETURN 


! 

t 

! 

t 

! 

I 

j 

! 

! 


Descriptions of the "Pre Run Menu" subroutines M3K1 , . . . , M3K8 : 

The eight subroutines M3K1 , . . . , M3K8 together implement the "Pre Run Menu". The following will 
be displayed at the top left of the CRT display when the "Pre Run Menu" is selected: 

M3K1: Return to MAIN menu 

M3K2: Enter Run & File Numbers 

M3K3: Enter Number of Samples 

M3K4; Select Traverse Axis for Profile 

M3K5: Print Coordinate Transformation Matrices 

M3K6: Setup Graphics 

M3K7: Tunnel Conditions 

M3K8: Traverse 


! M3K1 will change the current active menu from the "Pre Run Menu" to the "Main Menu". M3K2 has 

! the user enter a the Run and File numbers. A new run number should be assigned to each profile 

! while a new file number is assigned to each set of data. M3K3 has the user enter the desired 

! number of samples. M3K4 has the user select which axis to traverse in for the profiles. M3K5 

! prints the coordinate system transformation matrices for both traverse positions and velocities. 

1 M3K6 creates a new set of empty plots for new profiles. M3K7 will change the current active menu 

! from the "Pre Run Menu" to the "Tunnel Conditions Menu". M3K8 will change the current active menu 

! from the "Pre Run Menu" to the "Traverse Menu". 


! Change the current active menu from the "Pre Run Menu" to the "Main Menu". 
Menu=l 
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CALL Menu_disp (Menu, Menu? (*) ) 

RETURN 

! Have the user enter a the Run and File numbers. 

CALL Enter_value ( "Run", Run, "3D. 2D") 

CALL Enter_value ( "File*’ , File, "3D" ) 

RETURN 

! Have the user enter the desired number of samples. 

CALL Enter_value ( "Number of Samples “ , Nreads, "K" ) 

RETURN 

! Have the user select which axis to traverse in for the profiles. 

CALL Enter_string ("Traverse Axis for Profiles ", PaxisS, “K") 

SELECT Paxis$ 

CASE "X" 

Paxis=l 
CASE M Y“ 

Paxis=2 
CASE “2" 

Paxis=3 
CASE "A" 

Paxis=4 
CASE ELSE 

GOTO M3k4 
END SELECT 
GOSUB Fill 
RETURN 

! Prints the coordinate system transformation matrices for both traverse positions and velocities. 

GOSUB Read_calc_f ill 
OUTPUT PRT USING “#,2/" 

OUTPUT PRT USING "20X, K, /“; “TRAVERSE COORDINATE TRANSFORMATION MATRICES" 

OUTPUT PRT USING "20X, K, / , 4 ( 1 3X, 4 (8D . 5D) , / ) " ; “Transmitting side TCS8 to TUNNEL" # Tcs2t uni ( *> 

OUTPUT PRT USING "20X, K, /, 4 (13X, 4 (8D. 5D) “Receiving side TCS8 to TUNNEL" , Tcs2tun2 ( * ) 

OUTPUT PRT USING "20X, K, / , 4 ( 1 3X, 4 (8D. 5D) , /) "; “Transmitting side TUNNEL to TCS8" , Tun2tcsl ( * ) 

OUTPUT PRT USING "20X, K, /, 4 (13X, 4 (8D. 5D) "Receiving side TUNNEL to TCS8", Tun2tcs2 (*) 

OUTPUT PRT USING "20X, K, /, 3 (13X, 3 (8D. 5D) , /) "; "TUNNEL to MODEL ", Tun2mod (*) 

OUTPUT PRT USING “20X, K, / , 3 (13X, 3 (8D. 5D) , /) "; “MODEL to TUNNEL" , Mod2tun ( *) 

OUTPUT PRT USING "20X, K, / “ ; "VELOCITY COORDINATE TRANSFORMATION MATRICES" 

OUTPUT PRT USING “20X, K, / , 3 (13X, 3 (8D. 5D) , /} “LASER to TUNNEL", Ldv2tun <*) 

OUTPUT PRT USING "20X, K, / , 3 ( 1 3X, 3 (8D . 5D) , "TUNNEL to LASER", Tun21dv<*) 

OUTPUT PRT USING “20X, K, /, 3 (13X, 3 (8D. 5D) , /) "; "TUNNEL to MODEL", Tun2mod ( *) 

OUTPUT PRT USING “20X, K, / , 3 (13X, 3 (8D. 5D) , /) “MODEL to TUNNEL", Mod2tun < *) 

OUTPUT PRT USING 
RETURN 

! Display a new set of plots for new profiles. 

CALL Setup_graph (Array (*) , ImageS (*}, Paxis, Symbols {*) ) 

RETURN . 

! Change the current active menu from the "Pre Run Menu" to the "Tunnel Conditions Menu". 

Menu»4 

CALL Menu_disp (Menu, Menu$ (*) ) 

RETURN 

! Change the current active menu from the "Pre Run Menu" to the “Traverse Menu". 

Menu=5 

CALL Menu_disp (Menu, MenuS (*) ) 

RETURN 

! Descriptions of the "Tunnel Conditions Menu" subroutines M4K1 , . . . , M4K8 : 

! The eight subroutines M4K1 , . . . , M4K8 together implement the "Tunnel Conditions Menu". The 

! following will be displayed at the top left of the CRT display when the "Tunnel Conditions Menu" is 

! selected: 

t 

! M4K1: Return to PRE RUN menu 

! M4K2: Load Tunnel Conditions 

! M4K3: Save Tunnel Conditions 

! M4K4 : Print Tunnel Conditions 

! M4K5: Enter Tunnel Condition Data 

! M4K6: Enter Tunnel Condition Names 

! M4K7 : Enter Tunnel Condition Units 

! M4K8: Enter Tunnel Condition Images 

! 

! M4K1 will change the current active menu from the "Tunnel Conditions Menu" to the "Pre Run 

1 Menu". M4K2 loads the old tunnel conditions from a file on the disk. M4K3 saves the current 

! tunnel conditions to a file on the disk, M4K2 & M4K3 load and save default tunnel conditions from 

! the file "ARRAY" on the hard disk. The default values are not related to any particular run number. 

! M4K4 sends the current tunnel conditions to the printer. M4K5 has the user enter values for the 

! tunnel condition variables. M4K6 has the user enter names for the tunnel condition variables. 

! M4K7 has the user enter units for the tunnel condition variables. M4K8 has the user enter image 

! formats for the tunnel condition variables. 

! Change the current active menu from the "Tunnel Conditions Menu" to the "Pre Run Menu". 

Menu=3 

CALL Menu_disp (Menu, Menu$ (*) ) 

RETURN 
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! Load the old tunnel conditions from a file on the disk. This loads the default values. 

GOSUB Read_array 
GOSUB Read_calc_fill 
RETURN 

! Save the current tunnel conditions to a file on the disk. This updates the default values on the disk. 
GOSUB Read_calc_f ill 
GOSUB Save_array 
RETURN 

! Print the current tunnel conditions. 

GOSUB Read_calc_f ill 
GOSUB Print_header 
RETURN 

! Have the user enter values for the tunnel condition variables. 

GSTORE Gsave ( *) 

GOSUB Read_calc_f ill 

CALL Change ( "VALUES" , Array ( * ) , Name? (*) , Image$ ( *) , Units? (*) ) 

GOSUB Read_calc_f ill 
GLOAD Gsave (*) 

RETURN 

! Have the user enter names for the tunnel condition variables. 

GSTORE Gsave (*) 

GOSUB Read_calc_fill 

CALL Change ( "NAMES**, Array ( *) , Name$ (*} , Image$ ( *) , Units? (*) > 

GOSUB Rea d__ca 1 c_f ill 
GLOAD Gsave (*■) 

RETURN 

l Have the user enter units for the tunnel condition variables. 

GSTORE Gsave (*) 

GOSUB Rea d_ca 1 c_f ill 

CALL Change ( "UNITS'*, Array ( * ) , Name? (*) , Image? (*), Units? ( *) ) 

GOSUB Rea d_ca 1 c_f ill 
GLOAD Gsave (*) 

RETURN 

! Have the user enter image formats for the tunnel condition variables. 

GSTORE Gsave (*) 

GOSUB Read_ca lc_f il 1 

CALL Change ( " IMAGES", Array ( *) , Name? ( * ) , Image? (*) , Unit s? ( * ) ) 

GOSUB Rea d_ca 1 c_f ill 
GLOAD Gsave (*) 

RETURN 

! Descriptions of the "Traverse Menu" subroutines M5K1 , . . . , M5K8 : 

! The eight subroutines M5K1 , . . . , M5K8 together implement the "Traverse Menu". The following will 

! be displayed at the top left of the CRT display when the "Traverse Menu" is selected: 

i 

! M5K1; Return to PRE RUN menu 

! M5K2: View & Set TCS8 Positions 

! M5K3 : View & Set TCS8 Units 

! M5K4 : View & Set TCS8 Revolution 

! M5K5 : View & Set TCS8 Velocity • 

! M5K6; View & Set TCS8 Acceleration 

! M5K7: 

! M5K8 : 

! M5K1 will change the current active menu from the "Traverse Menu" to the "Pre Run Menu". M5K2 

! reads from the TCS8 the current positions and lets the user change them. The new positions are 

! then sent to the TCS8 . M5K3 reads from the TCS8 the current counts per unit length (inches) and 

! lets the user change them. The new counts per unit length are then sent to the TCS8. M5K4 reads 

! from the TCS8 the current counts per revolution and lets the user change them. The new counts per 

! revolution are then sent to the TCS8 . M5K5 reads from the TCS8 the current velocities and lets the 

l user change them. The new velocities are then sent to the TCS8 . M5K6 reads from the TCS8 the 

! current accelerations and lets the user change them. The new accelerations are then sent to the 

! TCS8. M5K7 does nothing. M5K8 does nothing. 

i 

I Change the current active menu from the "Traverse Menu" to the "Pre Run Menu". 

Menu=3 


CALL Menu_disp (Menu, Menu? (*) ) 
RETURN 

! Read current TCS8 positions, 
CALL Tcs8set ( "P", @Tcs8) 
GRAPHICS ON 

CALL Menu_disp (Menu, Menu? (*) ) 
RETURN 

! Read current TCS8 counts per 
CALL Tcs8set ("U", @Tcs8) 
GRAPHICS ON 

CALL Menu_disp (Menu, Menu? (*) ) 
RETURN 

! Read current TCS8 counts per 
CALL Tcs8set ( "R", QTcs8) 
GRAPHICS ON 


have the user update them, & then send the new values to the TCS8. 
! View and set TCS8 Positions. 


inch, have the user update them, & then send the new values to the TCS8. 
! View and set TCS8 counts per Unit length. 


revolution, have the user update them, & then send new values to the TCS8. 
! View and set TCS8 counts per Revolution. 
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CALL Menu_disp (Menu, Menu$ (*) ) 

RETURN 

! Read current TCS8 velocities, have the user update them,. & then send the new values to the TCS8. 

CALL Tcs8set ("V'*, @Tcs8) ! View and set TCS8 Velocities. 

GRAPHICS ON 

CALL Menu_disp(Menu,Menu$ (*) ) 

RETURN 

! Read current TCS8 accelerations, have the user update them, & then send the new values to the TCS8. 

CALL Tcs8set ( "A", @Tcs8) ( View and set TCS8 Accelerations. 

GRAPHICS ON 

CALL Menu__disp (Menu, MenuS ( * ) ) 

RETURN 

RETURN ! This subroutine does nothing. 

RETURN ! This subroutine does nothing. 

Quit=l ! Quit will be set during a multiple position scan (see M1K7) if any key is pressed on the 

RETURN ! keyboard during the scan. This indicates that the scan should be terminated. 

! This subroutine initializes the HP to LVDAS high speed parallel interface. A communications path named 
! "@Lvdas" is opened. Also, this subroutine creates the raw data to frequency conversion look up table. 
CALL Lvdas__init (QLvdas) 

CALL Table (Table (*) ) 

RETURN 

! This subroutine reads the initialization files from the disk. System$ tells the program where to read 
! system related files while Data$ tells the program where to read and store raw and reduced data. 

System$=" : , 1400, 0, 0" 

Data$=":, 1400,0,1" 

LOAD KEY " KEY S"& SystemS 
GOSUB Read_array 
GOSUB Rea d_ca 1 c_f ill 
GOSUB Save_ar ray 
CLEAR SCREEN 
RETURN 

I This subroutine initialized the HP to TCS8 serial interface. The communications path "@Tcs8" is opened. 
CALL Tcs8 init (@Tcs8) 

CALL Tcs8read (@Tcs8, Mod (* ) ,Tun(*) ,Tcsl(*) ,Tcs2(*) ,Tcs2tunl(*) ,Tcs2tun2(*) ,Tun2mod(M) 

GOSUB Calc 
GOSUB Fill 
RETURN 

! This subroutine defines the graphics symbols for plotting data points, clears and initializes the CRT, 

! and displays a new empty set of graphs for histogram and profile plotting. 

CALL Read symbols (Symbols (*) ) 

CALL Crt_init 

CALL Setup_graph (Array (*) , Image$ ( *) , Paxis, Symbols (*) ) 

RETURN 

l This subroutine defines the menu descriptors for all of the menus. The current menu is set to the "Main 
! Menu" and its menu is displayed at the top of the screen. 

CALL Menu_read (Menu$ ( *) ) 

CALL Menu_disp (Menu,Menu$ (*> ) 

GOSUB On key 
Busy=0 
Ready=l 
RETURN 

! This subroutine prints a header on the printers paper. The "header" is a formatted list of all of the 
! tunnel conditions, laser parameters, and graph scales. 

PRINTER IS PRT; WIDTH 144 

PRINT USING »#,5(K)";CHR$(27) &"&k2S"&CHR$ (27) &"&19D" 

CALL Arrayjprint (Array (*), Name$ (*), Image$ (*), Unit s$ (*) ) 

PRINT USING 5(K) ";CHR$ (27) &"E" 

PRINTER IS CRT 
RETURN 

! This subroutine extracts (reads) the tunnel conditions from the Array (*) . These values can be used to 
! calculate other tunnel conditions. The original tunnel conditions along with any calculated tunnel 
! conditions are then put back (filled) into the Array (*). 

GOSUB Read 
GOSUB Calc 
GOSUB Fill 
RETURN 

! This subroutine stores the header Array (*) and other arrays onto the disk. There will be one header 
! file for each run number. For example, if the run number equal 1, then the data will be stored in a 
! file named "Rl". This file will include an extensive list of tunnel conditions, laser parameters, graph 
! scales, traverse positions, coordinate system transformation matrices, etc. 

DISP "Storing Header" 

! Set File$ equal to the file name for the header file. Each run number will have a different file name. 
Fi le$="R" &VAL$ (Run) &Data$ 

! Check if the file already exists. If it does then, ask the user if he wants to overwrite the old file. 
ON ERROR GOTO 8960 
ASSIGN @Data TO File$ 

OFF ERROR 
FOR K«1 TO 10 
WAIT .2 
BEEP 
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8820 
8830 
8840 
88 50 
8860 
8870 
8880 
8890 
8900 
8910 
8920 
8930 
8940 
8950 
8960 
8970 
8980 
8990 
9000 
9010 
9020 
9030 
9040 

9050 Store_f lie: 

9060 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

9160 

9170 

9180 

9190 

9200 

9210 

9220 

9230 

9240 

92 50 

9260 

9270 

9280 

9290 

9300 

9310 

9320 

9330 

9340 . 

9350 

9360 

9370 

9380 

9390 

9400 

9410 Read_array: 

9420 

9430 

9440 

9450 

94 60 

9470 

9480 

94 90 

9500 

9510 

9520 

9530 

9540 

9550 

9560 

9570 

9580 

9590 

9600 

9610 


NEXT K 

INPUT "Over Write old file? (Y or N) ",L$ 

SELECT L$(l,l] 

CASE "Y", "y" ! If the user wants to overwrite the old file, then purge the old file. 

ASSIGN GData TO * 

PURGE File$ 

GOTO 8960 

CASE "N","n" ! If the user doesn't want to overwrite the old file, then have a new run# entered. 

CALL Enter_value ( "Run" , Run, "3D. 2D") 

CALL Enter_value ( "File", File, "3D") 

GOTO Store_header 
CASE ELSE 

GOTO Store_header 
END SELECT 
OFF ERROR 

Fsize^INT ( (3200+4000*3+128*4+72*4) /25 6*1 .05+1) ! 

CREATE BDAT File$,Fsize ! 

ASSIGN GData TO File$ ! 

OUTPUT GData; Array (*), Name$ ( * ) , ImageS ( *) , Unit s$ (* ) 

OUTPUT GData; Tun2tcsl (*) ,Tun2tcs2 (*) ,Mod2tun(*) ,Tun21dv(*) 

OUTPUT GData; Tcs2t uni (*) , Tcs2tun2 (*) , Tun2mod ( *) , Ldv2t un (*) 

ASSIGN GData TO * ! 

RETURN 

! This subroutine stores the header Array (*), the raw data, and the reduced data onto the disk. There 
! will be one data file for each data set. For example, if the run and file numbers equal 7 and 5 
! respectively, then the data will be stored in a file named "R7F5". 

GOSUB Calc ! Use the tunnel conditions to calculate and/or update other tunnel conditions. 

GOSUB Fill ! Fill Array(*) with the original tunnel conditions along with the updated tunnel conditions 

IF File=l THEN GOSUB Store_header 
DISP "Storing Data" 

File$="R"&VAL$ (Run) &"F"&VAL$ (File) &Data$ 


Calculate the headers file size. 
Create the header's file. 

Open the header's file. 


Close the header's file. 


! Check if the file already exists. If it does, then ask the user if he wants to overwrite the old file. 


ON ERROR GOTO 9340 
ASSIGN GData TO File$ 

OFF ERROR 
FOR K=1 TO 10 
WAIT .2 
BEEP 
NEXT K 

INPUT "Over Write old file? (Y or N) »,L$ 

SELECT L$ [ 1 , 1 ] 

CASE "Y", "y" ! If the user wants to overwrite the old file, then purge the old file. 

ASSIGN GData TO * 

PURGE File$ 

GOTO 9340 

CASE "N", ,, n" ! If the user doesn't want to overwrite the old file, then have a new run# entered. 

CALL Ente revalue ( "Run", Run, "3D. 2D") 

CALL Enter^value ("File", File, "3D") 


GOTO Store_file 
CASE ELSE 

GOTO Store_file 
END SELECT 
OFF ERROR 

Fsi ze=INT ( (3200+Nsam* 10*2+60+240) /256*1.05+1) 
CREATE BDAT File$,Fsize 
ASSIGN GData TO File$ 


! Calculate the data's file size. 
! Create the data's file. 

! Open the data's file. 


OUTPUT GData; Array (*), Raw (*) ,N(*) ,Sum(*) 

ASSIGN GData TO * ! Close the data's file. 

RETURN 

! This subroutine reads the header Array (*) off of the disk from a file named "ARRAY". The file will the 
! have default values for the tunnel conditions, laser parameters, graph scales, etc. This file is not 
! meant to be attached to any run number or profile scan. It is used to provide default values for the 
1 program so that the user will not have to enter a rather lengthy list of tunnel conditions. 

ON ERROR GOTO 9550 

! If the file already exists, then read the Array(*) from the disk. 

ASSIGN GFile TO " ARRAY"&$ystem$ 

ENTER QFi le; Array (*) , Name$ ( *) , Image$ (*) ,Units$(*) 

ENTER GFile;Tun2tcsl ( *) , Tun2tcs2 (*) ,Mod2tun (*) , Tun21dv (*) 

ENTER GFile;Tcs2tunl ( *) ,Tcs2tun2 ( *) , Tun2mod (*) ,Ldv2tun<*) 


ASSIGN GFile TO * 

OFF ERROR 
RETURN 

! If the file doesn't exist then create the file, read in default data, and store the Array(*) on disk. 

OFF ERROR 

ASSIGN GFile TO * 

ON ERROR GOTO 9590 
PURGE "ARRAY" fi System$ 

OFF ERROR 


CALL Array_init. (Name$ (*) , Array(*) , ImageS (*) ,Units$ (*) ) 
CREATE BDAT "ARRAY"iSystem$, 50 
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9620 

9630 

9640 Save_array: 

9650 

9660 

9670 

9680 

9690 

9700 

9710 

9720 

9730 

9740 Fills 


10220 

10230 

10240 

10250 

10260 

10270 

10280 

10290 

10300 

10310 

10320 

10330 

10340 

10350 

10360 

10370 

10380 

10390 

10400 

10410 


GOSUB Save_array 
RETURN 

! This subroutine saves the header Array (*) onto the disk in a file named "ARRAY". The file will then have 
! default values for the tunnel conditions# laser parameters, graph scales, etc. This file is not meant to 
! be attached to any run number or profile scan. It is used to provide default values for the program so 
! that the user will not have to enter a rather lengthy list of tunnel conditions. 

ASSIGN @ File TO " ARRAY" &System$ 

OUTPUT @File;Array (*) , Name$ ( * ) , Image$ ( *) , Unit s$ <*> 

OUTPUT @File;Tun2tcsl (*) , Tun2tcs2 ( *) , Mod2tun ( *) , Tun21dv ( *) 

OUTPUT @File;Tcs2tunl (*) , Tcs2tun2 (*) # Tun2mod ( *) # Ldv2tun(*> 

ASSIGN @File TO * 

RETURN 

! This subroutine fills the Array (*) with the current tunnel conditions, laser parameters, and histogram 


9750 

! & profile scales. 


9760 

Array (1, 1) -Date 

! Date. 

9770 

Array (1,2) -Mach 

! Mach Number. 

9780 

Array (1, 4) -Alpha ( 1) 

! Angle of Attack. 

9790 

Array (2, 1 ) =Time 

! Time. 

9800 

Array (2,2) =Temp 

! Room Temperature (deg. F) . 

9810 

Array (2, 4) =Alpha (2) 

! Angle of Yaw. 

9820 

Array (3, 1 ) -Run 

! Run Number. 

9830 

Array (3, 2) =Uedge 

! Freestream Velocity. 

9840 

Array (3, 4) =Alpha(3) 

! Angle of Roll. 

9850 

Array (4, l)=File 

! File Number. 

9860 

Array (4, 2) =U jet_ue 

! Jet exit velocity normalized by Uedge. 

9870 

Array (4, 4) =Theta 

! Tx Side Off Axis Angle. 

9880 

MAT Array (11:14,1)- Mod 

! Probe volume positions in MODEL coordinates. 

9890 

MAT Array (11:14,2)* Tun 

I Probe volume positions in TUNNEL coordinates. 

9900 

MAT Array (11:14, 3) = Tcsl 

! Tx side traverse positions in Tcs8 coordinates. 

9910 

MAT Array (11:14, 4)* Tcs2 

! Rx side traverse positions in Tcs8 coordinates. 

9920 

MAT Array (21, 1:3) = Index 

! Index of refraction of for laser light (eg: Nair,Nglass, Nwater) . 

9930 

MAT Array (22:24, 1:3)= Thetal 

I Angles between LASER & TUNNEL UVW laser beams in Air (N=Indexl) . 

9940 

MAT Array (25:27,1:3)= Theta3 

! Angles between LASER & TUNNEL UVW laser beams in Water (N=Index3) 

9950 

MAT Array (31, 1 : 3) = Beam spc 

! Beam spacing at lens. 

9960 

MAT Array (32, 1:3) = Focl_len 

! Focal length. 

9970 

MAT Array (33, 1: 3) = Beam_sep 

! Beam separation angle in degrees (full angle). 

9980 

MAT Array (34 , 1 : 3) = Wave len 

I Wave length. 

9990 

MAT Array (35, 1 : 3) = Frng_spc 

! Fringe spacing. 

10000 

MAT Array (36, 1 : 3} = Brg frq 

! Bragg frequency. 

10010 

MAT Array (37, 1:3) = Mix_frq 

! Mixing frequency. 

10020 

MAT Array (38, 1 : 3) = Mea_sgn 

! Sign of measured frequency in velocity equation. 

10030 

MAT Array (39, 1 : 3} = Brg_sgn 

! Sign of bragg frequency in velocity equation. 

10040 

MAT Array (40, 1 : 3) = Mix_sgn 

! Sign of mixing frequency in velocity equation. 

10050 

MAT Array (41, 1 ; 3} = Coin 

! Coincidence criteria 

10060 

Array (42, 1) =Umin 

! Frequency minimum for U calculation. 

10070 

Array (42, 2) =Vmin 

! Frequency minimum for V calculation. 

10080 

Array (42, 3} =Wmin 

! Frequency minimum for W calculation. 

10090 

Array (43,1) =Umax 

! Frequency maximum for U calculation. 

10100 

Array (43,2) =Vmax 

! Frequency maximum for V calculation. 

10110 

Array (43,3) =Wmax 

! Frequency maximum for W calculation. 

10120 

Array (51, 1) =Nreads 

1 Number of desired samples. 

10130 

Array (52, 1) =Nsam 

! Number of acquired samples. 

10140 . 

Array (51,2) =Atime 

! Acquisition time. 

10150 

Array (52, 2) =Ctime 

! Coincidence time. 

10160 

Array (51, 3) =At_exp 

! Acquisition time exponent. 

10170 

Array (52, 3) =Ct_exp 

! Coincidence time exponent. 

10180 

Array (51, 4) =Paxis 

! Axis for plots. 

10190 

Array ( 52, 4) =Clip 

! Clip: 1 turn histogram clipping on; 0 turns it off. 

10200 

RETURN 



& profile scales from the Array(*). 


Date=TIMEDATE 
Mach=Array (1, 2) 

Alpha (1) =Array (1, 4) 

Time=Date 
Temp-Array (2,2) 

Alpha (2) =Array (2, 4) 
Uedge=Array (3,2) 

Alpha (3) -Array (3, 4) 

U jet_ue=Array (4,2) 
Theta=Array (4, 4) 

MAT Mod= Array (11 : 14, 1) 

MAT Tun= Array (11 : 14 , 2) 

MAT Tcsl= Array (11:14,3) 

MAT Tcs2 = Array (11:14, 4) 

MAT Index= Array (21, 1 : 3) 

MAT Thetal= Array (22 : 24, 1 : 3) 
MAT Theta3= Array (25 : 27, 1 : 3) 
MAT Beam_spc= Array (31, 1 : 3) 
MAT Focl_len= Array (32, 1 : 3) 


! Date. 

! Mach Number. 

! Angle of Attack. 

! Time. 

! Room Temperature (deg. F) . 

! Angle of Yaw. 

! Freestream Velocity. 

! Angle of Roll. 

! Jet exit velocity normalized by Uedge. 

! Tx Side Off Axis Angle. 

! Probe volume positions in MODEL coordinates. 

! Probe volume positions in TUNNEL coordinates. 

! Tx side traverse positions in Tcs8 coordinates. 

! Rx side traverse positions in Tcs8 coordinates. 

! Index of refraction of for laser light (eg: Nair,Nglass, Nwater) . 

! Angles between LASER & TUNNEL UVW laser beams in Air (N=Indexl) . 

! Angles between LASER & TUNNEL UVW laser beams in Water (N=Index3) . 

! Beam spacing at lens. 

! Focal length. 
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10420 

10430 

1044Q 

10450 

10460 

10470 

10480 

10490 

10500 

10510 

10520 

10530 

10540 

10550 

10560 

10570 

10580 

10590 

10600 

10610 

10620 

10630 

10640 

10650 

10660 Calc: 
10670 
10680 
10690 
10700 
10710 
10720 
10730 
10740 
10750 
10760 
10770 
10780 
10790 
10800 
10810 
10820 
10830 
10840 
108.50 
10860 
10870 
10880 
10890 
10900 
10910 
10920 
10930 
10940. 

10950 
10960 
10970 
10980 
10990 
11000 
11010 
11020 
11030 
11040 
11050 
11060 
11070 
11080 
11090 
11100 
11110 
11120 
11130 
11140 
11150 . 
11160 
11170 
11180 
11190 
11200 
11210 


MAT Beam_sep= Array (33, 1 : 3) 
MAT Wave_len- Array (34, 1 ; 3) 
MAT Frng_spc= Array (35, 1 : 3) 
MAT Brg_f rq= Array (36, 1 : 3) 
MAT Mix_frq= Array (37 , 1:3) 
MAT Mea_sgn= Array (38, 1:3) 
MAT Brg_sgn= Array ( 39, 1 : 3) 
MAT Mix_sgn= Array (40, 1:3) 
MAT Coin== Array (41, 1:3) 
Umin=Array (42, 1) 

Vmin=Array (42, 2) 

Wmin=Array (42,3) 

Umax=Array (43, 1) 

Vmax=Array (43 , 2) 

Wmax=Array (43, 3) 
Nreads=Array (51, 1) 
Nsam=Array (52 ,1) 

At ime=Array (51,2) 
Ctime=Array (52, 2) 

At_exp=Ar ray (51,3) 
Ct_exp=Array ( 52, 3) 
Paxis=Array (51, 4) 
Clip=Array (52, 4) 


! Beam separation angle in degrees (full angle). 

! Wave length. 

! Fringe spacing. 

! Bragg frequency . 

\ Mixing frequency. 

! Sign of measured frequency in velocity equation. 

! Sign of bragg frequency in velocity equation. 

! Sign of mixing frequency in velocity equation. 

! Coincidence criteria. 

! Frequency minimum for U calculation. 

! Frequency minimum for V calculation. 

! Frequency minimum for W calculation. 

! Frequency maximum for U calculation. 

! Frequency maximum for V calculation. 

! Frequency maximum for W calculation. 

! Number of desired samples. 

! Number of acquired samples. 

! Acquisition time. 

! Coincidence time. 

! Acquisition time exponent. 

! Coincidence time exponent. 

I Axis for plots. 

! Clip: 1 turn histogram clipping on; 0 turns it off. 


RETURN 

! This subroutine uses the current tunnel conditions and laser parameters to calculate and/or update other 
! tunnel conditions and laser parameters. 

FOR K=1 TO 3 

IF K=2 THEN 

Beaml=Theta+ATN (Beam_spc (K) /2/Focl_len (K) ) ! Angles of the off axis beam pair in air. 

Beam2=Theta-ATN ( Beam_spc ( K) /2/Focl_len (K) ) 


ELSE 


Beaml=0+ATN<Beam_spc(K) /2/Focl_len (K) ) ! Angles of the on axis beam pairs in air. 

Beam2=0-ATN(Beam_spc(K) /2/Focl_len (K) ) 


END IF 

Beaml-ASN (Index (1) /Index (3) *SIN(Beaml) ) ! Angle of the beam pairs in water. 

Beam2=ASN (Index (1) /Index ( 3) *SIN(Beam2) ) 

Beam_sep (K) =Beaml-Beam2 ! Beam pair separation angle. 

Frng_spc (K) -Wave_len (K) /Index (3) / (2*SIN (Beam_sep (K) /2) ) /1000 ! Fringe spacing in water (urn). 

NEXT K 


MAT Array (33, 1:3) =. Beam_sep ! Beam separation angle in degrees (full angle). 

MAT Array ( 35 , 1 : 3) = Frng_spc ! Fringe spacing. 

! Calculate the TCS to TUNNEL (and visa versa) traverse coordinate system transformation matrices. 


Fs=Focl_len<l) 
Fr=Focl_len(l) 
Bs=Beam_spc ( 1 ) 
Br=3 

Ts=Theta 
Tr=17 .05 


Focal length of sending side lenses (inches). 
Focal length of receiving side lenses (inches). 


Beam spacing at sending lenses (inches) . 

Receiving side lens diameter (inches) . 

Sending side off axis angle (degrees) . 

Receiving side off axis angle (degrees) . 

! Ta is the offaxis sending side auxiliary rotation angle (degrees). 

CALL Ctm_tcs (Tcs2t uni ( * ) , Tcs2tun2 ( *) , Tun2tcsl (* ) , Tun2tcs2 (*),Fs,Fr,Bs,Br, Index (*) ,Ts,Tr,Ta) 

! Calculate the LASER to TUNNEL (and visa versa) velocity coordinate system transformation matrices. 
CALL Refract (Index(*) ,Thetal (*) ,Theta3(*) ) 

CALL Ctm_ldv (Theta3 (*) , Tun21dv (*) , Ldv2tun (*) ) 

! Calculate the TUNNEL to MODEL (and visa versa) coordinate system transformation matrices. 

CALL Ctnwnod ( Alpha (*) , Mod2tun (*) , Tun2mod ( *) ) 

1 Define the coincidence mask depending on the value of Coin(*). 

Cmask=Coin ( 1 ) *l+Coin(2) *2+Coin(3) *4 
! Define Paxis$ depending on the value of Paxis. 

SELECT Paxis 
CASE 1 

Paxis$="X" 


CASE 2 

Paxis$="Y" 
CASE 3 

Paxis$="Z" 


CASE 4 

Paxis$="A" 

CASE ELSE 
Paxis=2 
Paxis$="Y" 

GOSUB M3k4 
END SELECT 

! If the Run number or File number have not been defined then have the user enter their values. 
IF Run=0 OR File=0 THEN 

CALL Enter_value ( "Run Number ", Run, "3D. 2D") 

CALL Enter_value ("File Number ", File, "3D") 

GOTO 11150 
END IF 
RETURN 
END 
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11220 

11230 

11240 

11250 

11260 

11270 

11280 

11290 

11300 

11310 

11320 

11330 

11340 

11350 

11360 

11370 

11380 

11390 

11400 

11410 

11420 

11430 

11440 

11450 

11460 

11470 

11480 

11490 

11500 

11510 

11520 

11530 

11540 

11550 

11560 

11570 

11580 

11590 

11600 

11610 

11620 

11630 

11640 

11650 

11660 

11670 

11680 

11690 

11700 

11710 

11720 

11730 

11740. 

11750 

11760 

11770 

11780 

11790 

11800 

11810 

11820 

11830 

11840 

11850 

11860 

11870 

11880 

11890 

11900 

11910 

11920 

11930 

11940 

11950 

11960 

11970 

11980 

11990 

12000 

12010 


Do_nothing: 


Menu: 

Menu read: 


Menu_disp: 


SUB Do_nothing 

! Description: 

! This subprogram is called when the keys on the keyboard are pressed during TCS8 traverse 

! movements. This is done so that any STOP, PAUSE, or RESET keys will be ignored. This prevents 

! stopping the program while the HP and TCS8 are communicating with each other. Otherwise, they 

! might get out of sync while communicating resulting in system hang ups. 

! Variables: 

! K$ String used to flush the keyboard buffer. 

K$=KBD$ 

SUBEND 

!!!!!!!! M !!! J! !!!!!!!!!!!!!! I !!!!!!! I !!! 1 !!!!!! I !!!!!!!!! J! !!!!!!!!!!!!!!!!!! I I M !!!!!!! !t !! I !!!!!!!!!!! I 
SUB Menu_read (Menu$ (*) ) 

! Description: 

! This subprogram reads in the menu descriptors for each entry of the five menus. 

! Variables: 

! Menu Used as an index to the string array Menu$(*). 

! Key Used as an index to the string array Menu$(*>. 

! Menu$(*) String array where each element describes its corresponding menu subroutine's function. 

! L$ String use to read in the menu descriptor from the data statements. 

OPTION BASE 1 
DIM L $ {80] 

! Fill all of the menu entry's descriptions with "MxKx". 

FOR Menu=l TO SIZE (Menu$, 1) 

FOR Key=l TO 8 

Menu$ (Menu, Key) ="M"&VAL$ (Menu) & "K"&VAL$ (Key) 

NEXT Key 
NEXT Menu 

ON ERROR GOTO 11570 ! The following while loop will get error#36 when the data statements run out, 

! For each menu and key, enter the menu entry's description. 

WHILE 1*1 
READ L$ 


Menu=VAL (L$ [2, 2] ) 
Key=VAL (L$ (4, 4] ) 
Menu$ (Menu, Key) =L$ 
END WHILE 
SUBEXIT 
DATA "M1K1: 


DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 
DATA 

DATA "M1K3 : 
DATA "M1K4: 
DATA "M1K5 : 
DATA "M1K6: 
DATA "M1K7 : 
"M1K8: 


Laser Alignment" 

"M2K1: Return to main menu" 


Tx & Rx" 

MODEL" 

ABSOLUTE" 


"M2K2 : Sides : 

"M2K3: Coordinates: 

"M2K4 : Mode : 

"M2K5: Move X" 

"M2K6; Move Y» 

"M2K7: Move Z" 

"M2K8 : Move A" 

"M1K2 : Pre Run" 

"M3K1: Return to MAIN menu" 

"M3K2: Enter Run & File Numbers" 

"M3K3: Enter Number of Samples". 

"M3K4 : Select Traverse Axis for Profile" 

"M3K5: Print Coordinate Transformation Matrices" 
"M3K6: Setup Graphics" 

"M3K7: Tunnel Conditions" 

"M4K1: Return to PRE RUN menu" 

"M4K2 : Load Tunnel Conditions" 

"M4K3: Save Tunnel Conditions" 

"M4K4: Print Tunnel Conditions" 

"M4K5: Enter Tunnel Condition Data" 

"M4K6: Enter Tunnel Condition Names" 

"M4K7 : Enter Tunnel Condition Units" 

"M4K8: Enter Tunnel Condition Images" 

"M3K8: Traverse" 

"M5K1: Return to PRE RUN menu" 

"M5K2: View & Set TCS8 Positions" 

"M5K3: View & Set TCS8 Units" 

"M5K4: View & Set TCS8 Revolution" 

"M5K5 : View & Set TCS8 Velocity" 

"M5K6: View & Set TCS8 Acceleration" 


Post Run (Dump Graphics) " 
Set Auto Move Positions" 
Move traverse" 

Take data" 

Auto move and take" 
Display Histograms" 


DATA 
SUBEND 

SUB Menu_disp (Menu, MenuS (*) ) 

! Description: 

! This subprogram displays the current menu at the top of the CRT. 

! Variables: 

! Menu Used as an index to the string array Menu$ (*). 
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12020 

12030 

12040 

12050 

12060 

12070 

12080 

12090 

12100 

12110 

12120 Menu_status: 
12130 
12140 
12150 
12160 
12170 
12180 
12190 
12200 
12210 
12220 
12230 
12240 
12250 

12260 Enter: 

12270 Enter_value: 

12280 

12290 

12300 

12310 

12320 

12330 

12340 

12350 

12360 

12370 

12380 

12390 

12400 

12410 

12420 

12430 

12440 

12450 

12460 Enter_st ring : 

12470 

12480 

12490 

12500 

12510 

12520 

12530 

12540. 

12550 

12560 

12570 

12580 

12590 

12600 

12610 

12620 Array: 

12630 Array_init: 
12640 
12650 
12660 
12670 
12680 
12690 
12700 
12710 
12720 
12730 
12740 
12750 
12760 
12770 
12780 
12790 
12800 


12810 


! Key Used as an index to the string array Menu?!*). 

! Menu$(*) String array where each element describes its corresponding menu subroutine's function. 

PRINTER IS CRT 

PRINT CHR$(128); ! Turn off inverse video if it is on. 

IF Menu=0 THEN Menu=l 
FOR Key =1 TO 8 

Menu$ (Menu, Key) =Menu$ (Menu, Key) &RPT? ( " ", 50-LEN(Menu? (Menu, Key) ) ) 

PRINT TABXY{1, Key) ;Menu$ (Menu, Key) (3] 

NEXT Key 
SUBEND 

SUB Menu_status (Menu, Key, Pen, Menu? (*) ) 

! Description: 

! This subprogram displays the current menu selection in normal or inverse video. The inverse 

! video text style indicates that the subroutine for the current menu selection is busy. The 

! normal text style indicates that the subroutine for the current menu selection is has completed. 

! Variables: 

! Menu Indicates which of the menus has been selected as the current menu. 

! Key Indicates which one of eight menu subroutines in the menu is to be executed. 

! Pen Indicates Busy/Ready Status. Pen=0 for busy. Pen=l for ready. 

! Menu$ (*) String array where each element describes its corresponding menu subroutine's function. 

PRINTER IS CRT 

PRINT TABXY (1, Key) ;CHR$ (129-Pen) ;Menu$ (Menu, Key) [3] ;CHR$ (128) 

WAIT .1 
SUBEND 

!!!!!!! H !!!!!!!!!!!!!!!!!!!!! H !!!!!!! I !!!! I !!!! W !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1 !!!!!!! ! 
SUB Enter_val ue (Name$ , Val ue, Image$) 

! Description: 

! This subprogram displays the current value of a variable and then has the user enter its new 

! value. The old value will be kept if the RETURN key is pressed and no data is entered. 

! Variables: 

! Name$ Name of the variable. 

! Image$ Image format of the variable. Used for printing the variable with a format. 

! Value Contains the initial value and then the updated value for the variable. 

IF Name?="Date" OR Name?="Time" THEN SUBEXIT 
DISP CHR$ (129) ; 

DISP USING 12380; Name$ 

IMAGE # , " 01 d " , K, "=" 

IF Image?<>."» THEN DISP USING "Slmage? ; Val ue 

IF Image$=" " THEN DISP USING ”#,K";Value 
DISP USING 12420;Name? 

IMAGE #," Enter new " ,K 

INPUT " ? ", Value 
DISP CHR$ (128) ; 

SUBEND 

SUB Enter_string (Name$, Value$, Image$) 

! Description: 

1 This subprogram displays the current value of a string variable and then has the user enter its 

! new value. The old value will be kept if the RETURN key is pressed and no data is entered. 

! Variables: . 

! Name$ Name of the variable. 

! Value? Contains the initial value and then the updated value for the string variable. 

DISP CHR$ (129) ; 

DISP USING 12550;Name$ 

IMAGE #, "Old ", K, "=" 

DISP USING "#, "&Image?;Value? 

DISP USING 12580; Name? 

IMAGE Enter new ",K 

INPUT " ? ", Value? 

DISP CHR? (128 ) ; 

SUBEND 

j t j i i tj i i j s i i ! i j i ij i i i j M i t i t t t t | j ! j i j i i m i m j t M m m i r T ! i i i i i i j i ( | i j i i i t j i i ; i i | [ i j ! | t m t i ! s j ; s ? j i 1 1 j ! T ! i ! 
SUB Array_init (Name?!*) , Array (M , Image? (*) , Units? (*) ) 

! Description: 

! This subprogram reads in default data for each of the variable's names, values, image formats, 

! and units. These variables include, but are not limited to, the tunnel conditions, laser 

! parameters, graph scales, traverse positions, and coordinate system transformation matrices. 

! Variables: 

! Array(*) Array of tunnel conditions, laser parameters, graph scales, etc. 

! Name?!*) Names for the variables in Array(*). 

! Image?!*) Image formats for the variables in Array!*). 

! Units? (*) Units for the variables in Array(*). 

! X Used as an index to the above arrays and string arrays. 

! Y Used as an index to the above arrays and string arrays. 

! Before Number of digits before the decimal point in the image format. 

! After Number of digits after the decimal point in the image format. 

ON ERROR GOTO 12950 
READ Y 

FOR X=1 TO SIZE (Name?, 2) 

READ Name? (Y, X) , Array (Y,X) , Image? (Y, X) , Units? (Y, X) 

SELECT Image? (Y,X) 
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12820 

12830 

12840 

12850 

12860 

12870 

12880 

12890 

12900 

12910 

12920 

12930 

12940 

12950 

12960 

12970 

12980 

12990 

13000 

13010 

13020 

13030 

13040 

13050 

13060 

13070 

13080 

13090 

13100 

13110 

13120 

13130 

13140 

13150 

13160 

13170 

13180 

13190 

13200 

13210 

13220 

13230 

13240 

13250 

13260 

13270 

13280 

13290 

13300 

13310 

13320 

13330 

13340. 

13350 

13360 

13370 

13380 

13390 

13400 

13410 

13420 

13430 

13440 

13450 

13460 

13470 

13480 

13490 

13500 

13510 

13520 

13530 

13540 

13550 

13560 

13570 

13580 

13590 

13600 

13610 


CASE "0" 

Images (Y,X) ="9D" 

CASE "1" TO "7" 

After ~VAL (ImageS ( Y, X) ) 

Before=:8-After 

Images (Y, X) =VAL$ (Before) & "D. "&VALS (After) &"D" 
CASE "K" 

CASE "N" 

CASE ELSE 

Images <Y,X) = "9D" 

END SELECT 
NEXT X 
GOTO 12780 
SUBEXIT 


! 

Y 

*** ******x==i** ******* 

******* **X=2*** **** ** 

********* x=s 3 ********* 

******* **x=4*** **** * 

DATA 

1 # 

Date , 0,0,"" , 

Mach , 0,4,"" , 

"" , 0 , 0 ,"" , 

Alphal 

, 0,4,3 

DATA 

2 , 

Time , 0,0,'“' , 

Temp , 68.5, 4, 3F , 

"" , 0 , 0 ,"" , 

Alpha2 

, 0,4,3 

DATA 

3, 

Run , 0,2,"" , 

Uedge , . 0762, 4 , m/s, 

"" , 0 , 0 ,"" , 

Alpha3 

, 0,4,3 

DATA 

4 , 

File # 0,0, "" , 

Ujet/Ue , 0,4, m/s. 

"" , . 0 , 0 ,"" # 

Theta 

, 45,4,3 

! 

Y 

******* **X=1** A ****** 

*********X=2*** ****** 

******** **x~3** ******* 

******* **x= 4 ******** 

DATA 

11 , 

Xmod , 0,4, in , 

Xtun , 0,4, in # 

Xltcs , 0,4, in , 

X2tcs 

/ 0,4, in 

DATA 

12 # 

Ymod , 0,4, in , 

Ytun , 0,4, in , 

Yltcs , 0, 4, in , 

Y2tcs 

/ 0,4, in 

DATA 

13# 

Zmod # 0,4, in , 

Ztun , 0,4, in , 

Zltcs , 0,4#in , 

Z2tcs 

, 0,4, in 

DATA 

14, 

Amod , 0,4, in # 

Atun , 0,4, in , 

Altcs , 0,4, in , 

A2tcs 

, 0,4, in 

1 

Y 

******* **x=i*** ****** 

*********X=2 ********* 

*********X=3********* 

********* x- 4 ******** 

DATA 

21 # 

Indexl , 1.000, 3,"" , 

Index2 ,1.430,3,"" # 

Index3 #1.333,3,"" # 

II II 

, 0 , 0 ,"" 

DATA 

22 , 

ThetalAU, 0,4,3 , 

ThetalAV, 90,4,3 , 

ThetalAW# 90,4,3 , 

II II 

, 0 , 0 , »» 

DATA 

23, 

ThetalBU, 45,4,3 , 

ThetalBV, 135,4,3 , 

ThetalBW, 90,4,3 , 

II II 

, 0 , 0 ,"" 

DATA 

24, 

ThetalCU, 90,4,3 , 

ThetalCV, 90,4,3 , 

ThetalCW, 0,4,3 , 

II II 

/ 0 , 0 ,"" 

DATA 

25, 

Theta3AU, 0,4,3 , 

Theta3AV, 90,4,3 # 

Theta3AW, 90,4,3 , 

II II 

1 0 , 0 ,"" 

DATA 

26, 

Theta3BU, 45,4,3 , 

Theta3BV, 135,4,3 , 

Theta3BW, 90,4,3 , 

• 1 II 

# 0 , 0 ,"" 

DATA 

27# 

Theta3CU, 90,4,3 , 

Theta3CV, 90,4,3 , 

Theta3CW, 0,4,3 , 

II II 

, 0 , 0 ,"" 

j 

Y 

******* * *x=l * ******** 

******** *x= 2 * * * ****** 

*** ******x= 3 ********* 

******* **J(ss/J******* * 

DATA 

31, 

UBeamSpc, 2. 362, 3, in , 

VBeamSpc, 2 , 362, 3, in , 

WBeamSpc, 2 . 362# 3, in , 

It It 

1 0 , 0 #"" 

DATA 

32, 

UFoclLen, 19.413, 3, in. 

VFoclLen, 19.413,3, in, 

WFoclLen, 19. 413,3, in. 

It » 

, 0 , 0 ,"" 

DATA 

33, 

UBeamSep, 0.000, 3, 3 , 

VBeamSep, 0 . 000, 3, 3 , 

WBeamSep, 0. 000, 3, 3 # 

tt H 

, 0 , 0 ,"" 

DATA 

34, 

UWaveLen, 476.5, 3, nm , 

VWaveLen, 514 . 5, 3, nm , 

WWaveLen, 488. 0, 3, nm , 

II N 

, 0 , 0 ,"" 

DATA 

35, 

UFrngSpc, 00 . 00, 3, urn , 

VFrngSpc, 00 . 00, 3, urn , 

WFrngSpc, 00 .00, 3# um # 

M II 

, 0 , 0 ,"" 

DATA 

36, 

Ubrag , 40 . 00# 4, MHz, 

Vbrag , 40 .00, 4, MHz, 

Wbrag , 40 . 00# 4, MHz, 

II M 

, 0 , 0 ,"" 

DATA 

37# 

Umix ,39. 90, 4, MHz, 

Vmix , 39.90, 4, MHz, 

WMix #39.90# 4, MHz, 

it M 

, 0 , 0 ,"" 

DATA 

38, 

UmeaSgn , +1,0,"" , 

VmeaSgn , +1,0,"" , 

WmeaSgn # +1,0,"" , 

t« II 

, 0 , 0 ,"" 

DATA 

39, 

UbrgSgn , -1,0,"" , 

VbrgSgn # -1,0,"" , 

WbrgSgn , -1#0#"" , 

It M 

, 0 , 0 ,"" 

DATA 

40, 

UmixSgn , +1,0, "" , 

VmixSgn # +1,0,"" , 

WmixSgn , +1,0# »" , 

II It 

, 0 , 0 ,"" 

DATA 

41, 

U coin # 1,0, "" , 

V coin , 1 , 0 ,"" , 

W coin # 1 , 0 #"" , 

II II 

, 0 , 0 ,"" 

DATA 

42, 

UFreqMin, -99, 4, MHz, 

VFreqMin, -99,4,MHz, 

WFreqMin, -99# 4, MHz, 

II II 

, 0 , 0 ,"" 

DATA 

43, 

UFreqMax, 99, 4, MHz, 

VFreqMax, 99, 4, MHz, 

WFreqMax, 99# 4, MHz# 

II II 

, 0 , 0 ,"" 

j 

Y 

*********X=3.********* 

*** ******x= 2 *** ****** 

******* * *x=3* ******* * 

*********X=4 ******** 

DATA 

51, 

Nreads , 1000, 0, M " , 

Atime # 30, 6, s , 

ATexp , 10,0,"” , 

Paxis 

1 2 , 0 ,"" 

DATA 

52, 

Nsam , 1000,0,"" , 

Ct ime , . IE-2 , 6 , s , 

CTexp , 5,0#"" , 

Clip 

, 0 , 0 ,"" 

j 

Y 

******* **Xsl*** **** * * 

******* **x= 2 ********* 

********* x= 3 ********* 

*** ******x=4+ ******* 

DATA 

61# 

Xminl # 0.00, l,”" , 

Xmaxl # 1 .00, 1, "" , 

Yminl , 0,0,"" , 

Ymaxl 

, 100 , 0 ,"" 

DATA 

62, 

Xmin2 , 0.00,1,"" , 

Xmax2 # 1.00,1,"" , 

Ymin2 , 0,0,"" , 

Ymax2 

, 100 , 0 , "" 

DATA 

63, 

Xmin3 , 0.00,1,"" , 

Xmax3 , 1.00,1,"" , 

Ymin3 , 0,0,"" , 

Ymax3 

, 100 , 0 ,"" 

DATA 

64, 

Xmin4 , 0.00,1,"" , 

Xmax4 , 2.00,1,"" , 

Ymin4 , 0,0,"" , 

Ymax4 

, 100 , 0 ,"" 

DATA 

65, 

Xmin5 , 0.00,1,"" , 

Xmax5 , 2.00,1,"" # 

Ymin5 , 0,0,"" , 

Ymax5 

, 100 , 0 #"" 

DATA 

66 , 

Xmin 6 , 0,1,"" , 

Xmax 6 , 3,1,"" , 

Ymin 6 , -1.5,2#"" # 

Ymax 6 

, 1.5,2#"" 

DATA 

67, 

Xmin7 , 0,1,"" # 

Xmax7 , .5,1,"" , 

Ymin7 , -1.5,2,"" , 

Ymax7 

1 1.5,2,"" 

DATA 

68 , 

Xmin 8 ,-.025,3,"" , 

Xmax 8 , .025,3,"" , 

Ymin 8 , -1.5,2,"" , 

Ymax 8 

, 1.5,2,"" 

DATA 

69, 

Xmin9 , -.1,2,"" , 

Xmax9 , .1,2,"" # 

Ymin9 , -1.5,2,"" # 

Ymax9 

, 1.5,2#"" 

! 

Y 

*********x=x********* 

*********X =2 ********* 

*********X=:3********* 

******** *x=4* ******* 

DATA 

71, 

Xminl , 935, 0, pxl , 

Xmaxl , 1235, 0, pxl. 

Yminl , 725,0, pxl, 

Ymaxl 

, 825,0, pxl 

DATA 

72, 

Xmin2 , 935,0, pxl, 

Xmax2 , 1235,0, pxl, 

Ymin2 , 585,0, pxl, 

Ymax2 

, 685,0, pxl 

DATA 

73, 

Xmin3 , 935,0, pxl, 

Xmax3 , 1235,0, pxl, 

Ymin3 , 445, 0, pxl, 

Ymax3 

, 545,0, pxl 

DATA 

74, 

Xmin4 , 935,0, pxl. 

Xmax4 , 1235,0, pxl, 

Ymin4 , 305,0, pxl, 

Ymax4 

, 405,0, pxl 

DATA 

75, 

XminS , 935,0, pxl. 

Xmax5 , 1235,0, pxl, 

Ymin5 , 165,0, pxl, 

Ymax5 

, 265,0, pxl 

DATA 

76, 

Xmin 6 , 75,0, pxl, 

Xmax 6 , 325,0, pxl, 

Ymin 6 , 525,0, pxl. 

Ymax 6 

, 825,0, pxl 

DATA 

77, 

Xmin7 , 425,0, pxl. 

Xmax7 , 675,0, pxl. 

Ymin7 , 525,0, pxl. 

Ymax7 

, 825,0, pxl 

DATA 

78, 

Xmin 8 , 75,0, pxl, 

Xmax 8 , 325,0, pxl. 

Ymin 8 , 165,0, pxl, 

Ymax 8 

, 465# 0, pxl 

DATA 

79, 

Xmin9 , 425,0, pxl, 

Xmax9 , 675,0, pxl# 

Ymin9 , 165,0, pxl, 

Ymax9 

, 465,0, pxl 

! 

Y 

******* * *X — 1 * ******** 

******** *x =2 ********* 

*********X=3*** ****** 

******** *x=4 * ******* 

DATA 

81, 

Xdivl , 5,0,"" , 

Ydivl , 4, 0, "" , 

Xdiv 6 , 6,0,"" , 

Ydiv 6 

, 6 , 0 ,»» 

DATA 

82, 

Xdi v2 , 5,0,"" , 

Ydiv2 , 4,0,"" , 

Xdiv7 , 5# 0, »" , 

Ydiv7 

, 6 , 0 ,"" 

DATA 

83, 

Xdi v3 , 5,0,"" , 

Ydiv3 , 4,0,"" , 

Xdiv 8 , 2,0,"" , 

Ydiv 8 

/ 6 , 0 , "" 

DATA 

84, 

Xdiv4 , 5,0,"" , 

Ydiv4 , 4,0,"" , 

Xdiv9 , 4,0,"" , 

Ydiv9 

, 6 , 0 , "" 

DATA 

85, 

Xdiv5 , 5,0,"" , 

Ydiv5 , 4,0,"" , 

"" , 0 , 0 ,"" , 

"" 

, 0 , 0 ,"" 


SUBEND 

Array_print: SUB Array_print (Array (*) , NameS (*) , ImageS (*) ,Units$ (*) ) 

! Description: 

! This subprogram prints the values of each of the variables with their names# image formats, arid 

! units. These variables include, but are not limited to# the tunnel conditions, laser 
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13620 

13630 

13640 

13650 

13660 

13670 

13680 

13690 

13700 

13710 

13720 

13730 

13740 

13750 

13760 

13770 

13780 

13790 

13800 

13810 

13820 

13830 

13840 

13850 

13860 

13870 

13880 

13890 

13900 

13910 

13920 

13930 

13940 

13950 

13960 

13970 

13980 

13990 

14000 Change: 
14010 Change: 
14020 
14030 
14040 
14050 • 

14060 

14070 

14080 

14090 

14100 

14110 

14120 

14130 

14140, 

14150 

14160 

14170 

14180 

14190 

14200 

14210 

14220 

14230 

14240 

14250 

14260 

14270 

14280 

14290 

.14300 

14310 

14320 

14330 

14340 

14350 

14360 

14370 

14380 

14390 

14400 

14410 


! parameters, and graph scales. 

! Variables: 

! Array (*) Array of tunnel conditions, laser parameters, graph scales, etc. 

! Name$(*) Names for the variables in Array(*). 

! ImageS (*) Image formats for the variables in Array(*). 

! Units$(*) Units for the variables in Array (*). 

! X Used as in index to the above arrays and string arrays. 

2 Y Used as in index to the above arrays and string arrays. 

PRINT USING " # , 5 / " 

FOR Y=1 TO SIZE (Array, 1} 

MAT SEARCH Array (Y, *), #LOC (<>0) ;L1 
MAT SEARCH Name$ (Y, *), #LOC (<>""); L2 
IF Ll+L2=0 AND L3=0 THEN 13980 
L3=L1+L2 

PRINT USING "I, 28X" 

FOR X-l TO SIZE (Array , 2) 

SELECT Name$ ( Y, X) 

CASE "" 

PRINT USING " # , 28X" 

CASE "Date'* 

L$=DATE$ (Array (Y,X) ) 

L$=L$ (1,2) &L$ [4,6) &L$ [8,11] 

PRINT USING 8A,A, 9A,X, 3A, 6X";TRIM$ (NameS <Y,X) ) , "=",L$, UnitsS (Y,X) 

CASE "Time'* 

L$=" "&TIMES (Array ( Y, X) ) 

PRINT USING " # , 8A, A, 9A,X, 3A, 6X»;TRIM$ (Name$(Y,X) ) L$, Unit s$ (Y, X) 

CASE ELSE 

IF ImageS (Y,X)=»" THEN ImageS (Y, X) -"9D" 

ON ERROR GOTO 13930 

PRINT USING ** # , 8A, A, " &Image$ ( Y, X) *“,X,3A, 6X"; TRIMS (NameS (Y,X) ),»=», Array (Y, X) ,Unlts$(Y,X) 
GOTO 13950 
OFF ERROR 

PRINT USING *'#, 8A, A, K,X, 3A, 6X'*;TRIM$ (NameS (Y, X) ) , Array (Y, X) ,Units$(Y,X> 

END SELECT 
NEXT X 
PRINT 
NEXT Y 
SUBEND 

!!!!!!!!!!!!!!!!!!!!!!!! I !!!!!!!!!!!!!!!! n !!!!!!!!!!!!!!( 1 !!!!!!!!!!!!!!!!!!! 2 !!!!!!! ! I!!!!!!!!!!!!!!!! H 
SUB Change (TypeS, Array (*) , NameS (*> , ImageS (*) , UnitsS (*} > 

! Description: 

! This subprogram displays on the CRT the values of each of the variables with their names, 

1 image formats, and units. The user can select one of the variables and enter a new value, 

2 name, image format, or units. The user selects the particular variable by using the 

2 left, right, up, and down cursor keys. The selected variable will appear in inverse video. 

2 When it is not selected, it will appear in normal text. When the user has selected the 

2 appropriate variable he should then press the "Select" key on the keyboard. Then, depending on 

2 the value of TypeS he will be asked to enter a new value, name, image format, or units. To 

1 exit the change variables mode press. the "Escape" key. 

2 There are three types of data that are passed to the subprogram. The first type of data 

! includes, but is not limited to, the tunnel conditions, laser parameters, and graph scales. 

I With this first type the user is allowed to enter new variable values, names, image formats, and 

! units* The second type of data is the "Auto Move and Take" data. These data are for the pre 

! programed traverse positions used in a profile scan. The third type of data is the "View and 

! Set TCS8 parameters" data acquired from and then sent back to the TCS8. 

! Variables: 

! Array (*} Array whose values, names, image formats, or units are to be modified. 

! Name$ ( *) Names for the variables in Array(*). 

! Image$(*) Image formats for the variables in Array(*). 

! UnitsS (*) Units for the variables in Array(*). 

! Type$ Indicates which type of data is to be entered. 

1 Type$="VALUES" has the user enter a new value for the selected variable. 

2 Type$="NAMES" has the user enter a new name for the selected variable. 

2 Type$=" IMAGES " has the user enter a new image format for the selected variable. 

2 Type$="UNITS" has the user enter a new units for the selected variable. 

2 X Used as in index to the above arrays and string arrays. 

2 Y Used as in index to the above arrays and string arrays. 

PRINTER IS CRT 

FOR Y=1 TO SIZE (Array, 1) 

FOR Y1=Y TO SIZE (Array, 1) 

FOR X=1 TO SIZE (Array , 2) 

IF Name$ ( Yl, X) <>" " THEN 14380 
NEXT X 
NEXT Yl 
CLEAR SCREEN 
SUBEXIT 

FOR Y2=Y1 TO SIZE (Array, 1 ) 

FOR X=1 TO SIZE (Array, 2) 

IF Name$ ( Y2, X) <>" " THEN 14430 
NEXT X 
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14 420 
14430 
14440 
14450 
14460 
14470 
14480 
14490 
14500 
14510 
14520 
14530 
14540 
14550 

14560 Wait: 

14570 

14580 

14590 

14600 

14610 

14620 Kbd: 

14630 

14640 

14650 Display: 

14660 

14670 

14680 

14690 

14700 

14710 

14720 

14730 

14740 

14750 

14760 

14770 

14780 

14790 

14800 

14810 

14820 

14830 

14840 

14850 

14860 

14870 

14880 Select: 

14890 

14900 

14 910 

14920 

14930 

1494Q 

14950 

14960 

14970 

14980 

14990 

15000 

15010 

15020 

15030 

15040 

15050 

15060 

15070 

15080 

15090 

15100 

15110 

15120 

15130 

15140 

15150 

15160 

15170 

15180 

15190 

15200 

15210 


GOTO 14440 
NEXT Y2 

FOR Y2=Y2 TO SIZE (Array, 1) 

FOR X=1 TO SIZE (Array, 2) 

IF Name$ ( Y2, X) <>" M THEN 14490 
NEXT X 
NEXT Y2 
Y2=Y2-1 
CLEAR SCREEN 

CALL Display (Type$, Y1 , Y2, Array {*) , Name$ ( * ) , ImageS (*> , Units$ (* ) ) 

Done^O 

X=1 

Y=Y1 

ON KBD ALL, 15 GOSUB Kbd 
IF NOT Done THEN Wait 
OFF KBD 
CLEAR SCREEN 
Y-Y2 
NEXT Y 
SUBEXIT 

CALL Update (TypeS , X, Y, Y1 , Y2, Done, Array ( *) , Name$ (*) , Image$ (*) , UnitsS (*) > 

RETURN 

SUBEND 

SUB Display (Type$ , Yl, Y2, Array ( * ) , Name$ ( * ) , Image$ ( *) , Units$ (*) ) 

! Description: 

! This subprogram displays on the CRT the values of each of variables with their names, image 

! formats, and units. 

! Variables: 

! Array{*> 

! Name$ (*} 

! ImageS ( *) 

! UnitsS (*) 

! TypeS 


! X 

! Y 

FOR Y=Y1 TO Y2 

FOR X = 1 TO SIZE (Array, 2) 

CALL Select (TypeS , X, Y, Yl , Y2, 0, Array (* ) , NameS (*) , ImageS (*) , UnitsS (*) ) 

NEXT X 
NEXT Y 

CALL Select (TypeS , 1, Yl , Yl , Y2, 1, Array (*) , NameS (*) , ImageS (*) , UnitsS (*) ) 

SUBEND 

SUB Select (TypeS, X, Y, Yl, Y2, C, Array < *) , NameS (*) , ImageS (*) ,Units$<*) ) 

! Description: 

! This subprogram displays on the CRT .the value of one variable along with its names, image 

! format, and units. 

Array whose values, names, image formats, or units are to be modified. 

Names for the variables in Array {*) 

Image formats for the variables in Array (*) 

Units for the variables in Array (*) 

Indicates which type of data is to be entered. 

Type$=" VALUES" has the user enter a new value for the selected variable. 

Type$= "NAMES" has the user enter a new name for the selected variable. 

Type$=" IMAGES" has the user enter a new image format for the selected variable 

Type$="UNITS" has the user enter a new units for the selected variable. 

Used as in index to the above arrays and string arrays. 

Used as in index to the above arrays and string arrays. 

PRINT CHR$ (128+C) ;TABXY (26*X-24, 15+Y-Y1+1 ) / 

PRINT RPT$ ( " ", 23) ;TABXY (26*X-24, 15+Y-Y1+1) ; 

IF Name$ ( Y, X) ="" AND Array (Y,X)=0 THEN 15260 
Img$=Image$ (Y, X) 

Unt$=Units$ ( Y, X) 

IF ImageS (Y, X) 3 ="" THEN Img$="K" 

IF UnitsS (Y, X) ="" THEN Unt$=" 

SELECT TypeS 
CASE "VALUES" 

SELECT NameS (Y,X) 

CASE "Date" 

CASE "Time" 

CASE ELSE 

PRINT USING " # , 10A, A, "& Img$& ", X, 3A" ; NameS ( Y, X) Array (Y, X) , UntS 
END SELECT 
CASE "NAMES" 

PRINT USING " #, 10A, A, 8A"; NameS (Y, X) , " : " , NameS (Y, X) 

CASE "UNITS" 


Variables : 
Array ( *) 
NameS (*) 
ImageS (*) 
UnitsS (*) 
TypeS 


Array whose values, names, image formats, or units are to be modified. 

Names for the variables in Array (*) . 

Image formats for the variables in Array(*) . 

Units for the variables in Array(*) . 

Indicates which type of data is to be entered. 

Type$="VALUES" has the user enter a new value for the selected variable. 

Type$="NAMES" has the user enter a new name for the selected variable. 

Type$=" IMAGES" has the user enter a new image format for the selected variable 
Type$="UNITS" has the user enter a new units for the selected variable. 

Used as in index to the above arrays and string arrays. 

Used as in index to the above arrays and string arrays. 
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15220 

15230 

15240 

15250 

15260 

15270 

15280 Update: 

15290 

15300 

15310 

15320 

15330 

15340 

15350 

15360 

15370 

15380 

15390 

15400 

15410 

15420 

15430 

15440 

15450 

15460 

15470 

15480 

15490 

15500 

15510 

15520 

15530 

15540 

15550 

15560 

15570 

15580 

15590 

15600 

15610 

15620 

15630 

15640 

15650 

15660 

15670 

15680 

15690 

15700 

15710 

15720 

15730 

15740. 

15750 
15760 
15770 
15780 
15790 
15800 
15810 
15820 
15830 
15840 
x 15850 
15860 
15870 
15880 
15890 
15900 
15910 
15920 
15930 
15940 
15950 
15960 
15970 
15980 
15990 
16000 
16010 


PRINT USING " 10A, A, 8A"; NameS (Y, X) , " : ", Unit s$ (Y, X) 

CASE “IMAGES” 

PRINT USING »#,10A,A,8A";Name$(Y,X),":\Image$(Y,X> 

END SELECT 
PRINT CHR$ ( 128) ; 

SUBEND 

SUB Update (TypeS, X,Y, Yl, Y2, Done, Array (*) ,Name$<*) ,Image$(*) ,Units$(*) ) 

! Description: 

! This subprogram scrolls through the variables displayed on the CRT and has the user enter 

! updated values. The user can select one of the variables and enter a new value, name, image 

! format, or units. The user selects the particular variable by using the left, right, up, down 

! cursor keys. This subprogram will only have been called after a keyboard key has been pressed. 

! If a cursor key has been pressed then the previously selected variable will be redisplayed In 

! normal text and the new selected variable will appear in inverse video text. When the user has 

! selected the appropriate variable he will have pressed the “Select" key on the keyboard. Then, 

! depending on the value of the Type$ he will be asked to enter a new value, name, image format, 

! or units. To exit the change variables mode the user will have pressed the "Escape" key. 

! Variables: 

I Array (*) Array of tunnel conditions, laser parameters, graph scales, etc. 

! Name$(*) Names for the variables in Array(*). 

! Image$ ( * ) Image formats for the variables in Array(*). 

! UnitsS (*) Units for the variables in Array (*). 

! Type$ Indicates which type of data is to be entered. 

I Type$*"VALUES" has the user enter a new value for the selected variable. 

! Type$= "NAMES" has the user enter a new name for the selected variable. 

! Type$=“ IMAGES" has the user enter a new image format for the selected variable. 

! Type$="UNITS" has the user enter a new units for the selected variable. 

I X Used as in index to the above arrays and string arrays. 

! Y Used as in index to the above arrays and string arrays. 

DISABLE 

K$=KBD$ 

IF K$-"“ THEN 15990 
SELECT NUM (K$ [1,1] ) 

CASE 27 ! ESC 

Done-1 
CASE 255 

CALL Select (Type$, X, Y, Yl, Y2, 0, Array (*), Name$ ( *} , ImageS <*) ,Units$ (*) ) 

SELECT NUM (K$ [2,2]) 

CASE 73,80 ! Break, Stop' 

PAUSE 

CASE 124 ! Menu 

Done=l 

CASE 38 ! Select 

CALL Select (Type$ , X, Y, Yl, Y2, 1 , Array ( * ) , Name? ( *) , ImageS ( *) , UnitsS ( *) ) 

SELECT Type$ 

CASE "VALUES" 

IF Name$ (Y,X)="“ THEN CALL Enter__string ("Name for "&Name$ (Y, X) , NameS (Y, X) , "K" ) 

IF Images (Y,X) THEN CALL Enter_st r ing ( “Image for " &Name$ (Y, X) , ImageS (Y, X) , "K") 

CALL Enter_value (NameS ( Y, X) , Array (Y,X), ImageS (Y,X) ) 

CASE "NAMES" 

CALL Enter_string ("Name for "&Name$ (Y,X) , NameS (Y, X) , "K") 

CASE "UNITS" 

CALL Enter^string ("Units for "& NameS (Y, X) , UnitsS ( Y, X) , "K") 

CASE "IMAGES" 

CALL Enter_st ring ( "Image for "&Name$ ( Y, X) , ImageS ( Y, X) , "K" ) 

END SELECT 

CALL Select (TypeS , X, Y, Yl , Y2, 0 , Array ( * ) , NameS (*), ImageS (*) , UnitsS ( *) ) 

IF X=SIZE (Array, 2) THEN Y=Y+1 
X=X+1 


CASE 60 

! Left 

X=X-1 


CASE 62 

! Right 

X=X+1 


CASE 94 

! Up 

Y=Y-1 


CASE 86 

! Down 

Y=Y+1 


CASE 92 

! First 

X=1 


Y=1 


END SELECT 



X=(X-1) MOD SIZE (Array, 2) +1 
Y= (Y-Yl+l-l) MOD (Y2-Y1+1 ) +Y1 
IF X<1 THEN X=SIZE (Array, 2) 

IF Y<Y1 THEN Y=Y2 

CALL Select (TypeS, X, Y, Yl, Y2, 1 , Array (* ) , NameS ( *) , ImageS (*) , UnitsS ( *) ) 
END SELECT 
ENABLE 
SUBEXIT 
SUBEND 
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16020 

16030 

16040 

16050 

16060 

16070 

16080 

16090 

16100 

16110 

16120 

16130 

16140 

16150 

16160 

16170 

16180 

16190 

16200 

16210 

16220 

16230 

16240 

16250 

16260 

16270 

16280 

16290 

16300 

16310 

16320 

16330 

16340 

16350 

16360 

16370 

16380 

16390 

16400 

16410 

16420 

16430 

16440 

16450 

16460 

16470 

16480 

16490 

16500 

16510 

16520 

16530. 

16540 

16550 

16560 

16570 

16580 

16590 

16600 

16610 

16620 

16630 

16640 

16650 

16660 

16670 

16680 

16690 

16700 

16710 

16720 

16730 

16740 

16750 

16760 

16770 

16780 

16790 

16800 


Misc: 

Convert2words: 


Errors 


Scale : 


Table: 

Table: 


!!!!!!!!!!! I !!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! I !!! 1 !!! M !!!!!!! ! 
SUB Convert2words (Real# INTEGER High, Low) 

! Description: 

! This subprogram converts a single real precision variable into two 16 bit words. The initial 

! real precision variables is converted in to a 32 bit integer and then separated into high and 

! low 16 bit integers. The most significant 16 bits will be in the "High" variable while the 

! least significant 16 bits will be placed the the "Low" variable. The main purpose of this 

! subprogram is to provide a means to send a 32 bit integer to the LVDAS over the 16 bit high 

! speed interface. 

! Variables: 

! Real Initial real precision value for the variable. 

! Hex$ Hex value of "Real". String length will be 8 bytes for 32 bits. 

! High Most significant 16 bits of integerized "Real". 

! Low Least significant 16 bits of integerized "Real". 

Hex$=DVAL$ (Real, 16) 

High=IVAL (Hex$ [1, 41,16) 

Low«I VAL ( Hex$ [5,81,16) 

SUBEND 
SUB Error 

Description: ' 

This subprogram will print an error message when ever a program error occurs. The error message 
will be displayed at the top of the CRT and also printed on the printers paper. Such errors 
might occur when data to be printed will not fit in the image formats. Other errors will also 
generate a displayed and printed error message. 

BEEP 

DISP ERRM$ 

OUTPUT PRT; ERRM$ 

Prt=VAL ( SYSTEMS ( "PRINTER IS")) 

PRINTER IS CRT 

PRINT TABXY( 95,1) ;ERRM$ 

PRINTER IS Prt 
ERROR SUBEXIT 
SUBEND 
SUB Scale (G) 

Description: 

This subprogram selects one of nine histogram or profile plots. The plot's area of the CRT is 
selected and scaled to the appropriate scales. 

OPTION BASE 1 

COM /Graph/ Wndw ( * ) , Vwprt (*> ,Xdiv<*) , Ydiv <*) , Xlabel$ < *) , YlabelS (* ) , Title$ <*) , Ximage$ < *) , Yimage$ (*) , 
Legend$ ( *) 

VIEWPORT Vwprt (G, 1) /10 . 23 , Vwprt (G, 2) / 10 .23 , Vwprt (G, 3) /10. 23, Vwprt (G,4)/10.23 
WINDOW Wndw (G, 1) , Wndw (G, 2) , Wndw (G, 3) , Wndw(G, 4) 

SUBEND 

I !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!.!!!!!!!!!!!!! 1 !!!!!!!!!!!!!! ! !!!!!!!!!! !! !!!!!!!!!!! ! ! ! ! I ! ! ! ! ! ! !!!!!!! ! I ! ! ! 
SUB Table (Table(*) ) 

Description: 

This subprogram is used to create a lookup table array. The lookup table array facilitates 
the rapid conversion of raw encoded Macrodyne data into a usable frequency. Once the table 
has been filled, then the raw Macrodyne data can be used as an index to the table array. 
Variables: 

Lookup table of frequencies. 

The 10 bit mantissa part of the raw Macrodyne data (0 .. 1023) . 

The 1 bit Fringe Count part of the raw Macrodyne data (0:16, 1:8 fringes) . 

The 4 bit Exponent part of the raw Macrodyne data. 

An array of measurement times for a given number of Fringes and Exponent. 

An array of measured frequencies for a given number of Fringes and Exponent. 

! Bin Used to index Mantissa<*). 

Used as a subrange index for Table (*). 

Used as a subrange index for Table (*). 


Table ( *) 
Mantissa (*) 
Fringes 
Exponent 
Time (* ) 

Freq (* ) 

Bin 

Min 

Max 


OPTION BASE 1 

REAL Mantissa (0:1.023) , Time (0:1023 ), Freq (0:1023) 

! If the last entry in the table in not zero then the table has already been created. 
IF Table (32766) THEN SUBEXIT 

FOR Bin=0 TO 1023 ! Fill Mantissa array. 

Mantissa (Bin) =Bin 
NEXT Bin 
Mantissa (0) =1 
Min=0 

FOR Fringes=0 TO 1 

FOR Exponent=0 TO 15 
Max=Min+1023 
IF Max=32767 THEN 
Max=32766 

REDIM Mantissa (0:1022) , Time (0; 1022) , Freq(0:1022) 

END IF 

DISP Fringes, Exponent 

!MAT Time= Mantissa* (2 A (Exponent-1 ) /500000000) 

MAT Time= Mantissa* (2 A (Exponent-3 ) /500000000) 

MAT Freq= <2 A < 4-Fringes) ) /Time 
MAT Freq= Freq/ (1000000) 


0 indicates 16 fringes while 1 indicates 8 fringes. 


Maximum size of an array is 32766. 


! Use this line with new macrodynes. 
! Use this line with old macrodynes. 
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16810 

16820 

16830 

16840 

16850 

16860 Lvdas: 

16870 Lvdas_init: 

16880 

16890 

16900 

16910 

16920 

16930 

16940 

16950 

16960 


17020 

17030 

17040 


MAT Table (Min:Max)= Freq 
Min-Min+1024 
NEXT Exponent 
NEXT Fringes 
SUBEND 

!!!!!!!!!!!! 1 !!!!!!!!!!!!! 1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! I !!!!!!!!!!!!!!!!!!!!!!!!! 1 !!!!!!!!!!! 1 1 ! 
SUB Lvdas_init (GLvdas) 

! Description: 

! This subprogram is used to initialize the HP98622-66501 Rev B 16-bit General Purpose 

! Input Output (GPIO) interface. The subprogram also opens the LVDAS path on the HP computer 

! for command and data transfer. The I/O path is given the name "GLvdas". Data transferred 

! from the HP to the LVDAS will use the "OUTPUT GLvdas" statement. Data transferred to the HP 

! from LVDAS will use the "ENTER QLvdas" statement. 

! The I/O path has a select code of 12 and is initialized to perform unformatted word 

! transfers without any end of line designations. The DIP switches on the HP98622-66501 Rev B 

1 printed circuit board need to be set as shown below: 


PCTL=1 


16970 

! DIP 

switches 

for 

INT LVL : 

Bitl-0 

Bit0=0 




16980 

! DIP 

switches 

for 

Select Code : 

Bit4=0 

Bit3=l 

Bit2=l 

Bitl=0 

Bit0=0 

16990 

! DIP 

switches 

for 

DI15to08 elk: 

RDY =1 

BSY =0 

RD =1 



17000 

! DIP 

switches 

for 

DI07t oOO elk: 

RDY =1 

BSY =0 

RD' =1 



17010 

! DIP 

switches 

for 

Hndsk Levels: 

DOUT=0 

DIN =0 

HSHK=1 

PSTS=0 

PFLG-0 


ASSIGN GLvdas TO 12; WORD, FORMAT OFF, EOL 
OUTPUT GLvdas USING »#,AA";"HP" 


SUBEND 


17050 Lvdas_take: 

17060 

17070 

17080 

17090 

17100 

17110 

17120 

17130 

17140 

17150 

17160 

17170 

17180 

17190 

17200 

17210 

17220 

17230 

17240 

17250 

17260 

17270 

17280 

17290 

17300 

17310 

17320 

17330, 

17340 

17350 

17360 

17370 

17380 

17390 Data: 


SUB Lvdas_take (GLvdas , At ime, Ct ime, INTEGER At_exp, Ct_exp, Cmask, Nsam) 

! Description: 

! This subprogram samples the two analog, three digital, and two external trigger channels 

! from the LVDAS. The HP sends a "CS" to sample the LVDAS data with coincidence. Following the 

! "CS" the HP sends the LVDAS an additional eight words to specify the acquisition and 

! coincidence times, the interarrival and coincidence time exponents, the coincidence mask, and 

! the number of desired samples. After the desired number of samples is acquired or the desired 

! acquisition time expires then the LVDAS sends to the HP an updated number of samples (Nsam) . 

! The updated Nsam may be less that the original Nsam if the desired acquisition time expires 

! before the desired Nsam samples are realized. 

! Variables: 

! Atime The maximum desired acquisition time (seconds) . 

! Ctime The maximum desired coincidence time (seconds) . 

! Atl The upper word of integer of 1000000Q*At ime. 

! At2 The lower word of integer of 10000000*At ime. 

! Ctl The upper word of integer of 10000000*Ctime. 

! Ct2 The lower word of integer of 10000000*Ct ime. 

! At_exp Exponent for interarrival times. 

I Ct_exp Exponent for coincidence times. 

! Nsam Number of desired samples. 

! Cmask Coincidence Mask for U,V, W selection. 

! Raw(*) Array of raw data acquired LVDAS data. 

OPTION BASE 1 

COM /Data/ INTEGER Raw (*) , Valid (*), REAL Table (*), Ui (*), Vi (*), Wi (*) ,Ai (*>, Bi (*) , li (*) , Ci (*> 

INTEGER At 1 , At2, Ctl, Ct2 
DISP "Taking Data" 

CALL Convert2words (At ime*10000000, Atl, At2) 

CALL Convert2words (Ctime*100000Q0, Ctl, Ct2) 

OUTPUT GLvdas USING "AA, 8 (W) " ; "CS" , At 1 , At2, Ct 1, Ct2, At_exp, Ct_exp, Cmask, Nsam 
ENTER GLvdas USING "#,W";Nsam 
IF Nsam=0 THEN SUBEXIT 
REDIM Raw(l:Nsam, 1:10) 

ENTER GLvdas USING "#,W";Raw(*) 

SUBEND 

!' !! !!! !!!!!!!!!!! !!!!!! I !!!!!! 1 !!!!!!!!!!!!!!! !!!!!! I !!!!!!!!!!!!!!!!! !!! l l 1 1 ! !!! !!!!!!!! !!! ! I' I !! 1 !!!!!!! ! 


17400 Data_reduce: 

17410 

17420 

17430 

17440 

17450 

17460 

17470 

17480 

17490 

17500 

17510 

17520 

17530 

17540 

17550 

17560 

17570 

17580 

17590 

17600 


SUB Data_reduce (INTEGER At_exp, Ct_exp, Nsam) 

! Description: 

! This subprogram separates the ten by Nsam Raw(*) data array into multiple one by Nsam 

! arrays. The frequency arrays Ui,Vi,Wi are extracted from columns 6,7,8 of the Raw data array. 

! The voltage arrays Ai,Bi are extracted from columns 9,10 of the Raw data array. The 

! interarrival time array Ii is extracted from columns 1 of the Raw data array. The coincidence 

! time array Ci is extracted from columns 2 of the Raw data array. The validation word array 

! Valid is extracted from columns 5 of the Raw data array. If i'th sample acquired contains 

! valid data, then Valid (i) will be equal to one, and zero otherwise. All values for the Valid 

! array are initially set to one by the LVDAS. 

! The raw data from arrays Ui,Vi,Wi are converted into frequencies by using their initial 

! values as indexes to the frequency look up table array Table(*). The raw data from arrays 

! Ai , Bi are converted into voltages by multiplying their initial values by 5 volts over 2 A 15. 

! The raw data from array Ii are converted into interarrival times by multiplying their initial 

I values by 2 A At_exp over 10 to get us. The raw data from array Ci are converted into 

! coincidence times by multiplying their initial values by 2 A Ct_exp over 10 to get us. 

I Variables: 


! 


Table ( *) 
Raw (*) 

Ui (*) 
Vi(*) 


Lookup table of frequencies. 

Array of raw data acquired LVDAS data. 
Array of extracted raw U frequency data. 
Array of extracted raw V frequency data. 
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17610 
17620 
17630 
17640 
17650 
17660 
17670 
17680 
17690 
17700 
17710 
17720 
17730 
17740 
17750 
17760 
17770 
17780 
17790 
17800 
17810 
17820 
17830 
17840 
17850 
17860 
17870 
17880 
17890 
17900 
17910 
17920 
17930 
17940 
17950 
17960 
17970 
17 980 

17990 Data_clip: 
18000 
18010 
18020 
18030 
18040 
18050 
1806.0 
18070 
18080 
18090 
18100 
18110 
18120 
18130 . 

18140 

18150 

18160 

18170 

18180 

18190 

18200 

18210 

18220 

18230 

18240 

18250 

18260 

18270 

18280 

18290 

18300 

18310 

18320 

18330 

18340 

18350 

18360 

18370 

18380 

18390 

18400 


! Wi (*) Array of extracted raw W frequency data. 

! Ai (*) Array of extracted raw A voltage data. 

! Bi (*) Array of extracted raw B voltage data. 

! Ii (*) Array of extracted raw interarrival time data. 

! Ci (*) Array of extracted raw coincidence time data. 

! Valid (*) Array of extracted raw validation words. 

I At_exp Exponent of interarrival times. 

! Ct_exp Exponent of coincidence times. 

! Nsam Number of samples acquired. 

OPTION BASE 1 

COM /Data/ INTEGER Raw(*) , Valid (*), REAL Table (*), Ui (*), Vi (*), Wi (*) ,Ai (*) , Bi <*), Ii (*>, Ci (*) 

REDIM Ui (Nsam) , Vi (Nsam) , Wi (Nsam) , Ai (Nsam) , Bi (Nsam) , Ii (Nsam) , Ci (Nsam) , Valid (Nsam) 

DISP "Reducing Data" 

MAT Ii* Raw(*,l) 

MAT Ci* Raw ( * , 2) 

MAT Valid* Raw(*,5) 

MAT Ui* Raw ( * , 6) 

MAT Vi* Raw ( * , 7) 

MAT Wi* Raw ( * , 8) 

MAT Ai* Raw { * , 9) 

MAT Bi* Raw ( *, 10) 

FOR K*1 TO Nsam 

Ui (K) =Table (Ui (K) ) 

Vi (K) *Table (Vi (K) ) 

Wi (K) -Table (Wi (K) ) 

NEXT K 

MAT Ai= Ai* (5/32768) 

MAT Bi* Bi* (5/32768) 

MAT Ii* Ii* (2 A At_exp/10) 

MAT Ci* Ci* (2 A Ct_exp/10) 

MAT Ui* Ui . Valid 
MAT Vi* Vi . Valid 
MAT Wi* Wi . Valid 
MAT Ai* Ai . Valid 
MAT Bi* Bi . Valid 
MAT Ii* Ii . Valid 
MAT Ci* Ci . Valid 
SUBEND 

SUB Data_clip (INTEGER Nsam, REAL Umin, Umax, Vmin, Vmax, Wmin, Wmax) 
l Description: 

! This subprogram compares each of the instantaneous U,V, and W frequencies with user 

! selectable minimum and maximum frequencies. If the instantaneous value is less than the 

! desired minimum, then the validation word is set to zero. Also, if the instantaneous value is 

! greater than the desired maximum, then the validation word is set to zero. The setting of the 

! validation words to zero will have the net effect of discarding the data samples from the data 

I set. In other words, the data is weighted as zero for the average, sdv, normal and shear 

i stress calculations. 

! Variables: 

! Nsam Number of samples acquired. 

! Ui(*) Array of instantaneous U frequencies (MHz). 

! Vi (*) Array of instantaneous V frequencies (MHz). 

! Wi (*) Array of instantaneous W frequencies (MHz) . 

! Valid (*) Array of sample validation words. 

I Umin The minimum acceptable U frequency (MHz) . 

! Umax The maximum acceptable U frequency (MHz) . 

! Vmin The minimum acceptable V frequency (MHz). 

! Vmax The maximum acceptable V frequency (MHz) . 

! Wmin The minimum acceptable W frequency (MHz) . 

! Wmax The maximum acceptable W frequency (MHz) . 

OPTION BASE 1 

COM /Data/ INTEGER Raw ( * ) , Val id ( * ) , REAL Table (*), Ui (*), Vi (*) , Wi (*) ,Ai (*), Bi (*), Ii (*), Ci (*) 

DISP "Clipping Histograms" 

FOR K=1 TO Nsam 

MAT SEARCH Ui ( * ) , LOC ( <Umi n) ; L, K 
IF L<Nsam THEN Valid (L)=0 
K=L 

NEXT K 

FOR K=1 TO Nsam 

MAT SEARCH Ui (* ) , LOC ( >Umax) ; L, K 
IF L<Nsam THEN Valid (L)*0 
K=L 

NEXT K 

FOR K=1 TO Nsam 

MAT SEARCH Vi (*) , LOC (<Vmin) ;L, K 
IF L<Nsam THEN Valid(L)=0 
K=L 

NEXT K 

FOR K=1 TO Nsam 

MAT SEARCH Vi ( *) , LOC (>Vmax) ;L, K. 

IF L<Nsam THEN Valid(L)*0 
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18410 

18420 

18430 

18440 

18450 

18460 

18470 

18480 

18490 

18500 

18510 

18520 

18530 

18540 

18550 

18560 

18570 

18580 

18590 

18600 

18610 Data_f convert 

18620 

18630 

18640 

18650 

18660 

18670 

18680 

18 690 

18700 

18710 

18720 

18730 

18740 

18750 

18760 

18770 

18780 

18790 

18800 

18810 

18820 

18830 

18840 

18850 

18860 

18870 

18880 

18890 

18900 

18910 

18 920 

18930, 

18940 

18950 

18960 

18970 

18980 

18990 

19000 Data_sum: 

19010 

19020 

19030 

19040 

19050 

19060 

19070 

19080 

19090 

19100 

19110 

19120 

19130 

19140 

19150 

19160 

19170 

19180 

19190 

19200 


K=L 

NEXT K 

FOR K=1 TO Nsam 

MAT SEARCH Wi ( *) , LOC ( <Wmin) ;L, K 

IF L<Nsam THEN Valid (L>=0 

K=L 

NEXT K 

FOR K=1 TO Nsam 

MAT SEARCH Wl ( *) , LOC ( >Wmax) ;L, K 
IF L<Nsam THEN Valid (L)=0 
K=L 

NEXT K 

MAT Ui = Ui . Valid 
MAT Vi= Vi . Valid ' 

MAT Wi= Wi . Valid 
MAT Ai = Ai . Valid 
MAT Bi= Bi . Valid 
MAT 11= Ii . Valid 
MAT Ci= Ci . Valid 
SUBEND 

SUB Data_fconvert (Array (*) ) 

! Descriptions 

! This subprogram takes the frequency values from the arrays Ui,Vi,Wi and replaces them with 

! velocities after doing the frequency to velocity conversion. 

! Variables: 

! Array (*) An array containing relevant LDV laser and tunnel condition parameters 

! Frng_spc(* Fringe Spacings extracted from Array(*). 

! Brg_frq(*) Bragg Frequencies extracted from Array(*). 

! Mix_frq(*) Mixing Freqs. extracted from Array(*) . 

! Mea_sgn(*) Measured Freq's. Signs extracted from Array (*) 

! Brg_sgn ( *) Bragg Freq's. Signs extracted from Array (*) . 

! Mix__sgn(*) Mixing Freq's. Signs extracted from Array(*). 

! Ui (*) Array of instantaneous U data. 

! Vi (*) Array of instantaneous V data. 

! Wi(*) Array of instantaneous W data. 

! Equations: .• 

! The following equations are used to convert the frequencies to velocities 

I Velocity = Fs * Ftotal 

! Ftotal = MeaSgn*Fmeas+BrgSgn*Fbrag+MixSgn*Fmix 

OPTION BASE 1 

COM /Data/ INTEGER Raw<*) , Valid (*), REAL Table <*), Ui (*), Vi <*), Wi <*) ,Ai <*), Bi (*), Ii (*) ,Ci (*) 

DIM Frng_spc (3) ,Brg_frq(3) ,Mix_frq(3} ,Mea_sgn (3) , Brg_sgn(3) ,Mix_sgn(3) 

DISP "Converting Data" 

MAT Frng_spc= Array (35, 1 : 3) 

MAT Brg_frq= Array ( 36, 1 : 3 ) 

MAT Mix_f rq= Array ( 37 , 1 : 3 ) 

MAT Mea_sgn= Array ( 38 , 1 : 3 ) 

MAT Brg_sgn= Array (39, 1 : 3 ) 

MAT Mix_sgn= Array ( 40 , 1 : 3 ) 

MAT Ui= Ui* (Mea_sgn(l) ) 

MAT Vi = Vi* (Mea_sgn(2) ) 

MAT Wi= Wi * (Mea_sgn (3 ) ) 

MAT Ui= Ui+ (Brg_sgn (1 ) *Brg_frq (1) +Mix_sgn (1) *Mix_frq ( 1) ) 

MAT Vi= Vi+(Brg_sgn<2) *Brg_frq(2) +Mix_sgn<2) *Mix_frq<2) ) 

MAT Wi= Wi + (Brg_sgn (3) *Brg_frq (3) +Mix_jsgn (3) *Mix__frq ( 3) ) 

MAT Ui* Ui* (Frng__spc ( 1) ) 

MAT Vi= Vi* (Frng_spc (2) ) 

MAT Wi= Wi* (Frng_spc (3) ) 

SUBEND 

SUB Data_sum(Sum(*) , INTEGER N(*),Nsam) 

! Description: 

! This subprogram performs the summations on the instantaneous LDV and analog data. Data 

! will be weighted as zero in the summations if the value of the validation word is set to zero. 

! Intermediate arrays will be made so that summations of the products of the LDV and analog data 

I can be determined. 

! Variables: 

! Nsam Number of samples acquired. 

! Valid (*) Array of sample validation words. 

! Ui (*) Array of instantaneous U frequency or velocity samples, 

! Vi(*) Array of instantaneous V frequency or velocity samples. 

! Wi(*) Array of instantaneous W frequency or velocity samples. 

! Ai ( *) Array of instantaneous A voltage samples. 

! Bi ( *) Array of instantaneous B voltage samples. 

! Ii(*) Array of interarrival times. 

! Ci ( * > Array of coincidence times. 

[ Uu(*) Instantaneous product of the instantaneous Ui & Ui. 

! Vv ( * ) Instantaneous product of the instantaneous Vi & Vi. 

! Ww ( *) Instantaneous product of the instantaneous Wi & Wi. 

! Aa(*) Instantaneous product of the instantaneous Ai & Ai. 

! Bb(*) Instantaneous product of the instantaneous Bi & Bi. 
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19210 

19220 

19230 

19240 

19250 

19260 

19270 

19280 

19290 

19300 

19310 

19320 

19330 

19340 

19350 

19360 

19370 

19380 

19390 

19400 

19410 

19420 

19430 

19440 

19450 

19460 

19470 

19480 

19490 

19500 

19510 

19520 

19530 

19540 

19550 

19560 

19570 

19580 

19590 

19600 

19610 

19620 

19630 

19640 

19650 

19660 

19670 

19680 

19690 

19700 

19710 

19720 

19730. 

19740 

19750 

19760 

19770 

19780 

19790 

19800 

19810 

19820 

19830 

19840 

19850 

19860 

19870 

19880 

19890 

19900 

19910 

19920 

19930 

19940 

19950 

19960 

19970 Data_calc: 

19980 

19990 


! I2(*) Instantaneous product of the instantaneous Ii & Ii. 

! C2(*) Instantaneous product of the instantaneous Ci & Ci. 

! Uv(*) Instantaneous product of the instantaneous Ui & Vi. 

! Vw (*) Instantaneous product of the instantaneous Vi & Wi. 

! Wu(*) Instantaneous product of the instantaneous Wi & Ui. 

! Ab(*) Instantaneous product of the instantaneous Ai & Bi. 

! Ua {*) Instantaneous product of the instantaneous Ui & Ai. 

! Va(*) Instantaneous product of the instantaneous Vi & Ai. 

! Wa(*) Instantaneous product of the instantaneous Wi & Ai. 

! Sum(l,l) Summation of the array Ui. 

! Sum(2,l) Summation of the array Vi. 

! Sum(3,l) Summation of the array Wi. 

! Sum (4, 1) Summation of the array Ai. 

! Sum(5,l) Summation of the array Bi. 

! Sum(6,l) Summation of the array Ii. 

! Sum(7,l) Summation of the array Ci. 

! Sum(l,2) Summation of the array Uu. 

! Sum(2,2) Summation of the array Vv. 

! Sum(3,2) Summation of the array Ww. 

! Sum(4,2) Summation of the array Aa. 

! Sum (5,2) Summation of the array Bb. 

! Sum (6,2) Summation of the array 12. 

! Sum (7,2) Summation of the array C2. 

! Sum(l,3) Summation of the array Uv. 

! Sum(2,3) Summation of the array Vw. 

! Sum(3,3) Summation of the array Wu. 

! Sum(4, 3) Summation of the array Ab. 

! Sum(5,3) Summation of the array Ua. 

! Sum(6,3) Summation of the array Va. 

! Sum(7,3) Summation of the array Wa. 

OPTION BASE 1 

COM /Data/ INTEGER Raw ( *> , Val id < * > , REAL Table (*) , Ui <*) , Vi (*), Wi (*), Ai <*), Bi (*), Ii (*), CI (* ) 

REAL Uu (1000) , Vv( 1000) ,Ww< 1000) ,Aa<1000) , Bb(1000) , 12 <1000) , C2 (1000) 

REAL Uv (1000) , Vw( 1000) ,Wu< 1000) ,Ab (1000) , Ua (1000) ,Va<1000) , Wa (1000) 

REDIM Uu (Nsam) , Vv (Nsam) , Ww (Nsam) , Aa (Nsam) , Bb (Nsam) , 12 (Nsam) , C2 (Nsam) 

REDIM Uv (Nsam) , Vw (Nsam) , Wu (Nsam) , Ab (Nsam) , Ua (Nsam) , Va (Nsam ) , Wa (Nsam) 

DISP "Summing Data" 

I 

MAT Uu= Ui . Ui 
MAT Vv= Vi . Vi 
MAT Ww= Wi . Wi 
MAT Aa= Ai . Ai 
MAT Bb= Bi . Bi 
MAT Uv= Ui . Vi 
MAT Vw= Vi . Wi 
MAT Wu= Wi . Ui 
MAT Ab= Ai . Bi 
MAT Ua= Ui . Ai 
MAT Va= Vi . Ai 
MAT Wa- Wi . Ai 
MAT 12= Ii . Ii 
MAT C2= Ci . Ci 

t 

Sum (1,1) =SUM ( Ui ) 

Sum (2 , 1 ) =SUM ( Vi ) 

Sum(3,l)=SUM(Wi) 

Sum (4,1) =SUM ( Ai ) 

Sum ( 5 , 1 ) =SUM ( Bi ) 

Sum ( 6, 1) =SUM ( Ii) 

Sum (7 , 1) =SUM (Ci ) 

Sum(l, 2) =SUM(Uu) 

Sum (2, 2) =SUM(Vv) 

Sum (3, 2) =SUM(Ww) 

Sum (4 , 2) =SUM (Aa) 

Sum (5,2) =SUM ( Bb) 

Sum(6, 2) =SUM( 12) 

Sum (7 , 2) =SUM(C2) 

Sum(l, 3) =SUM(Uv) 

Sum(2, 3) =SUM(Vw) 

Sum(3, 3) =SUM(Wu) 

Sum(4, 3) =SUM(Ab) 

Sum (5,3) =SUM(Ua) 

Sum (6, 3) =SUM(Va) 

Sum (7,3) =SUM ( Wa ) 

MAT N= (SUM (Valid) ) 

SUBEND 

SUB Data_calc (INTEGER N ( * ) , REAL Sum<*) # U, V, W, A, B, I,C, Ul,Vl,Wl,Al f BI, II, Cl , Ulvl, Vlwl, Wlul, Albl, 
Ulal,Vlal,Wlal) 

! Description: 

! This subprogram uses the summations on the instantaneous LDV and analog data as well as the 
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20000 

20010 

20020 

20030 

20040 

20050 

20060 

20070 

20080 

20090 

20100 

20110 

20120 

20130 

20140 

20150 

20160 

20170 

20180 

20190 

20200 

20210 

20220 

20230 

20240 

20250 

20260 

20270 

20280 

20290 

20300 

20310 

20320 

20330 

20340 

20350 

20360 

20370 

20380 

20390 

20400 

20410 

20420 

20430 

20440 

20450 

20460 

20470 

20480 

20490 

20500 

20510 

20520 

20530 

20540 

20550 

20560 

20570 

20580 

20590 

20600 

20610 

20620 

20630 

20640 

20650 

20660 

20670 

20680 

20690 

20700 

20710 

20720 

20730 

20740 

20750 

20760 

20770 

20780 

20790 


! summations of the products of the LDV and analog data. The subprogram takes these summations 

! and calculates the averages, standard deviations, and shear stresses. 

! Variables: 

! Nsam The number of valid samples. 

! Sum(l,l) Summation of the array Ui. 

! Sum(2,l) Summation of the array Vi. 

! Sum(3,l) Summation of the array Wi. 

! Sum(4,l) Summation of the array Ai. 

! Sum(5, 1) Summation of the array Bi. 

! Sum(6,l) Summation of the array Ii. 

! Sum(7,l) Summation of the array Ci, 

! Sum(l,2) Summation of the array Uu. 

! Sum(2,2) Summation of the array Vv. 

! Sum(3,2) Summation of the array Ww. 

! Sum(4,2) Summation of the array Aa. 

! Sum(5,2) Summation of the array Bb. 

! Sum (6, 2) Summation of the array 12. 

! Sum(7,2) Summation of the array C2. 

! Sum(l,3) Summation of the array Uv. 

! Sum(2,3) Summation of the array Vw. 

! Sum(3,3) Summation of the array Wu. 

! Sum(4,3) Summation of the array Ab. 

! Sum(5,3) Summation of the array Ua. 

! Sum (6, 3) Summation of the array Va. 

! Sum(7,3) Summation of the array Wa. 

! U Average U frequency or velocity. 

! V Average V frequency or velocity. 

! W Average W frequency or velocity. 

! A Average A voltage. 

! B Average B voltage. 

! I Average interarrival time. 

! C Average coincidence time. 

! Ul Standard deviation for U frequency or velocity. 

! VI Standard deviation for V frequency or velocity. 

I Wl Standard deviation for W frequency or velocity. 

I Al Standard deviation for A voltage. 

I Bl Standard deviation for B voltage. 

! II Standard deviation for interarrival time. 

■ Cl Standard deviation for coincidence time. 

! Ulvl Velocity : Velocity Shear Stress. 

! Vlwl Velocity : Velocity Shear Stress. 

! Wlul Velocity : Velocity Shear Stress. 

! Albl Voltage :Voltage Shear Stress. 

! Ulal Velocity : Voltage Shear Stress. 

! Vial Velocity : Voltage Shear Stress. 

! Wlal Velocity : Voltage Shear Stress. 

DISP "Calculating Results" 

Nsam=N(l, 1) 

IF Nsam>0 THEN 

U=Sum(l, 1) /Nsam 
V=Sum (2, 1) /Nsam 
W=Sum(3,l) /Nsam 
A=Sum ( 4 , 1 ) /Nsam 
B=Sum (5, 1 ) /Nsam 
I=Sum (6, 1) /Nsam 
C=Sum(7,l) /Nsam 
U1=SQR (ABS ( Sum (1 , 2)/Nsam-U*U) ) 

V1=SQR (ABS (Sum (2,2) /Nsam-V*V) ) 

Wl=SQR(ABS<Sum(3,2)/Nsam-W*W) ) 

A1=SQR (ABS (Sum(4, 2) /Nsam-A*A) ) 

B1=SQR (ABS (Sum (5, 2) /Nsam-B*B) ) 

I1=SQR (ABS (Sum (6, 2) /Nsam-I*I) ) 

C1=SQR (ABS (Sum<7, 2) /Nsam-C*C) ) 

Ulvl=Sum (1, 3) /Nsam-U*V 
Vlwl=Sum (2, 3) /Nsam-V*W 
Wlul=Sum( 3, 3) /Nsam-W*U 
Albl=Sum (4,3) /Nsam-A*B 
Ulal=Sum (5,3) /Nsam-U*A 
Vlal=Sum (6,3) /Nsam-V*A 
Wlal=Sum < 7, 3) /Nsam-W*A 

ELSE 

U=0 

V=0 

w=o 

A~0 
B=0 
1=0 
C=0 
U1=0 
VI =0 
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20800 

20810 

20820 

20830 

20840 

20850 

20860 

20870 

20880 

20890 

20900 

20910 

20920 

20930 

20940 

20950 

20960 

20970 

20980 

20990 

21000 

21010 

21020 

21030 

21040 

21050 

21060 

21070 

21080 

21090 

21100 

21110 

21120 

21130 

21140 

21150 

21160 

21170 

21180 

21190 

21200 

21210 

21220 

21230 

21240 

21250 

21260 

21270 

21280 

21290 

21300 

21310 

21320. 

21330 

21340 

21350 

21360 

21370 

21380 

21390 

21400 

21410 

21420 

21430 

21440 

21450 

21460 

21470 

21480 

21490 

21500 

21510 

21520 

21530 

21540 

21550 

21560 

21570 

21580 

21590 


W1=0 
A1*0 
B1=0 
11=0 
C1=0 
Ulvl=0 
Vlwl=0 
Wlul=0 . 

Albl=0 
Ulal=0 
Vlal=0 
Wlal=0 
END IF 
SUBEND 

Data_trnsfrm: SUB Data_trnsfrm(REAL K<*) , U, V, W, Ul, VI, Wl, Ulvl, Vlwl, Wlul, Ulal, Vial, Wlal) 

! Description: 

! This subprogram performs a coordinate system transformation on the averages, standard 

! deviations, and shear stresses. The coordinate system transformation to be applied is passed 

! through the "K3X3" array. If a LASER to TUNNEL coordinate system transformation is to be 

! performed, then the array ,, Ldv2tun" array will be passed -to the "K3X3" array. If a TUNNEL to 

! MODEL coordinate system transformation is to be performed, then the array ,, Tun2mod" array will 


be passed to the ”K3X3" array. 


J 

Variables: 




I 

U 

Average U velocity. 



I 

V 

Average V velocity. 



J 

W 

Average W velocity. 



1 

U1 

Standard deviation for U 

velocity . 


t 

VI 

Standard deviation for V 

velocity. 


I 

Wl 

Standard deviation for W 

velocity. 


t 

Ulul 

Velocity : Velocity Normal 

Stress. 


J 

Ulvl 

Velocity : Velocity Shear 

Stress. 


J 

Ulwl 

Velocity : Velocity Shear 

Stress. 


t 

VI ul 

Velocity : Velocity Shear 

Stress. 


1 

Vlvl 

Velocity : Velocity Normal 

Stress. 


I 

Vlwl 

Velocity : Velocity Shear 

Stress . 


J 

Wlul 

Velocity : Velocity Shear 

Stress. 


1 

Wlvl 

Velocity ;Velocity Shear 

Stress . 


J 

Wlwl 

Velocity: Velocity Normal 

Stress . 


1 

Ulal 

Velocity: Voltage Shear 

Stress. 


J 

Vial 

Velocity : Voltage Shear 

Stress . 


1 

Wlal 

Velocity: Voltage Shear 

Stress . 


1 

R ( *) 

Original U,V,W. 



J 

F ( *) 

Original Ulal, Vial, Wlal . 



1 

E>(*) 

Original stress terms Ulul, Ulvl, Wlwl. 


J 

K3X3 

Coordinate system transformation matrix for average and Velocity : Voltage 

shear stress 

1 


conversions. 



I 

K9X9 

Coordinate system transformation matrix for Velocity : Velocity normal and 

shear stress 

{ 


conversions . 



1 

S(*) 

Transformed U,V,W. 



J 

H ( *) 

Transformed Ulal, Vial , Wlal . 


1 

Q ( *) 

Transformed stress terms 

Ulul, Ulvl, . . .,Wlwl. 


OPTION BASE 1 





REAL R(3) ,$(3),F(3),H(3),P(9) , Q{9) , K3x3 (3, 3) , K9x9 (9, 9) 

DISP "Transforming Results" 

! Calculate Ulul, Vlvl , Wlwl using U1,V1,W1. 

Ulul-Ul*Ul 

Vlvl=Vl*Vl 

Wlwl=Wl*Wl 

! Set Ulwl,Vlul,Wlvl equal to Wlul, Ulvl, Vlwl . 

Ulwl=Wlul 

Vlul=Ulvl 

Wlvl=Vlwl 

! Fill the matrix R with U,V,W. 

R(1)=U 
R (2) =V 
R (3) =W 

! Fill the matrix F with Ulal , Vial, Wlal . 

F (1) =Ulal 
F (2) =Vlal 
F ( 3) =Wlal 

! Fill the matrix P with Ulul, Ulvl, Ulwl,Vlul, Vlvl, Vlwl, Wlul, Wlvl, Wlwl . 
P (1) =Ulul 
P (2) =Ulvl 
P (3) =Ulwl 
P (4) =Vlul 
P < 5 > =Vlvl 
P (6) =Vlwl 
P (7) =Wlul 
P (8) =Wlvl 
P (9) =Wlwl 
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21600 

21610 

21620 

21630 

21640 

21650 

.21660 

21670 

21680 

21690 

21700 

21710 

21720 

21730 

21740 

21750 

21760 

21770 

21780 

21790 

21800 

21810 

21820 

21830 

21840 

21850 

21860 

21870 

21880 

21890 

21900 

21910 

21920 

21930 

21940 

21950 

21960 

21970 

21980 

21990 

22000 Print: 

22010 Data_printj 

22020 

22030 

22040 

22050 

22060 

22070 

22080 

22090 

22100 

22110. 

22120 

22130 

22140 

22150 

22160 

22170 

22180 

22190 

22200 

22210 

22220 

22230 

22240 

22250 

22260 

22270 

22280 

22290 

22300 

22310 

22320 

22330 

22340 

22350 

22360 

22370 

22380 


! Define the matrix K9x9 using products of the elements from then matrix K3x3 . 

FOR X=1 TO 9 

FOR Y =1 TO 9 

Yl = ( ( Y — 1 ) DIV 3) +1 
X1«((X-1) DIV 3) +1 
Y2 = ( ( Y-l) MOD 3) +1 
X2=((X-1) MOD 3) +1 
K9x9 (Y, X) =K3x3 (Yl, XI) *K3x3 (Y2, X2) 

NEXT Y 
NEXT X 

! Transform matrix R to S using K3x3. 

MAT S= K3x3*R 

! Transform matrix F to H using K3x3, 

MAT H= K3x3*F 

! Transform matrix P to Q using K9x9. 

MAT Q= K9x9*P 

! Extract the transformed U,V,W from the matrix S. 

U=S(1) 

V-S (2) 

W=S (3) 

! Extract the transformed Ulal, Vial, Wlal from the matrix H. 

Ulal=H(l) 

Vlal=H (2) 

Wlal=H {3} 

! Extract the transformed Ulul, Ulvl, Ulwl, VI ul, Vlvl,Vlwl, Wlul, Wlvl , Wlwl from the matrix Q. 

Ulul=Q (1) 

Ulvl=Q (2) 

Ulwl=Q<3> 

Vlul=Q ( 4 ) 

Vlvl=Q{5) 

Vlwl=Q<6) 

Wlul=Q ( 7) 

Wlvl=Q(8) 

Wlwl=Q ( 9) 

! Calculate U1,V1,W1 using Ulul , Vlvl, Wlwl . 

U1=SQR (ABS (Ulul) ) 

V1=SQR(ABS (Vlvl) ) 

W1=SQR <ABS (Wlwl) ) 

! Return transformed U, V, W, Ul, VI, W1 , Ulvl, Vlwl , Wlul, Ulal, Vial, Wlal to main program. 

SUBEND 

!!!!!! I !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1 !!!!!!!!! 2 I!!!!!!!!!!!!!!!!!!!!!!!!!!! 2 !!!!!!!!!> M I ! 1 !!!!!! ! 
SUB Data_print (Axis, Pos, INTEGER Nsam, C$, REAL U,V,W,A, B, I , C, Ul , VI, W1 , Al, B1 , II, Cl, Ulvl, Vlwl , Wlul , Albl, 

Ulal, Vial, Wlal, Uedge) 

! Description: 

! This subprogram prints the averages, standard deviations, and shear & normal stress in 

! tabular form. This subprogram may be called several times. The first call might print the 

! reduced velocity data when their units are in frequency (MHz) . Subsequent calls will print the 

! reduced data when their units are in velocity (m/s) . These subsequent calls will print the 

! the data in one of three coordinate systems: LASER, TUNNEL, and MODEL. 

! Variables: 

! U Average U velocity. 

! V Average V velocity. 

! W Average W velocity. 

! A Average W velocity. 

! B Average B voltage. 

! I Average interarrival time. 

! C Average coincidence time. 

! Ul Standard deviation for U velocity. 

! VI Standard deviation for V velocity. 

! W1 Standard deviation for W velocity. 

! A1 Standard deviation for A voltage. 

! Bl Standard deviation for B voltage. 

! II Standard deviation for interarrival time. 

! Cl Standard deviation for coincidence time. 

! Ulvl Velocity : Velocity Shear Stress. 

! Vlwl Velocity :Velocity Shear Stress. 

! Wlul Velocity : Velocity Shear Stress. 

! Albl Voltage :Voltage Shear Stress. 

! Ulal Velocity:Voltage Shear Stress. 

! Vial Velocity ; Voltage Shear Stress. 

! Wlal Velocity :Voltage Shear Stress. 

! Axis$ Indicates one of the three axes X,Y,Z being traversed. 

! Pos Current Traverse Position. 

! Nsam Number of samples acquired. 

! C$ Indicates units and/or coordinate system of data printed. 

DISP "Printing Results" 

ON ERROR CALL Error 
PRINTER IS PRT; WIDTH 144 
Axis$=CHR$ (NUM ( "X") +Axis-1 ) 

PRINT USING 22490;L$, Pos, U,U1, Ulvl 
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22390 

PRINT 

USING 22530; A, Al, Albl, Ulal 


22400 

PRINT 

USING 22500 ; "N" , Nsam, V, VI, Vlwl 


22410 

PRINT 

USING 22540; B, Bl, 11 

, Vial 


22420 

PRINT 

USING 22510 ; C$ [ 1 , 3 ) 

, W, Wl , Wl ul 


22430 

PRINT 

USING 22550;C, I , Cl, 

Wlal 


22440 

PRINT 

USING 22520; Uedge, Ul/Uedge, Vl/Uedge 


22450 

PRINT 

USING 22560 ; Wl/ Uedge, (Ulvl) /Uedge A 2 

(Vlwl) /Uedge A 2, (Wlul ) /Uedge A 2 

22460 

IF C$O"M0D" THEN PRINT 



22470 

PRINTER IS CRT 



22480 

OFF ERROR 



22490 

IMAGE 

#, 8X, A, "=" , 3D . 4D, " 

U=" , 3D. 5D, " 

Ul=" , 3D . 5D, " Ulvl=",4D.6D 

22500 

IMAGE 

#,8X, A, , 8D, " 

V=", 3D. 5D, " 

Vl=" , 3D . 5D, " Vlwl=",4D.6D 

22510 

IMAGE 

#,8X,3A, 7X, " 

W=", 3D. 5D, " 

Wl=" , 3D . 5D, " Wlul=", 4D. 6D 

22520 

IMAGE 

#,18X, 

UE=", 3D. 4D, " U1/UE=" , 3D. 5D, " V1/UE=» , 3D. 4D 

22530 

IMAGE 

" A =" , 3D. 5D, " 

Al =" , 3D. 5D, " 

W/UE =", 3D . 5D, " Ulal=",2D.6D 

22540 

IMAGE 

" U/UE=" , 3D. 5D, " 

V/UE =", 3D.5D, " 

IAT1=", 9D, " Vlal=", 2D. 6D 

22550 

IMAGE 

» CT=",9D, " 

IAT=", 9D, " CT1 

=", 9D, " Wlal=", 2D . 6D 

22560 

IMAGE 

" W1/UE=", ID. 4D, " Ulvl/UE2=", 2D. 4D, " Vlwl/UE2= H , 2D . 4D, " Wlul/UE2=", 

22570 

SUBEND 





2D.4D 


22580 Datajplot: 
22590 
22600 
22610 
22620 
22630 
22640 
22650 
22660 
22670 
22680 
22 690 
22700 
22710 
22720 
22730 
22740 
22750 
22760 
22770 
22780 
22790 
22800 
22810 
22820 
22830 
22840 
22850 
22860 
22870 
22880 
22890 
22900 
22 910. 

22920 

22930 

22940 

22950 

22960 

22970 

22980 

22990 

23000 

23010 

23020 

23030 

23040 

23050 

23060 

23070 

23080 

23090 

23100 

23110 

23120 

23130 

23140 Tcs8 : 

23150 Tcs8init : 

23160 

23170 

23180 


SUB Data_plot (Array (*) , Symbols ( *) , G, Y, PI, P2,P3, Scale, INTEGER N1,N2,N3) 

! Description: 

! This subprogram plots the averages, standard deviations, or shear & normal stress in 

! the 4 profile plots on the CRT. This subprogram will typically be called 4 times. The first 

! call will plot the average velocities normalized by Uedge. The second call will plot the 

! normalized standard deviations of the velocities. The third call will plot the normalized 

! velocity shear stresses. The forth and last call will plot the average and standard deviations 

! of the analog data in the forth and last plot. Data points outside the plot boundaries will 

! be plotted at the plot boundary. 

! Variables: 

! Array ( *) 

! Symbols (*) 

! G 

! Y 

I PI 

! 92 

! P3 

! Scale 

! N1 

! N2 

! N3 

! Wndw(*) 

! Vwprt ( * ) 

! Symbol (*) 

OPTION BASE 1 
DIM Wndw ( 4 ) , Vwprt (4) , Symbol (20, 3) 

DISP “Plotting Results" 

AREA PEN -1 
PEN 1 

MAT Wndw= Array ( 60+G, * ) 

MAT Vwprt = Array ( 70+G, * ) 

VIEWPORT Vwprt (1) /10. 23, Vwprt (2) /10 . 23, Vwprt ( 3) /I 0.23, Vwprt (4) /10 .23 
WINDOW Wndw (1 ) ,Wndw(2) , Wndw (3 ) , Wndw (4 ) 

CLIP ON 
FOR 1=0 TO 2 

IF 1=0 AND N1=0 THEN 23120 
IF 1=1 AND N2=0 THEN 23120 
IF 1=2 AND N3=0 THEN 23120 
Sy=I+l 

Noc=Symbols (Sy, 0, 1) 

REDIM Symbol (Noc, 3) 

MAT Symbol= Symbols (Sy, 1 : Noc, *) 

SELECT I 
CASE 0 

X=Pl*Scale 
CASE 1 

X=P2*Scale 
CASE 2 

X=P3*Scale 
END SELECT 

Xm=MIN (MAX (X, Wndw (1) ) , Wndw (2) ) 

Ym=MIN (MAX (Y, Wndw (3) ) ,Wndw(4) ) 

MOVE Xm, Ym 

SYMBOL Symbol (*), FILL, EDGE 
NEXT I 
SUBEND 

t i i t i !!!!!!!!!!!!!!!!!!!!!( 1 ! I !!!!!!!!!!!!!!!!!!! I !! I !!!!!!!!!!!!!!!!!!!!!!!!!! J !!!!!!!!!!!!! I !!!!!!!!! ! 
SUB Tcs8init (0Tcs8) 

! Description: 

! This subprogram is used to initialize this computer's internal RS232 serial interface. 

! The subprogram also opens the TCS8 path on the Hewlett Packard series 9000 model 3 XX computer 


Array containing the plot positions and scales. 

Array of Symbol arrays. Each symbol array contains a distinct geometric symbol. 
Indicates which plot that the normalized P1,P2,P3 will be plotted against Y in. 
Vertical position of the normalized data points in the plot. 

Horizontal position of the 1st data point (PI will be normalized by Scale) . 

Horizontal position of the 2nd data point (P2 will be normalized by Scale) . 

Horizontal position of the 3rd data point ( P 3 will be normalized by Scale) . 

The value that the horizontal positions will be normalized by. 

The number of samples contributing to Pi's value. PI will be plotted if N1>0. 

The number of samples contributing to P2's value. P2 will be plotted if N2>0. 

The number of samples contributing to P3's value. P3 will be plotted if N3>0. 

Array containing the plot's scales. 

Array containing the plot's CRT position. 

Array containing a distinct geometric symbol. 
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1 for command and data transfer. The I/O path Is given the name "@Tcs8 M . Data transferred 

! from the HP to the TCS8 will use the "OUTPUT @Tcs8" statement. Data transferred to the HP 

! from TCS8 will use the "ENTER @Tcs8" statement. 

! The I/O path has a select code of 9 and is initialized to perform unformatted byte 

! transfers without any end of line designations. 

ASSIGN 0Tcs8 TO 9; BYTE, FORMAT OFF, EOL "" 

CONTROL 9,0;1 ! Reset interface, 

CONTROL 9, 3; 9600 ! Select a baud rate of 9600. 

CONTROL 9,4;31 ! Select even parity, enable parity, 2 stop bits, 8 bits per character. 

CONTROL 9, 12; IVAL ("EF", 16) ! Enable Carrier Detect. Disable Data Set Ready. Disable Clear To Send. 

CONTROL 9, 13; 9600 ! Default baud rate of 9600. 

CONTROL 9, 14; 31 ! Default character format: Even parity enabled, 2 stop, 8 bits/ char. 

SUBEND 

SUB TcsSset (Commands, @Tcs8) 

! Description: 

! This subprogram allows the user to view and then set the various initialization parameters 

! of each channel of the TCS8. These parameters are the current position, counts per inch, 

! counts per revolution, motor velocity, motor acceleration, plus and minus limit switches, 

! home switch, and motor stall indication. All of these parameters can be viewed and set except 

! the limit and home switches and the stall indication. They can only be viewed. 

! Variables: 

! Commands A TCS8 command string which indicates which parameter we want to view 4 set. 

! View(*) Array of old TCS8 parameters viewed (received from TCS8) . One for each channel. 

! Set ( *) Array of new TCS8 parameters to be set (sent to TCS8) . One for each channel. 

! NameS (*) String array of TCS8 parameter names. 

! ImageS (*) String array of image formats. 

! UnitsS (*) String array of units. 

! Channel Indicates the TCS8 channel number. Used to index the above arrays. 

OPTION BASE 1 

DIM V i ew ( 8 , 1 ) , Set (8,2) , Name$ (8,1) [10] , Image$ (8, 1) [10} ,Units$ (8, 1) [10] 

OUTPUT @Tcs8 USING “K, /"; "V"4Command$4"0" ! Tell the TCS8 we want to View a parameter. 

ENTER @Tcs8 USING "8 (K) M ;View (*) ! Enter the parameter specified by Commands. 

! Initialize the NameS, ImageS, UnitsS and Set arrays. 

READ Name$ <*) 

MAT Image$= ("6D.3D") 

DATA X1,X2,Y1,Y2,Z1,Z2,A1,A2 
FOR Channel-1 TO 8 

Set (Channel, 1 ) =Channel 
SELECT Commands 

CASE "P" ! Command$="P" indicates we want to view the encoder Positions in inches. 

NameS (Channel , 1) =Name$ (Channel, 1) 4" (pos) " 

Units$ (Channel, 1) -"in" 

CASE "U" ! Command$="U" indicates we want to view the Units in counts per inch. 

Name$ (Channel, 1) =Name$ (Channel, 1) 6" (cpi) " 

UnitsS (Channel, 1) ="cnt" 

CASE "R" ! Command$="R" indicates we want to view the number counts per Revolution^ 

NameS (Channel, 1) =Name$ (Channel, 1) 4" (cpr) " 

UnitsS (Channel, 1) ="cnt" 

CASE "V" ! Command$-"V" indicates we want to view the Velocity in revolution per second. 

NameS (Channel , 1) =Name$ (Channel, 1)4" (vel) " 

UnitsS (Channel, 1) =" rev" 

CASE "A" ! Command$-"A" indicates we want to view the Acceleration in revolution per second A 2. 

NameS (Channel, 1) -NameS (Channel, 1) 4" (acc) " 

UnitsS (Channel, 1) -"rev" 

CASE " + " ! Command$="+" indicates we want to view the current + direction limit switches. 

NameS (Channel, 1) =Name$ (Channel, 1) 4" ( +LS) " 

UnitsS (Channel, 1) =" " 

CASE ! Command$="- M indicates we want to view the current - direction limit switches. 

NameS (Channel, 1) =Name$ (Channel, 1) 4" (-LS) " 

UnitsS (Channel, 1) =" " 

CASE "S" ! Command$s="S" indicates we want to view the current motor Stall indication status. 

NameS (Channel, 1) =Name$ (Channel, 1) 4" (STALL)" 

UnitsS (Channel, 1) " 

CASE "H" ! Command$="H" indicates we want to view the current Home limit switches. 

NameS (Channel, 1) =Name$ (Channel, 1) 4" (HS) " 

UnitsS (Channel, 1) =" " 

END SELECT 
NEXT Channel 

! The "Change" subprogram allows the user to see and then change the values of the viewed parameters 
CALL Change ( "VALUES", View ( *) ,Name$(*) ,Image$<*) ,Units$(*) ) 

! The "Set" parameters command is now sent to the TCS8 

SELECT Commands 

CASE "P», »U", "R", "V", "A" 

MAT Set ( * , 2) = View(*,l) 

OUTPUT @Tcs8 USING 23940; "S"4Command$ , Set (* ) 

IMAGE K, 8 (D,":",M6D.4D, ","),/ 

END SELECT 
SUBEND 

SUB Tcs8read (@Tcs8,Mod(*) ,Tun (*) , Tcsl (*) , Tcs2 <*) , Tcs2tunl (*) , Tcs2tun2 (*) ;Tun2mod(M ) 

! Description: 
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! This subprogram reads the current TCS8 positions. The 8 positions are read in TCS 

! coordinates with the units being in inches. Four of the eight positions (XI, Yl, Zl, Al) which 

! are the transmitting side traverse positions are entered into the Tcsl array. The other four 

! positions (X2, Y2, Z2, A2) which are the receiving tide traverse positions are entered into the 

! Tcs2 array. The Tcsl & Tcs2 arrays are converted from TCS to TUNNEL to MODEL coordinates. 

! The current updated positions in the three coordinate systems are printed on the top of the 

! CRT. They are also returned to the main program. 

! Variables: 

! Tcsl (*) TCS8 transmitting side traverse positions (XI , Yl, Zl, Al) in TCS coordinates. 

! Tcs2 ( * ) TCS8 receiving side traverse positions (X2, Y2, Z2, A2) in TCS coordinates. 

! Tun(*) Traverse positions (X,Y, Z) in TUNNEL coordinates. 

! Mod(*) Traverse positions (X,Y,Z) in MODEL coordinates. 

! Tcs2tunl(*) Coordinate system transformation matrix for converting Tcsl(*) to Tun(*). 

! Tcs2tun2 ( * ) Coordinate system transformation matrix for converting Tcs2(*) to Tun(*). 

! Tun2mod(*) Coordinate system transformation matrix for converting Tun(*) to Mod(*). 

OUTPUT @Tcs8 USING "K,/";"VP0" 

ENTER @Tcs8 USING "8 (K) Tcsl (1) , Tcs2 (1) , Tcsl (2) , Tcs2 (2) , Tcsl (3) , Tcs2 (3) , Tcsl (4) , Tcs2 (4) 

MAT Tun- Tcs2tunl *Tcsl 
REDIM Tun (1:3) ,Mod(l:3) 

MAT Mod= Tun2mod*Tun 
REDIM Tun (1:4), Mod (1:4) 

Mod (4 ) =0 
Tun (4 ) =0 

CALL Tcs8print (Mod(*) , Tun (*) , Tcsl (*) , Tcs2 (*) ) 

SUBEND 

SUB Tcs8print (Mod (*) , Tun <*) , Tcsl (*) , Tcs2 (* ) ) 

! Description: 

! This subprogram prints the current updated TCS8 positions at the top of the CRT. The 

! positions are printed in TCS , TUNNEL , and MODEL coordinates. 

I Variables: 

! Tcsl(*) TCS8 transmitting side traverse positions (XI , Yl , Zl, Al) in TCS coordinates. 

! Tcs2(*) TCS8 receiving side traverse positions (X2 , Y2, Z2, A2) in TCS coordinates. 

! Tun(*) Traverse positions (X,Y,Z) in TUNNEL coordinates. 

! Mod(*) Traverse positions (X,Y,Z) in MODEL coordinates. 

PRINT CHR$ (128) ; 

PRINT TABXY (50,1) ; " 

PRINT TABXY (50,2) ; " MOD TUN TCS1 TCS2 M 

PRINT TABXY (50, 3) 

PRINT TABXY (50, 4); 

PRINT USING "#, K, 4 (M3D.4D) ,X" ; " X : Mod (1 ), Tun (1) , Tcsl ( 1) , Tcs2 (1) 

PRINT TABXY (50, 5); 

PRINT USING "#,K, 4 (M3D.4D) ,X";" Y : Mod (2 ) , Tun (2) , Tcsl (2) , Tcs2 (2) 

PRINT TABXY (50, 6) ; 

PRINT USING "#,K, 4 (M3D.4D) ,X";" Z : " , Mod ( 3 ) , Tun (3) , Tcsl ( 3) , Tcs2 (3) 

PRINT TABXY (50, 7) ; 

PRINT USING "#,K, 4 (M3D.4D) ,X";" A:", Mod (4 > , Tun (4) , Tcsl (4) , Tcs2 (4) 

PRINT TABXY (50, 8) ; " " 

SUBEND 

SUB Tcs8move ( @Tcs8, Mod ( *) ,Tun(*) , Tcsl (*) , Tcs2 ( *) , Mod2tun ( *) , Tun2tcsl ( * ) , Tun2t cs2 ( *) , Side$, Coor$, Mode$, 

K, Movement) 

! Description: 

! This subprogram allows for the movement of the probe volume and collecting optics in one of 

! three coordinate systems. The three coordinate systems implemented are the TSC, TUNNEL and 

! MODEL coordinate systems. Two movements modes are available. The first movement mode makes 

! moves relative to the current position. The second movement mode makes moves to an absolute 

! fixed position. Both the transmitting side and receiving side traverses can be moved in tandem 

! or separately. 

! Variables: 

! Tcsl ( * ) TCS8 transmitting side traverse positions (XI , Yl , Zl, Al) in TCS coordinates. 

! Tcs2(*) TCS8 receiving side traverse positions (X2 , Y2, Z2, A2) in TCS coordinates. 

! Tun(*) Traverse positions (X,Y,Z) in TUNNEL coordinates. 

! Mod (* ) Traverse positions (X,Y,Z) in MODEL coordinates. 

! Mod2tun(*) Coordinate system transformation matrix for converting Tcsl(*> to Tun(*). 

! Tun2tcsl(*) Coordinate system transformation matrix for converting Tcs2(*) to Tun(*). 

! Tun2tcs2 ( * ) Coordinate system transformation matrix for converting Tcs2(*) to Tun(*). 

! Side$ Indicates which sides are to be moved: 

! Tx : Transmitting side only. 

! Rx : Receiving side only. 

! Tx & Rx : Both sides together. 

! Coor$ Indicates which coordinate system the movement is to be made: 

! MODEL : MODEL coordinates. 

! TUNNEL : TUNNEL coordinates. 

! TCS : TCS coordinates. 

! Mode$ Indicates which movement mode is to be completed: 

! RELATIVE: Movements are relative to current positions. 

! ABSOLUTE: Movements are to absolute positions. 

! K Indicates which axis of the four axes is to be moved. 

! Movement Indicates the desired movement for the selected axis. 

! I(*) Array of viewed TCS8 "Initialized" parameters. 

! C ( *) Array of viewed TCS8 "Currents On" parameters. 


89 


24780 

24790 

24800 

24810 

24820 

24830 

24840 

24850 

24860 

24870 

24880 

24890 

24900 

24910 

24920 

24930 

24940 

24950 

24960 

24970 

24980 

24990 

25000 

25010 

25020 

25030 

25040 

25050 

25060 

25070 

25080 

25090 

25100 

25110 

25120 

25130 

25140 

25150 

25160 

25170 

25180 

25190 

25200 

25210 

25220 

25230 

25240 

25250 

25260 

25270 

25280 

25290 

25300 

25310 

25320 

25330 

25340 Ctm: 

25350 Refract: 

25360 

25370 

25380 

25390 

25400 

25410 

25420 

25430 

25440 

25450 

25460 

25470 

25480 

25490 

25500 

25510 

25520 

25530 

25540 

25550 Ctm_ldv: 

25560 

25570 


OPTION BASE 1 
DIM L$ [100] 

REAL Move (8,2),I(8),C(8) 

! If all of the channels have not yet been Initialized, then do so now. 

OUTPUT @Tcs8 USING "K,/";"VI0" 

ENTER 0Tcs8 USING "8(K)";I(*> 

IF SUM (I) 08 THEN OUTPUT @Tcs8 USING "K, /"; "SI0" 

! If all of the channels do not have their currents turned on, then do so now. 

OUTPUT STcs8 USING "K,/";"VC0" 

ENTER GTcs8 USING “8 (K) "; C (*) 

IF SUM (C) 08 THEN OUTPUT GTcs8 USING "K, /"; "SCO: 1 , M 

! If the movement mode is to be RELATIVE, then clear all of the previously read positions. 
IF Mode$=" RELATIVE" THEN 
MAT Mod= (0) 

MAT Tun= (0) 

MAT Tcsl- (0) 

MAT Tcs2= (0) 

END IF 

! Set the new Tcsl(*) and Tcs2(*) position arrays. 

SELECT Coor$ 

CASE "MODEL" 

Mod (K) =Movement 
REDIM Tun ( 1 : 3 ) , Mod (1:3) 

MAT Tun= Mod2tun*Mod 
REDIM Tun (1:4) , Mod (1:4) 

IF POS (Side$, "Tx") THEN MAT Tcsl= Tun2tcsl*Tun 
IF POS (Side$, "Rx") THEN MAT Tcs2= Tun2tcs2*Tun 
CASE "TUNNEL" 

Tun (K) =Movement 


IF POS(Side$, "Tx") THEN MAT Tcsl- Tun2tcsl*Tun 
IF POS (Side$, "Rx") THEN MAT Tcs2= Tun2tcs2*Tun 
CASE "TCS" 

IF P0S(Side$, "Tx") THEN Tcsl (K) ^Movement 
IF P0S(Side$, "Rx") THEN Tcs2 ( K) ^Movement 
END SELECT 

! File the move array. 

FOR Channel=l TO 8 

Move (Channel, 1) =Channel 
NEXT Channel 


Moved, 2) =Tcsl (1) 

Move (2,2) =Tcs2 (1) 

Moved, 2) =Tcsl (2) 

Move (4,2) =Tcs2 (2) 

Move (5,2) =Tcs 1(3) 

Move (6, 2) =Tcs2 (3) 

Move (7, 2) =Tcsl (4) 

Move (8,2) =Tcs2 (4) 

! Initiate the start of the move. 

IF Mode$= "ABSOLUTE" THEN OUTPUT @Tcs8 USING 2 5280; "MA", Move ( * ) 

IF Mode$=" RELATIVE" THEN OUTPUT @Tcs8 USING 25280? "MR" , Move (* ) 

IMAGE K, 8 (D, " : S2D. 5D, " , " ) , / 

! The TCS8 will return the new updated positions only after the move is complete. 

ENTER @Tcs8 USING "8 (K) "?Tcsl (1) , Tcs2 (1) , Tcsl (2) , Tcs2 (2) , Tcsl (3) , Tcs2 (3) , Tcsl (4) , Tcs2 (4) 
! Turn off the motor drive currents. 

OUTPUT @Tcs8 USING "K, /" ? "SCO : 0, " 

SUBEND 
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SUB Refract (REAL Index (*), Thetal (*), Theta3 (*) ) 

Description: 

This subprogram uses the laser beam angles outside the tunnel to compute the angles inside 
the water tunnel. This requires the knowledge of the indexes of refraction for the various 
! media that the beams go through. The Mediums are air, glass, and water. 

! Variables: 

! Index (*) Array of indexes of refraction. 

! Index (1) : Index of refraction for Air. 

Index(2) : Index of refraction for Glass. 

Index<3) : Index of refraction for Water. 

Thetal (*) Laser beam angles outside the water tunnel. 

Theta3(*) Laser beam angles inside the water tunnel. 

OPTION BASE 1 

! Correct Theta for angles in water. 

MAT Theta3= Thetal 
IF Index (1) <> Index (3) THEN 

Theta3 (2, 1) =ASN ( Index (1) /Index (3) *SIN (Thetal (2, 1) ) ) 

Theta 3 (2,2) =ASN (Index (1) /Index (3) *SIN (Thetal (2, 2) ) ) +90 
END IF 
SUBEND 

SUB Ctm^ldv (Theta (*) ,Tun21dv(*) ,Ldv2tun(*) ) 

! Description: 

I This subprogram computes directly the TUNNEL to LASER coordinate system transformation 
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! matrix "Tun21dv (*) *' . However, the desired coordinate system transformation matrix M Ldv2tun M is 

! required. It is the matrix inverse of "Tun21dv , ‘. 

! Variables: 

! Theta(*) Laser beam angles inside the water tunnel. 

t Tun21dv (*) Coordinate system transformation matrix for converting from TUNNEL to LASER. 

! Ldv2tun (*) Coordinate system transformation matrix for converting from LASER to TUNNEL. 

OPTION BASE 1 

! Tun21dv converts TUNNEL coordinates to LASER coordinates. 

Tun21dv (1, 1) =COS (Theta (1, 1) ) 

Tun2 1 dv ( 1 , 2 ) =C0S ( Thet a (1, 2) ) 

T un2 1 d v < 1 , 3 ) = COS ( Thet a ( 1 , 3 > ) 

Tun21dv (2, 1) =C0S (Theta (2,1)) 

Tun21dv (2, 2) =C0S (Theta (2, 2) ) 

Tun21dv (2 , 3) ^COS (Theta (2, 3) ) 

Tun21dv (3,1) -COS (Theta (3, 1) ) 

Tun21dv (3, 2) =COS (Theta (3, 2) ) 

Tun21dv (3,3) *COS (Theta (3, 3) ) 

! Ldv2tun converts LASER coordinates to TUNNEL coordinates. 

MAT Ldv2tun= INV(Tun21dv) 

SUBEND 

SUB Ctm_mod (Alpha ( * ) , Mod2tun ( * ) , Tun2mod ( * ) ) 

! Description: 

! This subprogram computes directly the MODEL to TUNNEL coordinate system transformation 

! matrix "Mod2tun (*) M . However, the desired coordinate system transformation matrix M Tun2mod" is 

I required. It is the matrix inverse of "Mod2tun". 

! Variables: 

! Alpha (*) Angles of attack, yaw, and roll. 

i T1 (*) Partial coordinate system transformation matrix for converting from MODEL to 

! TUNNEL coordinates. Takes into account a model at angle of attack. 

i t 2 ( * ) Partial coordinate system transformation matrix for converting from MODEL to 

i TUNNEL coordinates. Takes into account a model at angle of yaw. 

I T3 ( * ) Partial coordinate system transformation matrix for converting from MODEL to 

! TUNNEL coordinates. Takes into account a model at angle of roll. 

J Mod2tun(*) Coordinate system transformation matrix for converting from MODEL to TUNNEL. 

! Tun2mod(*) Coordinate system transformation matrix for converting from TUNNEL to MODEL. 

OPTION BASE 1 

REAL Tl (3,3) , T2 (3,3) , T3 (3, 3) , Temp (3,3) , 

! Define 1st coordinate transformation matrix for Mod2tun. 

J Rotation in the x-y plane about the z-axis. 

! Used when model is at an angle of attack. 

T1(1,1)=C0S (Alpha (1) ) 

Tl (1, 2) *SIN (Alpha (1) ) 

Tl (1, 3) =0 

Tl (2, 1) --SIN (Alpha (1) ) 

Tl (2, 2) =C0$ (Alpha (1) i 
Tl (2, 3} =0 
Tl (3, 1) =0 
Tl (3, 2) =0 
Tl (3, 3) =1 

! Define 2nd coordinate transformation matrix for Mod2tun. 

! Rotation in the x-z plane about the y-axis. 

! Used when model is at an angle of yaw. 

T2 (1, 1) =C0S (-Alpha (2) ) 

T2 (1, 2) -0 

T2(l,3)= s -SIN (-Alpha (2) ) 

T2 (2, 1) =0 
T2 (2, 2) =1 
T2(2, 3) =0 

T2 (3, 1) =SIN (-Alpha (2) ) 

T2 ( 3, 2) -0 

T2 (3, 3) =COS (-Alpha (2) ) 

! Define 3rd coordinate transformation matrix for Mod2tun. 

! Rotation in the y-z plane about the x-axis. 

! Used when model is at an angle of roll. 

T3 (1, 1) =1 
T3 (1, 2) =0 
T3 (1, 3) =0 
T3 (2, 1) =0 

T3 (2, 2) =C0S (-Alpha (3) ) 

T3 (2, 3) =SIN (-Alpha (3) ) 

T3 (3, 1) =0 

T3 (3, 2) --SIN (-Alpha (3) ) 

T3 (3, 3) =COS (“Alpha (3) ) 

! Mod2tun converts MODEL coordinates to TUNNEL coordinates. 

MAT Temp= T2*T1 
MAT Mod2t un= T3*Temp 

! Tun2mod converts TUNNEL coordinates to MODEL coordinates. 

MAT Tun2mod= INV(Mod2tun) 

SUBEND 

SUB Ctm_tcs (Tcs2tunl (*) ,Tcs2tun2(*) ,Tun2tcsl(*) ,Tun2tcs2(*) , Fs, Fr, Bs, Br, Index <*) , Ts, Tr, Ta) 
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! Description: • 

t This subprogram computes the TUNNEL to TCS coordinate system transformation matrices 

! "Tun2tcsl (*) M and "Tun2tcs2 {*) The coordinate system transformation matrices "Tcs2tunl“ and 

! "Tcs2tunl M are the matrix inverses of "Tun2tcsl ( * ) “ and "Tun2t cs2 ( *) " respectively. 

! Variables: 

! Tcs2tunl(*) Sending side coordinate transformation matrix converting Tcs(*) to Tun(*) . 

! Tun2tcsl (*) Sending side coordinate transformation matrix converting Tun(*) to Tcs(*) . 

! Tcs2tun2 (*) Receiving side coordinate transformation matrix converting Tcs<*) to Tun(*). 

! Tun2tcs2 (*) Receiving side coordinate transformation matrix converting Tun<*) to Tcs(*). 

! Fs Focal length for sending side onaxis and offaxis lenses. 

! Fr Focal length for receiving side offaxis lens. 

! Bs Beam spacings for sending side onaxis and offaxis beam pairs. 

I Br Beam spacing for receiving side offaxis. 

! Index(*) Array of indexes of refraction for air, glass, and water. 

! Ts Angle of offaxis sending side beam pair. 

! Tr Angle of offaxis receiving side beam pair. 

! Ta Sending side offaxis auxiliary traverse angle. Returned to main program. 

! Xs_on,Ys_on Starting coordinates of onaxis sending side lens. 

! Xs_of f s, Ys_of fs Starting coordinates of offaxis sending side lens. 

! Xs_off r, Ys_offr Starting coordinates of offaxis receiving side lens. 

! Xc,Yc The common point in air of two beam path equations, 

! Ba , Bb The Y intercepts of, two beam path equations. 

! Theta (*) Array of angles in which each beam contacts the window. 

! X ( * ) Array of X coordinates for the points in which each beam contacts the window. 

! Yposition The Y coordinate of the point where all beams cross in the water. 

! Y1 , X2, Y2, X3, Y3 Temporary variables to hold the results of the first call to Findstart. 

! Thickness The thickness of the window. 

! Beam Subscript used while determining the X{*) array. 

OPTION BASE 1 

REAL Xs_on, Ys_on, Xs_offs, Ys__offs, Xs_of fr, Ys_offr, Xc, Yc, Ba, Bb, Theta (6) , X (6) 

REAL Yposition, Y1,X2,Y2,X3,Y3, Thickness 

INTEGER Beam 

Thickness^l .25 

Yposi tion~0 

GOSUB Findstart 

Yl=Ys_on 

X2=Xs_offs 

Y2=Ys_of f s 

X3=Xs_offr 

Y3=Ys_of f r 

Yposit ion=l 

GOSUB Findstart 

MAT Tun2tcsl- IDN 

MAT Tun2tcs2= IDN 

Tun2t csl (2,2) =-Ys_on+Yl 

Tun2tcsl (4,2) »-SQRT < <Xs_of fs-X2) A 2+(Ys_offs- (Ys_on-Yl+Y2) ) A 2) 

Tun2tcsl ( 4, 4) =0 

Ta=ATN ( (Xs_of fs-X2) / <Ys_offs- (Ys_on-Yl+Y2 ) j ) 

Tun2tcs2(l,2) =Xs_offr-X3 
Tun2tcs2 (2, 2) =-Ys_of f r+Y3 
Tun2tcs2 ( 4, 4 ) =0 
MAT Tcs2tunl= INV (Tun2tcsl) 

MAT Tcs2tun2= INV (Tun2tcs2) 

Tcs2tunl {4,2} =0 
Tcs2tun2 (4,2) =0 
SUBEXIT 

! This subroutine finds the starting coordinates for the onaxis (Xs_on, Ys_on) and offaxis (Xs_offs, 

! Ys_offs) sending side lenses and the offaxis (Xs_of fr, Ys_of fr) receiving side lens given the point 
! at which all beams cross in the tunnel. The crossing point is given to be (Opposition) . The 
! method in which the starting coordinates are found involves solving simultaneously the equations for 
! the path of each beam pair in air yielding the common point (Xc,Yc) . Given the focal length of the 
I lens, the starting coordinate can be calculated. The equation for each beam path in air is obtained 
! by determining the angle and the point a beam contacts the window. 

! 

! These six equations find the six angles. 

Theta (1) =ATN ( Bs/ (2*Fs) ) 

Theta (2) =-ATN (Bs/ (2*Fs) ) 

Theta (3) =Ts+ATN (Bs/ (2*Fs) ) 

Theta (4) =Ts-ATN (Bs/ (2*Fs) ) 

Theta (5) =Tr+ATN (Br/ <2*Fr) ) 

Theta (6) -Tr-ATN (Br/ (2*Fr) ) 

! This equation finds the X coordinate of the six points. The Y coordinate is equal to -Thickness. 

FOR Beam=l TO 6 

X (Beam) =-Yposition*TAN (ASN (Index (1) /Index (3) *SIN (Theta (Beam) ) ) )- 
Thickness*TAN ( ASN ( Index ( 1 ) /Index (2 ) * SIN (Theta ( Beam) ) ) ) 

NEXT Beam 

! Determine the Y intercepts for the onaxis beam paths. 

Ba=-Thickness-X(l) /TAN (Theta (1) ) 

Bb=-Thickness-X (2) /TAN (Theta (2) ) 

! Solve for the common point. 
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27170 

27180 

27190 

27200 

27210 

27220 

27230 

27240 

27250 

27260 

27270 

27280 

27290 

27300 

27310 

27320 

27330 

27340 

27350 

27360 

27370 

27380 

27390 

27400 

27410 

27420 

27430 

27440 

27450 

27460 

27470 

27480 

27490 

27500 

27510 

27520 

27530 

27540 

27550 

27560 

27570 

27580 

27590 

27600 

27610 

27620 

27630 

27640 

27650 

27660 

27670 

27680 

27690 

27700 

27710 

27720 

27730 

27740 

27750 

27760 

27770 

27780 

27790 

27800 

27810 

27820 

27830 

27840 

27.850 

27860 

27870 

27880 

27890 

27900 

27910 

27920 

27930 

27940 

27950 


Xc= (Bb-Ba) / (1/TAN (Theta (1) ) -1/TAN (Theta (2) ) ) 

Yc=Xc/TAN (Theta (2) ) +Bb 

! Calculate the onaxis lens starting coordinate using the focal length and the onaxis angle (=0 deg). 
Xs_on~Xc-Fs*SIN (0) 

Ys_on=Yc-Fs*COS (0) 

! Determine the Y intercepts for the offaxis sending side beam paths. 

Ba=-Thickness-X (3) /TAN (Theta (3) ) 

Bb=-Thickness-X ( 4 ) /TAN (Theta (4) ) 

! Solve for the common point. 

Xc= ( Bb-Ba ) / ( 1 /TAN (The t a ( 3 ) > -1 /TAN (Theta ( 4 ) ) ) 

Yc=Xc/TAN (Theta (4) ) +Bb 

! With the focal length and the offaxis angle calculate the starting coordinate 
! of the offaxis sending side lens. 

Xsof fs=Xc-Fs*SIN (Ts) 

Ys_of fs~Yc*“Fs*COS (Ts) 

! Determine the Y intercepts for the offaxis receiving side beam paths. 

Ba«-Thickness-X (5) /TAN (Theta ( 5) ) 

Bb=-Thickness-X (6) /TAN (Theta ( 6) ) 

! Solve for the common point. 

Xc= (Bb-Ba) / (1/TAN (Theta (5) ) -1/TAN (Theta (6) ) ) 

Yc=Xc/TAN (Theta ( 6) ) + Bb 

1 With the focal length and the offaxis angle calculate the starting coordinate 
! of the offaxis receiving side lens. 

Xs_of f r=Xc-Fr *SIN (Tr) 

Ys_of fr=Yc-Fr *C0S (Tr) 

RETURN 

SUBEND 


Graph: 

Crt init: 


!!!!!!!!!!!!!!!!!!!!!!!!!! !!! !!!!!!! i !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 1 !!!!!!!!!!!!!! 1 !!!!! 1 !!!! ! 
SUB Crt_init 

Description: 

This subprogram initializes the CRT as the plotting device and clears both the alpha 
numerics and graphics part of the CRT. 


Select the CRT as the plotting device. 

Select black for area fills. 

Select white for line drawing and labeling. 
Select the CRT as the printing device. 

Send ERROR and DISP message to CRT. 

Hide the special function key labels for fl..f8. 


Read_symbols : 


PLOTTER IS CRT, "INTERNAL" 

AREA PEN 0 
PEN 1 

PRINTER IS CRT 
PRINTALL IS CRT 
KEY LABELS OFF 
SUBEND 

SUB Read_symbols (Symbols (*) ) 

Description: 

This subprogram defines 5 geometric symbols to be used with the SYMBOL statement. The 
symbols provided are as . follows: Square, Octagon, Diamond, and Triangles (upwards & downwards 

pointing triangles) . All of the symbols have a dot added to their center. 

Variables: 

Symbols(*) Array of Symbol arrays. Each symbol arrays contains a distinct geometric symbol. 
Symbol (*) Array of coordinates which when connected produce a distinct geometric symbol. 

Dot ( * ) Array of coordinates which produce a dot. The dot symbol is added to all symbols. 

Noc The number of coordinates in a symbol. 

S Used to index the Symbols array. 

OPTION BASE 1 
REAL Symbol (20, 3) , Dot (2, 3) 

READ Dot ( * ) 

FOR S=1 TO 5 
READ Noc 

REDIM Symbol (Noc, 3) 

READ Symbol (*) 

MAT Symbols (S, 1 :Noc, *) = Symbol 
MAT Symbols (S,Noc+l :Noc+2, *) = Dot 
Symbols (S, 0, 1 ) =Noc+2 


6.5,11.5,-!, 


Dot : 

NEXT 

DATA 

S 

4.5, 7.5, -2, 

4.5, 7. 5,-1 




Square: 

DATA 

5, 

0.5, 3.5, -2, 

8.5, 3. 5,-1, 

8. 5, 11.5, -1, 

0.5, 11. 5,-1, 

0.5, 3. 5,-1 

Octagon: 

DATA 

9, 

0.5, 5. 5, -2, 

2. 5, 11. 5,-1, 

2.5, 3. 5,-1, 
f 0.5, 9. 5,-1, 

6.5, 3.5, -1, 
0.5, 5. 5,-1 

8.5, 5. 5,-1, 

8.5, 9. 5,-1 

Diamond: 

DATA 

5, 

-0.5, 7.5, -2, 

4.5, 2. 5,-1, 

9.5, 7.5, -1, 

4. 5, 12. 5,-1, 

-0.5, 7. 5,-1 

Utriangle: 

DATA 

4, 

0.5, 4.5, -2, 

8.5, 4. 5,-1, 

4.5, 13. 5,-1, 

0.5, 4. 5,-1 


Dtriangle: 

DATA 

SUBEND 

4, 

0.5, 10. 5, -2, 

8.5,10.5,-!, 

4.5, 1.5, -1, 

0.5, 10. 5,-1 



Setup_graph: SUB Setup_graph (Array (*) , Image$ (*) , Paxis, Symbols (*) ) 


! 


Description: 

This subprogram sets up nine empty plots on the CRT screen. Four plots are profile plots 
while the other five plots are histogram plots. The profile and histogram plots provided are 
as follows; Graph# Type Description 

1 Histogram #1 U frequency data in MHz. 

2 Histogram #2 V frequency data in MHz. 

3 Histogram #3 W frequency data in MHz. 

4 Histogram #4 Analog Channel #1 data in volts. 

5 Histogram #5 Analog Channel #2 data in volts. 

6 Profile Plot #1 Velocity Averages versus Traverse Position. 
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27960 

27970 

27980 

27990 

28000 

28010 

28020 

28030 

28040 

28050 

28060 

28070 

28080 

28090 

28100 

28110 

28120 

28130 

28140 

28150 

28160 

28170 
28180 
28190 
28200 
28210 
28220 
28230 
28240 
28250 
28260 
28270 
28280 
28290 
28300 
28310 
28320 
28330 
28340 
28350 
28360 
28370 
28380 
28390 
28400 
28410 
28420 
28430 
28440 
28450 
28460 
28470 
28480 
28490 
28500 
28510 
' 28520 
28530 
28540 
28550 
28560 
28570 
28580 
28590 
28600 
28610 
28620 
28630 
28640 
28650 
28660 
28670 

28680 

28690 

28700 

28710 

28720 

28730 


X axis. 
Y axis. 


Array containing the plot positions and scales. 

String array containing image formats for the axes labeling. 

Array containing the plot's scales. 

Array containing the plot's CRT position. 

Array containing the number of X divisions for the plot's 

Array containing the number of Y divisions for the plot's 

String array containing labels for the X axis. 

String array containing labels for the Y axis. 

String array containing labels for the Plots. 

String array containing image formats for the X axis labeling. 

String array containing image formats for the Y axis labeling. 

String array containing labels for each symbol in a profile plot. 

Array of Symbol arrays. Each symbol arrays contains a distinct geometric symbol. 
Used as an index to the above arrays. Specifies one of nine plots. 

Used an an index to the LegendS array. 


7 Profile Plot #2 Velocity SDVs versus Traverse Position. 

8 Profile Plot #3 Velocity Shear Stresses versus Traverse Position. 

9 Profile Plot #4 Average voltages & SDVs versus Traverse Position. 

Variables: 

! Array (*) 

! Image$ (*) 

! Wndw ( * ) 

! Vwprt (*) 

! Xdiv(*) 

! Ydiv (*) 

! XlabelS (*) 

! YlabelS (*) 

! Title$ (*) 

! XimageS (*) 

! YimageS (*) 

! Legends {*) 

! Symbols {*) 

! G 

! I 

OPTION BASE 1 

COM /Graph/ Wndw (*) ,Vwprt <*), Xdiv (*), Ydiv (*) , XlabelS (*), YlabelS <*) , TitleS (*), XimageS < *) , YimageS (*) , 
Legends (*) 

MAT Wndw= Array ( 61 : 69, * ) 

MAT Vwprt = Array (71:79, *) 

MAT Xdiv (1 :5) - Ar ray < 81 : 85, 1 > 

MAT Xdiv (6:9)= Array ( 81 : 84 , 3) 

MAT Ydiv (1:5) — Array (81:85,2) 

MAT Ydiv (6:9)= Array ( 81 : 84, 4) 

MAT Ximage$= ImageS ( 61 : 69, 1 ) 

MAT Yimage$= ImageS ( 61 : 69, 3) 

FOR G=1 TO 9 

READ G, Xlabel $ (G) 

FOR 1=1 TO SIZE(Legend$,2) 

READ Legends (G, I) 

NEXT I 
SELECT G 
CASE 1 TO 5 

YlabelS (G)= M " 

CASE 6 TO 9 

YlabelS (G)=CHR$ (NUM("X") +Paxis-1) 

END SELECT 

CALL Set_up (G, Symbols {*) ) 

NEXT G 
SUBEXIT 


I 

G, 

X axis Label 


, Symbol #1 

. . . 5 labels 

DATA 

1, 

II II 


II II 

9 9 

II II 

9 

II II II II II II 
/ 9 

DATA 

2, 

II II 


II H 

9 9 

II II 

/ 

II II II H II II 
9 9 

DATA 

3, 

II II 


II II 

9 9 

<1 II 

9 

II II II II II M 
/ / 

DATA 

4, 

II It 


II II 

9 9 

II II 

9 

II II II II M M 
9 9 

DATA 

5, 

II II 


II II 

II II 

9 

II II II II IP II 
9 9 

DATA 

6, 

"Velocities / 

Uedge" 

, "U", 

"V", 

11 W I 0 11 11 11,1 

DATA 

7, 

"RMS / Uedge" 


, "Ul", 

"VI", 

"Wl", "" 

DATA 

8, 

"Shear Stress 

/ Uedge A 2" 

, "Ulvl ", 

"Vlwl" , 

“Wlul", " ", "" 

DATA 

9, 

"Fluorescence : 

r C, Cl (volts) " 

, "C", 

It II 

9 

"Cl","", "" 


SUBEND 

Set_up: SUB Set_up (G, Symbols (*) ) 

Description: 

This subprogram clears and then redraws one of nine empty plots on the CRT screen. 

Variables: 

Wndw(*) Array containing the plot's scales. 

Vwprt (*) Array containing the plot's CRT position. 

Xdiv(*) Array containing the number of X divisions for the plot's X axis. 

Ydiv (* ) Array containing the number of Y divisions for the plot's Y axis. 

XlabelS (*} String array containing labels for the X axis. 

Ylabel$(*) String array containing labels for the Y axis. 

Title$ (*) String array containing labels for the Plots. 

Ximage$(M String array containing image formats for the X axis labeling. 

Yimage$(*) String array containing image formats for the Y axis labeling. 

Legend$(*) String array containing labels for each symbol in a profile plot. 

Symbols (*) Array of Symbol arrays. Each symbol arrays contains a distinct geometric symbol. 

G Used as an index to the above arrays. Specifies one of nine plots. 

OPTION BASE 1 

COM /Graph/ Wndw(*) , Vwprt (*) , Xdiv (*) , Ydiv <*) , XlabelS ( *) , YlabelS (*) , Title$ <*)., XimageS <*), YimageS (*) , 
Legends (*) 

DIM L$ ( 80 ] 

ON ERROR CALL Error 
PLOTTER IS CRT, "INTERNAL" 

! Define the pen numbers for the colors black and white. 

Black=-1 
White=l 


! 
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28740 

28750 

28760 

28770 

28780 

28790 

28800 

28810 

28820 

28830 

28840 

28850 

28860 

28870 

28880 

28890 

28900 

28910 

28920 

28930 

28940 

28950 

28960 

28970 

28980 

28990 

29000 

29010 

29020 

29030 

29040 

29050 

29060 Background: 

29070 

29080 

29090 

29100 

29110 

29120 

29130 

29140 Axes: 

29150 

29160 

29170 

29180 

29190 

29200 

29210 

29220 

29230 

29240 

29250 

29260 

29270 

29280 Grid: 

29290 

29300 

29310 

29320 

29330 

29340 

29350 Plot_area: 

29360 

29370 

29380 

29390 Xlabel: 

29400 

29410 

29420 

29430 

29440 

29450 

29460 

29470 

29480 

29490 Ylabel: 

29500 

29510 

29520 

29530 


CSIZE 100*15/1023 ! Select a character labeling size of 15 pixels high. 

! Define the values for the left, right , bottom, top ends of the horizontal and vertical scales. 
Xmin=Wndw(G,l) 

Xmax=Wndw (G, 2) 

Ymin=Wndw(G, 3) 

Ymax=Wndw (G, 4 ) 

I Define the values for the left , right, bottom, top pixel locations for the plot. 

Xpixl=Vwprt (G, 1) 

Xpix2=Vwprt (G, 2) 

Ypixl^Vwprt (G, 3) 

Ypix2=Vwprt (G, 4) 

! Define the step size between grid lines, axis tick marks, and axis labels. 

Xstep= (Xmax-Xmin) /Xdi v (G) 

Ystep= (Ymax-Ymin) /Ydiv (G) 

! Define the amount of scale X and Y which equals the size of one pixel {picture element) • 
Xpixel- (Xmax-Xmin) / <Xpix2-Xpixl) 

Ypixel= (Ymax-Ymin) / (Ypix2-Ypixl) 

! Clear the plots back ground & plot area and also draw the plots borders, grids, and axes. 
AREA PEN Black 
PEN White 
GOSUB Background 
GOSUB Axes 
! GOSUB Grid 
GOSUB Plot_area 

! Draw the X and Y axis labels. 

CLIP OFF 
GOSUB Ylabel 
GOSUB Xlabel 

! Create a legend to define which symbol is used with which data. 

CALL Legend (G, Symbols {*) ) 

OFF ERROR 
SUBEXIT 

I This subroutine clears the plot’s background. 

VIEWPORT (Xpixl-75) /10.23, <Xpix2+25) /10.23, (Ypixl-33) /10.23, (Ypix2+6) /10.23 

WINDOW -1 .E+9,l.E+9,-l.E+9, l.E+9 

MOVE 0,0 

WINDOW 0,1,0, 1 

MOVE 0,0 

RECTANGLE 1,1, FILL 
RETURN 

! This subroutine draws the plot's X and Y axes. 

VIEWPORT (Xpixl-1) /10.23, (Xpix2+1> /10 .23, (Ypixl-6) /10 .23, (Ypixl-1) /10 .23 
WINDOW Xmin, Xmax, 1,0 
AXES Xstep, 2, Xmin, 0, 1,1,1 

VIEWPORT (Xpixl-1) /10 .23, <Xpix2+l ) /10 . 23, (Ypix2+1) /10 .23, (Ypix2+6) /10 .23 
WINDOW Xmin, Xmax, 0,1 
AXES Xstep, 2, Xmin, 0,1, 1,1 

VIEWPORT (Xpixl-6) /10.23, (Xpixl-1 ) /10 . 23, (Ypixl-1) /10 .23, (Ypix2+1) /10 .23 
WINDOW 1, 0, Ymin, Ymax 
AXES 2, Ystep, 0, Ymin, 1,1,1 

VIEWPORT (Xpix2+1) /10.23, (Xpix2+6) /10 . 23, (Ypixl-1 ) /10 .23, (Ypix2+1) /10 .23 
WINDOW 0, 1, Ymin, Ymax 
AXES 2, Ystep, 0, Ymin, 1,1,1 
RETURN 

! This subroutine draws the plot’s X and Y grid lines. 

VIEWPORT (Xpixl-1) /10. 23, (Xpix2+1 ) /10 . 23, (Ypixl-1) /10 .23, (Ypix2+1) /10 .23 
WINDOW Xmin, Xmax, Ymin, Ymax 
LINE TYPE 4 

GRID Xstep, Ystep, Xmin, Ymin 

LINE TYPE 1 

RETURN 

! This subroutine selects part of the CRT plot area and give it scales for the X and Y axes. 
VIEWPORT Xpixl/10 . 23, Xpix2/10 .23,Ypixl/10.23, Ypix2/10 .23 
WINDOW Xmin, Xmax, Ymin, Ymax 
RETURN 

! This subroutine labels the X axis and also names the X axis. 

LORG 5 

FOR X=Xmi n TO Xmax+Xstep/100 STEP Xstep 
MOVE X, Ymin-12*Ypixel 
OUTPUT L$ USING Ximage$ (G) ;X 
LABEL TRIMS (L$) 

NEXT X 

MOVE (Xmin+Xmax) /2, Ymin-25*Yplxel 
LABEL XlabelS (G) 

RETURN 

! This subroutine labels the Y axis and also names the Y axis. 

LORG 8 
Len=0 

FOR Y=Ymi n TO Ymax+Ystep/100 STEP Ystep 
MOVE Xmin-5*Xpixel, Y 
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29540 

29550 

29560 

29570 

29580 

29590 

29600 

29610 

29620 

29630 Legend: 

29640 

29650 

29660 

29670 

29680 

29690 

29700 

29710 

29720 

29730 

29740 

29750 

29760 

29770 

29780 

29790 

29800 

29810 

29820 

29830 

29840 

29850 

29860 

29870 

29880 

29890 

29900 

29910 

29920 

29930 

29940 

29950 

29960 

29970 

29980 

29990 

30000 

30010 

30020 

30030 

30040 

30050 

30060 

30070 

30080 

30090 

30100 

30110 

30120 

30130 Hlsto: 

30140 Rt_histo: 

30150 

30160 

30170 

30180 

30190 

30200 

30210 

30220 

30230 

30240 

30250 

30260 

30270 

30280 

30290 

30300 

30310 

30320 


OUTPUT L$ USING Yimage$ (G) ; Y 
LABEL TRIM$(L$) 

Len«MAX<Len,LEN(TRIM$(L$) ) ) 

NEXT Y 

MOVE Xmin- (5+7*Len) *Xpixel, (Ymin+Ymax) /2 
LABEL Ylabel$ (G) 

LORG 5 
RETURN 
SUBEND 

SUB Legend (G, Symbols (*) ) 

! Description: 

! This subprogram produces a legend within one of the nine plots on the CRT screen. 

! Variables: 

! Wndw(*) Array containing the plot's scales. 

! Vwprt(*) Array containing the plot's CRT position. 

! Xdiv(M Array containing the number of X divisions for the plot’s X axis. 

! Ydiv(*) Array containing the number of Y divisions for the plot’s Y axis. 

! Xlabel$(*) String array containing labels for the X axis. 

! Ylabel$(*) String array containing labels for the Y axis. 

J Title$(*) String array containing labels for the Plots. 

! Ximage$(*) String array containing image formats for the X axis labeling. 

! Yimage${*) String array containing image formats for the Y axis labeling. 

! Legend$ (* ) String array containing labels for each symbol in a profile plot. 

! Symbols (*) Array of Symbol arrays. Each symbol arrays contains a distinct geometric symbol. 

! Symbol (*) Array of coordinates which when connected produce a distinct geometric symbol. 

! G Used as an index to the above arrays. Specifies one of nine plots. 

! S Used to index the Legend$ array. 

! Noc The number of coordinates in a symbol. 

I Len Length of a Legend$ array element. 

OPTION BASE 1 

COM /Graph/ Wndw( *) , Vwprt (*) , Xdiv (*) , Ydiv (*) , XlabelS ( *) , Ylabel$ (*) ,Title$ (*) , Ximage$ <*) , Yimage$ (*) , 
Legend$ ( * > 

DIM Symbol (20, 3) 

VIEWPORT Vwprt (G, 1) /10. 23, Vwprt (G,2) /10. 23, Vwprt (G, 3) /10. 23, Vwprt (G,4) /10.23 
WINDOW Vwprt <G,1) , Vwprt (G,2) , Vwprt (G, 3} , Vwprt (G, 4) 

! Define the pen numbers for the colors black and white. 

Black=-1 

White=l 

! Define the colors for symbol filling and edge drawing. 

AREA PEN Black 
PEN White 

CSIZE 100*15/1023 ! Select a character labeling size of 15 pixels high. 

LORG 2 

! Calculate the maximum length of all of the symbol labels. 

Len=0 

FOR S=1 TO SIZE (Legends, 2) 

Len=MAX (LEN (Legend$ (G, S) ), Len) 

NEXT S 

! For each symbol put up a sample symbol and its label. 

FOR S=1 TO SIZE (Legend$, 2) 

IF LEN (Legend$ (G, S) ) =0 THEN 30110 
Noc=Symbols (S, 0, 1) 

REDIM Symbol (Noc, 3) 

MAT Symbol = Symbols (S, 1 : Noc, * ) 

MOVE Vwprt (G, 2) -7*Len-23, Vwprt (G, 4) -15*S+5 
SYMBOL Symbol (*) , FILL, EDGE 
MOVE Vwprt (G, 2) -7 *Len-10, Vwprt (G, 4)-15*S+4 
LABEL Legends (G, S) 

NEXT S 
SUBEND 

t i ! i I j m j j it i i j j j j m i i i ; i ! t t j t j j ; t t t j j i i t i tt i i i.m i i 1. 1 j i ! j i r t t i it m | i i j i m j t t i j M j j j j j j m i j i j i ! t j j j i i i ! j ! i i 
SUB Rt_histo ( @Lvdas, Symbols (*) , Repeat) 

! Description: 

! This subprogram plots real time histograms within five of the nine plots on the CRT screen. 

! The histogram data is acquired from the LVDAS over a specified acquisition time. 

! Variables Defined in Main Program: 

! Wndw (*) , Vwprt (*) , Xdiv (*), Ydiv (* ) ,Xlabel$(*> ,Ylabel$ (*) , Title$ ( * ) , Ximage$ ( *) , Yimage$ ( *) , Legend$ ( *) 

! Local Variables: 

! Histo ( * ) Array of bin numbers, old histogram bin heights, and new histogram bin heights. 

! Nbins Number of bins in the Histo<*) . 

! Bin 2*Bin is the bin width of individual histogram vertical bars. 

! Min Minimum value for histogram. Left side of histogram scale. 

! Max Maximum value for histogram, right side of histogram scale. 

! FI Upper 16bits of integerized Min. 

t F2 Lower 16bits of integerized Min. 

! A1 Upper 16bits of integerized histogram acquisition time. 

! A2 Lower 16bits of integerized histogram acquisition time. 

! Nnew Number of samples in the most up to date histogram. 

! Nold Number of samples in the previous histogram. 

! N ( * ) Number of samples for each histogram of the five separate channels. 
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! Channel Used to select the LVDAS channel that will be sampled for a histogram. 

! Kw Converts Hz to MHz or raw data to volts. 

! Ww Window width of each vertical histogram bar. 

! Old Histogram height of previous histogram at a particular bin. 

! New Histogram height of current histogram at a particular bin. 

t XI Horizontal position of histogram rectangle. 

J Y1 Vertical position of histogram rectangle. 

! X2 Horizontal width of histogram rectangle. 

! Y2 Vertical width of histogram rectangle. 

! I Used as an index to the Histo (*). Specifies one of Nbins bins. 

! G Used as an index to the graphics arrays. Specifies one of nine plots. 

OPTION BASE 1 

COM /Graph/ Wndw(*) ,Vwprt <*) , Xdiv (*) , Ydiv (*) , Xlabel$ (*) , Ylabel$ (*) ,Title$ <*) , Ximage$ <*) ,Yimage$ (*) , 


Legend$ (*) 

INTEGER Histo (1000, 3) , Nbins, FI, F2, Al, A2 
REAL Nnew, Nold, N ( 5) 

! Clear all of the histogram data within the LVDAS. 

OUTPUT QLvdas USING "AA'V'CA” 

! Draw new plots for the five histograms. 

FOR Channel=l TO 5 

CALL Set_up (Channel, Symbols (*) ) 

NEXT Channel 

! Calculate the acquisition time. 0.1*10000000 will give an acquisition of 0,1 seconds. 
CALL Convert2words (. 1*10000000, A1,A2) 

! Enable the keyboard to terminate histogram plotting. 

ON KBD GOSUB Hdone 
REPEAT 


FOR Channel=l TO 5 
G=Channel 


SELECT Channel 
CASE 1,2,3 

Kw=1000000 
Min=Kw*Wndw (G, 1) 

Max=Kw*Wndw (G, 2) 

Bin-INT (LGT ( (Max-Min) /100 ) /LGT (2) )+l 
Ww=2 A Bin 

CALL Convert2words (Min, FI, F2) 

CASE 4,5 

Kw=32768/5 
Min=Kw*Wndw(G, 1> 

Max=Kw*Wndw (G, 2) 

Bin=INT (LGT ( (Max-Min) /100) /LGT (2) ) +1 
Ww=2 A Bin 

CALL Convert2words (Min, FI, F2) 


! Channels 1,2,3 are for LDV frequency data. 

! Converts Hz to MHz. 

! Minimum frequency for left histogram scale. 

! Maximum frequency for right histogram scale. 

! 2 A Bin is the window width of each vertical bars. 
! Window width of each vertical histogram bar. 

! Channels 4,5 are for analog voltage data. 

! Converts raw data to volts. 

! Minimum voltage for left histogram scale. 

! Maximum voltage for right histogram scale. 

! 2 A Bin is the window width of each vertical bars. 
! Window width of each vertical histogram bar. 


CASE ELSE 

GOTO 31100 
END SELECT 

l Tell the LVDAS to Take a Histogram. 

OUTPUT QLvdas USING "AA, 6 (W) ",* "TH", FI , F2, Bin, A1 , A2, Channel 
! Enter number of bins in the histogram. 

ENTER GLvdas USING W";Nbins 

! Redimension the Histo(*} and the enter the histogram data. 

IF Nbins>0 THEN 

REDIM Histo (Nbins, 3) 

ENTER QLvdas USING "#,W"; Histo(*) 

END IF 

! Enter the number of samples for the previous and current histogram. 

ENTER QLvdas USING »#, W"; Nnew, Nold 

! Scale part of the CRT for the histogram plotting. 

VIEWPORT Vwprt (G, 1) /10.23,Vwprt (G,2) /10.23,Vwprt (G,3) /10.23,Vwprt (G,4) /10.23 
WINDOW Kw*Wndw (G, 1) , Kw*Wndw (G, 2) , Wndw (G, 3) , Wndw (G, 4) 

Xpixel=Kw* (Wndw (Channel, 2} -Wndw (Channel, 1 ) ) / (Vwprt (Channel, 2) -Vwprt (Channel, 1) ) 


N1=N (Channel) 

N2=N (Channel) -Nold+Nnew 
N (Channel) =N(Channel) -Nold+Nnew 
FOR 1=1 TO Nbins 


01d=MIN (Histo ( 1 , 3 ) ,Wndw (Channel, 4) ) 
New=MIN (Histo (1,2) , Wndw (Channel, 4) ) 


Positive pens will plot while negative histograms erase. 
Calculate histogram bar horizontal position. 

Calculate histogram bar horizontal width. 

Calculate histogram bar vertical position. 

Calculate histogram bar vertical width. 

If XKXmin then set Xl=Xmin 
If Xl>Xmax then set Xl=Xmax 


AREA PEN SGN (New-Old) 

Xl=Histo (1,1) *Ww+Min 
X2=Ww 
Yl=01d 

Y2=New-01d ! 

IF Xl<Kw*Wndw (G, 1 ) THEN Xl=Kw*Wndw (G, 1) 

IF Xl>Kw*Wndw (G, 2) -X2 THEN Xl=Kw*Wndw (G, 2) -X2 
MOVE XI, Y1 

RECTANGLE X2-Xpixel, Y2, FILL . ! Draw the rectangle representing one bar of the bar graph. 

NEXT I 

NEXT Channel 

UNTIL KBD$<>" " OR NOT Repeat ! Quit if any key on the keyboard has been pressed. 


! 
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SUBEXIT 


Done=l 

RETURN 

SUBEND 

!!!! n !!!!!!!!!!!!!!!!!!!! I I !!!!£!!!!!! M !!!!!!!!!!!!! t !!!!!!!!!! t !!!! 1 !!!!!!!!!! M !!! M !!!! M !! t !!!!!!!! ! 
SUB Pt_histo (Symbols (*), Run, File, Pos, INTEGER Nsam) 


! 


! 


Description: 

This subprogram plots post time histograms within five of the nine plots on the CRT screen. 

The histogram data is acquired from the LVDAS over a specified acquisition time. 

Variables Defined in Main Program: 

Wndw (*) ,Vwprt (*> ,Xdiv(*> ,Ydiv(*) ,Xlabel$(*) , Ylabel$ (*) , Title$ (*) ,Ximage$(*) , Yimage$ (*) ,Legend$ <*) 
Ui(*) , Vi (*) ,Wi(*) ,Ai(*) ,Bi(*) 

Local Variables: 

Histo {*) Array of histogram bin heights indexed by bin number. 

Array of instantaneous U,V, W velocity or A, B voltage data. 

Number of samples acquired. 

Left side of histogram scale, 
right side of histogram scale, 
histogram bar. 
arrays. 

Used as an index to the Histo (*). Specifies one of 100 bins. 

Horizontal length of one picture on the CRT in scale units. 

Selects one of the 5 channels of Ui (*) , Vi (*) , Wi (*) , Ai (*) , Bi (*) data. 

Used as an index to the graphics arrays. Specifies one of nine plots. 


Data (*) 

Nsam 

Xmin 

Xmax 

Xwin 

K 

L 

Xpixel 

Channel 

G 


Minimum value for histogram. 
Maximum value for histogram. 
Window width of each vertical 
Used as an index to the above 


OPTION BASE 1 

COM /Data/ INTEGER Raw ( *) , Val id < * ) , REAL Table <*) , Ui <*) , Vi <*> , Wi <*) ,Ai (*) , Bi (*), Ii <*> , Ci <*) 

COM /Graph/ Wndw ( *) , Vwprt <*) , Xdiv (*) , Ydiv (*> , XlabelS (*) , Ylabel$ (*) f Title$ <*) , Ximage$ ( *) , Yimage$ (*> , 
Legend$ (*) 

INTEGER Histo (0:100) 

REAL Data (1000) 

REDIM Data (Nsam) 

FOR Channel=l TO 5 

! Fill the data array with Ui (*) , Vi (*) , Wi (*) , Ai (*) , or Bi(*> depending on Channel. 

G-Channel 

IF Channel=l THEN MAT Data= Ui 
IF Channels THEN MAT Data= Vi 
IF Channels THEN MAT Data= Wi 


IF Channels THEN MAT Data= Ai 
IF Channels THEN MAT Data= Bi 
! Draw a new empty histogram plot. 

CALL Set__up (Channel, Symbols (*) ) 

Xmin=Wndw (Channel , 1) 

Xmax-Wndw (Channel , 2) 

Xwin= (Xmax-Xmin) /100 
! Sort the data into a histogram. 

MAT Data= Data- (Xmin) 

MAT Data= Data/ ( (Xmax-Xmin) /100) 

MAT Hlsto= (0) 

FOR K=1 TO Nsam 

L=MAX (MIN (Data (K) , 100) , 0) 

Histo (L) =Histo (L) +1 
NEXT K 

! Scale part of the CRT for histogram plotting. 

VIEWPORT Vwprt (G, 1) /10 . 23, Vwprt (G, 2) /10.23, Vwprt (G,3> /10.23 f Vwprt (G,4)/10.23 
WINDOW 0, 100 / Wndw(G # 3) # Wndw(G, 4) 

Xpixel= (100-0) / (Vwprt (Channel, 2) -Vwprt (Channel, 1) ) 

! Draw the histogram. 

FOR K=0 TO 100 

IF Histo ( K) THEN 
MOVE K- . 5 , 0 

AREA PEN SGN(l) ! Positive pens will plot while negative histograms erase. 
RECTANGLE 1-Xpixel , Histo ( K) , FILL 
END IF 
NEXT K 
NEXT Channel 
SUBEXIT 
SUBEND 
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